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A FONCTION VECTORIELLE 



SES APPLICATIONS A LA PHYSIQUE 
PAR U. B. ÉLIE 

INTRODUCTION 

La théorie de la fonction potentielle présente un exemple de 

ce fait, qu'à l'aide d'un petit nombre de formules et par 

l'emploi d'opérations toujours les mêmes, des théories très 

diverses peuvent être exposées d'une façon uniforme. Un 

nombre restreint de théorèmes abstraits peuvent alors recevoir 

les interprétations concrètes les plus importantes et les plus 

variées. 

> Nous nous proposons, dans ce travail, d'exposer une théorie 

ÎJJ permettant de simplifier l'étude des phénomènes physiques oii 

Centre en première ligne la notion de direction dans l'espace, 

'^la notion de vecteur. 

^5 La fonction vectorielle paraît appelée à jouer en physique 
lun rôle analogue à celui de la fonction potentielle. Mais il est 



[■e, pour que son emploi devienne général, que les 
transformations qui s'y présentent à chaque pas, puissent être 
effectuées à l'aide d'un système simple d'opérations; de là la 
nécessité de l'introduction de notations et de symboles, dont 
les propriétés sont exposées dans notre première partie. 

L'étude des propriétés des fonctions vectorielles dans le cas 
des axes rectangulaires est faite longuement dans la première 
partie. Nous n'avons pas cru devoir abréger cette exposition, 
et cela pour deux motifs : d'abord, la plupart des formules 
auxquelles nous arrivons seront utilisées par la suite; en 



431485 



lauyGoogle 



3 B. ÉL1E. 

second lieu, nous estimons que des notations usuelles ne 
doivent pas être présentées sous une forme par trop concise. 

Le deuxième chapitre de cette première partie a été réservé 
aux coordonnées obliques. Jusqu'ici, leur emploi a été très 
restreint, trop restreint selon nous. Cela est dû probablement 
à ce que les milieux possédant des axes obliques ont été relati- 
vement peu étudiés, et aussi à ce que les formules relatives à 
ce cas sont naturellement de forme un peu plus compliquée. 
Déjà en 1853, Rankine avait appelé l'attention des physiciens 
sur les avantages que présente l'emploi de ces coordonnées. 
Dans un mémoire sur les constantes d'élasticité publié en 
-1886 dans les Mémoires de la Société des Sciences physiques 
et naturelles de Bordeaux, nous avons développé cette idée 
de Rankine. Cela nous permet ici d'être plus bref et d'arriver 
plus rapidement aux applications. 

Dans un tableau donné à la fin du § 2 du chapitre \^'(["' partie), 
nous avons résumé les principales applications à la physique 
de la notion lîe fonction vectorielle. Nous ne pouvions songer 
à en développer les différentes parties. Nous avons dû nous 
restreindre, et dans les quatre chapitres de la seconde partie, 
nous avons montré l'emploi des notations et fornmies établies 
précédemment à quelques problèmes particuliers. 

En fait, c'est l'étude de ces questions ou des questions de 
même nature qui nous a conduit à considérer comme néces- 
saire l'introduction en physique de l'idée de fonction vectorielle, 
et à essayer de la présenter sous une forme systématique. 

Dans un premier chapitre, nous considérons la déformation 
des milieux et nous nous occupons spécialement de la réduc- 
tion du nombre des constantes d'élasticité. 

Le second chapitre traite de la détermination en grandeur 
et en direction de l'électricité dégagée par un cristal, lorsqu'il 
est soumis à un système quelconque de pressions. Les résultats 
auxquels nous arrivons se prêtent à des vérifications expéri- 
mentales faciles à réaliser. 

Il n'a été question jusqu'ici, dans les applications précé- 
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FONCTION VECTOBIELLE- 3 

dentés, que de la distribution d'un vecteur autour d'un point. 
Nous abordons dans le chapitre III, avec les tensions électro- 
statiques, la considération de la distribution d'un vecteur aux 
divers points de l'espace, ce qui nous permet de déterminer 
les tensions développées dans un milieu électrique anisotrope. 

L'exemple d'une distribution de vecteurs, auxquels nous 
avons donné le nom de vecteurs gauches, dans la première 
partie, se présente dans le quatrième chapitre. Ce dernier 
chapitre est consacré aux phénomènes vibratoires auxquels on 
attribue la lumière. Les équations différentielles dont ils dé- 
pendent peuvent se déduire soit des formules usitées en 
électroinagnétisme, soit de celles employées en élasticité, mais 
dans ce cas, îi la condition d'imaginer un milieu doué de 
propriétés spéciales. L'emploi des coordonnées obliques nous 
conduit rapidement à la détermination de la surface des ondes 
électromagnétiques que M. lleaviside avait obtenue le premier 
et il laquelle les expériences de M. Hertz ont donné une impor- 
tance nouvelle. 

On aura une idée plus précise de l'ensemble du travail en 
jetant un coup d'ceil sur la suite des différents paragraphes 
qu'il contient. 

PREMIÈRE PARTIE 

PropriâtéB dei fonctions vectorielles linéaires. 

Ch. I. — Fonctions vectorielles relalives àdes axes rectangulaires. . 
g 1. — Scalaîfes et vecleurs. 
S 2, — Fondions veclorielles limiaires. Notations. 
S ^. ~ Fonction orthogonale. Formules de Rudrigues. 
S 4. — Ouadriques et ellipsoïdes. Surface des ondes. 
S 5. — Décomposition en somme. Déformation générale, 
g r>. — Changement d'aies de coordonnées. Fonctions sexlinômes. 
§ T. — Invariants. Cubiques des surfaces. 
S 8. — Décomposition en produit. Cas limites. 
S 9. — Théorèmes généraui sar le calcul 9jmbolif(ue. 

Ch. II. — Fonctions vectorielles relatives à des axes obliques. 

g 1 , — Relations entre tes coefticients. Contra ordonnées. 
§ 2. — Changement d'aies. Rotations. Invariants. 
§ 3. — Formes canoniques. 
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DEUXIÈME PARTIE 
Applications . 
Ch, I. — Énergie de déformation. Élasticité. 

§ 1. — État acluel de la Ihéorie. 

g 2 — «tymélrie rechngle. 

g 3 — Dérormdlions et pressions dans les milieux à symétrie oblique. 

^4 — f-neigie ehstique dans les milieux à symétrie ubiique. Dis- 

Ch. II. — Piezoelectnciie et pyroélectricité. 

gl — H^iioUiese et classilication qui en résulte pour les cristaux 

piézoélectriques, 
g 2. — Ponction vectorielle liant le déplacement électrique aux tensions, 
g 3. ~ Application de la méthode, 
g 4. — Phénomènes ilo réciprocité, 
g 5. — Pyroéleclricité. 

Ch. III. — Surfaces équipotentielles et tensions électrostatiques dans 
les milieus anisotropes. 

g 4. — Généralisation des vecteurs solénoidaux et lamellaires. 

g 2. — Lignes de force dans les milieux anisotropes. 

§ a. — Définitions des déplacements électriques. 

g 4. — Tensions électrostatiques dans un diélectrique aiiîsotrope. 

Ch. IV. — Ondes éieciro magné tiques dans les milieux anisotropes. 

Si. — Énergies potentielle et cinétique. 

g 2 et 8. — Equations différentielles ilu mouvement, 

g 4. — Surtice des vitesses. 

g 5. — Surface des ondes. 

S G. — Polarisation rolatoire. 

g 7. — Examen d'un cas particulier. 
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FONCTION VECTOBIELLE. 



PREMIERE PARTIE 

PROPRIÉTÉS DES FONCTIONS VECTORIELLES LINÉAIRES 



CHAPITRE I 
Fonctions rapportées t\ des axes rectangles. 



§ 1. — Scalaires et vecteum. 

Certaines grandeurs pliysiques, telles que la densité, la 
fonction potentielle, la température, ne comportent dans leui 
définition aucune idée de direction dans Tespace. Elles sonf 
déterminées en un point par un nombre p (|ui leur sert de 
mesure, et en tout point par une relation unique, 

entre p et les coordonnées du point considéré. Ces grandeurs 
non directives sont dites scalaires. 

A d'autres grandeurs, telles que la vitesse, la force électro- 
motrice, le flux calorifique, s'attache, au contraire, l'idée de 
direction. Elles sont alors déterminées en chaque point par la 
longueur et la direction d'un segment de droite porté ti partir 
du point où on les considère, et, pour les définir analytiquement 
en tout point, il faudi'a donner trois relations entre les projec- 
tions du segment qui les représente et les coordonnées du 
point. Les grandeurs de cette catégorie sont appelées des 
vecteurs. 

Il peut arriver qu'en un point une grandeur physique de 
nature vectorielle dépende de la valeur d'une autre grandeur 
aussi de nature vectorielle. Cette dépendance est exprimée 
analytiquement par trois relations entre les trois projections 
de chacun dos deux vecteurs. 
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) B. ËLIE. 

j'adopterai désormais la notation à indices et désignerai 
r w„ w„ w,; u„ v„ v, les projections des vecteurs. Ces 
ojections sont alors liées par trois équations de la forme 

ç(«!„ «?,, w,, I',, lî„ V,) = 0. 

Supposons maintenant que dans les développements, par la 
rie de Taylor, de iy,„ on ne conserve que les termes du 
emier degré en v,„; les équations précédentes deviennent 
js simples. C'est leur étude qui va faire l'objet de notre 
emière partie. 

§ 2. — Foncti07t vectorielle linéaire (•). Notations. 

Lorsque les projections d'un vecteur w sont des fonctions 
léaires de ceux d'un autre v, c'est-à-dire lorsque l'on a 

( U», = 6„t', + 6,j!î, -h ftist'a, 

) ) w» = fin"! + ''ii^'î + ('«"i. 

' w, = btiVi ■+■ frjjf'j + bjjî'j, 

dirai que le vecteur w est une fonction vectorielle linéaire 
! V, et je remplacerai ces trois équations par des symboles 
us simples; j'écrirai par exemple : 

!r„, = I 6,, fr,j (>[, I r,ij, 

bu &„ &« 

! K, b,, b,, I 

■rsqu'il y aura lieu de mettre les coefficients en évidence; 

10= \b\v, 
ins le cas contraire. 
De cette fonction, que j'appellerai directe, on déduit deux 

.itres ; 



(<) C'est IlamiUon qui, le premier, a introduit cette dénomination de fonction 
«Eorielle Ilnéaii'e. Le traité élémentaire des qaatemions de M. Tait contient les 
lancés de tliéorèmes l'clalirs à cette fonction. 

Sir W. Thomson et M. Tait l'ont iilili^éc dans le premier volume de leur 
ivrage: Natural Philoaophij, à propos de questions de ciné roaliqne. 
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FONCTION VECTORIELLE. 7 

L'une, dite inverse, résulte de la résolution des équations (i) 
par rapport à u,„. Je l'écrirai 

v=\b\-'w, 

ou explicitement, en représentant par Dj le déterminant formé 
avec le coefficient 6, 

^ |&„fr„ — b„îi„ Kbi,— b,,b„ 6,,(ï„— 6„ftiJ 

Je désignerai par &"' les coefficients de la nouvelle fonction, 
c'est-à-dire les différences inscrites divisées par le détermi- 
nant D». .l'appellerai ces nouveaux coefficients les inverses 
des b. On aura donc 

L'autre fonction, dite transverse (ou conjuguée), sera repré- 
sentée par \ft I et se déduira de la fonction directe en renversant 
ses coefficients relativement î'i la diagonale du déterminant ou, 
ce qui revient au même, en permutant les indices. 

Il y a aussi lieu d'adopter des dénominations et des symboles 
pour les cas où les coefficients sont soumis îi certaines con- 
ditions. 

Ainsi, le déterminant D^ peut être symétrique; sa l'onction 
est dite alors symétrique (') (ou autotransverse ou autoconju- 
guée) et s'écrira \&/. 

Si les coefficients de la diagonale sont nuls et les autres 
égaux et de signes contraires deux à deux, c'est-à-dire si 

^1, = ~ K, K = — fc... ''.. = — fr.i. 
la fonction est dite gauche et s'écrit : \t/. 

Si les coefficients satisfont aux mêmes relations que les 
neuf cosinus des angles de deux systèmes d'axes rectangulaires, 
la fonction est dite orthpgonale et représentée par \b\. 

Enfin, nous verrons que par des changements d'axes on 

<>) V. Tait, Qualeridoits. 
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■ansformer une fonction symétrique en une autre ne 
lant que les termes de la diagonale, et la fonction 
e en une autre où les coefficients vérifient les relations 
à propos de la définition. des fonctions gauches. Ces 
, dites canoniques, s'écPÎPOnt \&/ et |&|. ' 
îsuraé, les huit types de fonctions à distinguer sont les 
;s : 
I b I fonction directe générale, 



\b\-' 


-— inverse. 


W 


— transverse, 


W 


— symétrique, 


\s; 


— gauclie, 


w 


— orthogonale. 


\S; 


— symétrique canonique 


Ul 


— générale canonique. 



définitions, on peut tirer immédiatement quelques 
uences ; 

équations (2) montrent que, si une fonction est syraé- 
son inverse l'est aussi ; 

nverses \b\~' et \fcî""' d'une fonction et de sa transverse 
ansverses l'une de l'autre; ou encore, l'inverse d'une 
n transverse est égale à la transverse de son inverse, 
t par là que l'on pourra faire porter le signe \ de la 
wsion, soit sur les coefficients, soit sur la fonction 
me, absolument comme le signe — 1 de l'inversion. 
3ition des vecteurs est commutative. La somme d'une 
n et de sa transverso 

\b\ + \b\ 
■ fonction symétrique. 

différence : ]&| — \6| est une fonction gaucho. 
1 élimine v, u..., entre les fonctions 

t» = iH», . = |t|«, .., u=ld\l, 
itat 

«> = |6| lt| ... |(J1 î 
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FONCTION VËCTORIBLLE. 9 

sera appelé le produit des fonctions j6j, \c\, .... jd|. Remar- 
quons dès maintenant que, des équations 

w^\b\v, v = \b\-'w, 
on déduit 

«. = |6||6|-'|.; 
d'où 

IHI*|-' = l. 

Je désignerai par une grande lettre affectée des indices 1, i, 
3 et accentuée ou non, les vecteurs dont les constituants sont 
les termes des rangées verticales ou liorizontales du détermi- 
nant d'une fonction. Si l'on considère deux fonctions \b\ et \c\, 
on pourra donc former les tableaux suivants : 

b; b; b„ c; c; c;, 

B, &,, fe„ &,„ c, e„ c„ c,„ 

B, b„ 6„ (»„, C, c., c„ c„, 

Bj tj, fc„ fcj„ G, c„ c„ c„. 

Les valeurs explicites des tonnes du produit |6| \c\ seront 

alors 

[6,,c„+&„Cj,-+-(*„c„ 6|(C,,+&,,c,j+6|,c„ b,,c,,4-&,,c„+fc,jC„ 
(3) ! 6,iC„ +fc,.c„ +fc,,c„ 6uC„+6,,c„+fc,j(;„ 6,,c„+fc„c,,+fc,,c„ 

ou encore (') 

B, c; S B, c, B, c; 
B, c; B, c; B, c, 
S B, c; S B, g; B, G,. 

il en résulte que, si les fonctions \b\ et \c\ sont transverses 
l'une de l'autre, c'est-à-dire si 

c; = B„ c; = B„ c; = b^, 

les fonctions \b\\b\ et \b\\b\ sont symétriques sans être 

(<) \ous einpruntons â la théorie des qualemions quelques dénominations et 
rguelques symboles abréviatiFs : nous posons 

x, y. + x.y, + x .ii. = Sxj 
Kx; + x; + x; =«x. 

9 et ^ sont respectivement le scalaire et te tenseur de Hamillon. 
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10 B. ftUE. 

égales entre elles. Cela n'aurait lieu que si la fonction \b\ 
était symétrique. Il est facile de vépifier les identités suivantes, 
qui seront utiles par la suite : 



(4) 



i D^ SBr'Br' = Sb^Sb',' — [Sb,b;]S 

i D? S Bt B7' = S b; b;3 b; B, — 3 b;» b;b;. 



Si les scalaii-es formées ayec (ioux colonnes verticales d'une 
fonction sont nulles, il en est par conséquent de même des 
scalaires formées avec les lignes horizontales de son inverse. 

Il ne m'a pas paru nécessaire de développer les calculs en 
supposant quelcon(iue le nombre des variables. Les formules 
auraient été tout aussi symétriques et n'auraient pas occupé 
une place bien plus grande; mais, il n'aurait pas été possible 
de les donner sous une forme plus explieito, sous une forme 
pratique. Les développements n'auraient pas pu être effectués 
et ils ne sont utiles pour nous que dans le cas de trois variables 
et parfois seulement de six ou de neuf variables. 

Relativement à ces derniers cas, les pages suivantes et notre 
travail sur l'élasticité, paru dans les Mémoires de la Société 
des Sciences physiques et naturelles de Bordeaux, contiennent 
tout ce dont nous avons besoin. Nous avons voulu nous borner 
à ce qui est absolument indispensable au physicien. S'il s'était 
agi de considérations mathématiques pures, nous accordons 
qu'il aurait été préférable de traiter la question à un point de 
vue plus général. Nous n'avions point à le faire ici. 

Toutes les définitions et les notations précédentes, ainsi tjue 
les conséquences qui en découlent, ont été données dans le 
cas d'une fonction à trois variables; mais elles s'appliquent au 
cas d'un nombre quelconque de variables, et en particulier de 
six ou neuf quantités, fonctions linéaires de six ou neuf autres. 

Il ne sera pas inutile de montrer dès à présent t'usage des 
fonctions vectorielles. Je laisserai de côté leurs applications à 
la géométrie, telles que les formules de transformation de 
coordonnées ainsi que celles qui relient les coefficients d'une 
piême quadrique rapportée à deux systèmes d'axes, pour me 
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limiter exclusivement à celles relatives à la physique. Je 
donne dans le tableau suivant l'énoncé des principaux phéno- 
mènes physiques pour lesquels la fonction vectorielle peut 
être utilisée, en indiquant poui- chaque phénomène ta signifi- 
cation des vecteurs u et w et des constantes 6. 

U peut paraître, au premier abord, assez étrange de voir 
dans un môme tableau et comme applications d'une même 
théorie, des phénomènes aussi distincts que la chaleur, l'élec- 
tricité, la lumière et le magnétisme. 

Ce n'est pas que par l'emploi des fonctions vectorielles, 
nous soyons conduits à assimiler ces phénomènes ou à les 
identifier. Les fonctions vectorielles se présentent simplement 
comme un instrument puissant avec lequel nous pouvons 
essayer de pénétrer dans la nature de chacun de ces phéno- 
mènes. Les fonctions vectorielles nous permettent de raisonner 
d'une façon uniforme dans des cas qui, en fait, peuvent être 
de nature essentiellement différente. L'ensemble des applica- 
tions où il y a lieu d'employer cet instrument semble indiquer 
toute son importance. 
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NATURE DU PHÉNOMÈNE 


SIGNIFICATION DU VECTEUR v 




Foactions générales 


Élasticité, 


Distance de deux points voisins : dx, 




Vitesse en un point : u. 


El ectrod y namiqiie , 

■' 1 milieux 


Force électro motrice, 

Variation de la température : -r— i 


Chaleur, i anisotropes, 




Fonction 


Élasticité, \ 


Normale à un élément, 


Électrostatique, 1 


Force électrostatique. 


Magnétisme, ( ^ju^u^ 


Force magnétique, 


Chaleur, anisotropes, 


Vecteur d'un point, 


LiiEniëre, 

(théorie de Fresnel) 


Déplacement de l'éther, 


Thermoélectricité, 


Flux de chaleur, 




Fonclioti 


Électromagnétisine, 


Intensité d'un courant, 


Id. 


Vitesse d'un point du conducteur, 




Fonction 


Élasticité, 


Déformations -;— dans un système d'axes, 
dx 


Id. 


Pressions X,,..., 


Id. 


Déformations dues à des pressions. 


Piézoélectricité, 


Pressions dues à l'électrisalion, 




Dnitz.dbvGoogle 



FONCTION VECTORIELLE. 



SIGNIFICATION DU VECTEUR w SIGNIFICATION des COEFFICIENTS b 



't (trindmes). 

éplacement relatif de ces points : du, 

K^élération en ce point : ii, 
nlensité du courant, 



Vitesse des déformations : -r- • 

dx 

Coefficients de conductibilité. 



lux calorifique. 

Dans les deux derniers cas, les inverses de h sont les résistances. 



^'métriques : \b/. 

'ression sur cet élément, 

léplacement électrique, 
nduction magnétique, 
ecteurde ce point après dilatation, 
'orce correspondante, 

léplacementélectrique. 

auches : \('/. 

orce ponde romotrice d'Ampère, 
'orce électromotrice d'induction, 

sxtlnûmes : |I>{4. 

•éformations dans un autre système, 
'ressions dans un autre système, 



'ressions correspondantes, 
léplacemeut électrique, 



( pressions sur 3 
ments rectangles. 

Constantes diélectriques. 

Coefficients de perméabilité. 

Coefficients de dilatation. 

Constantes d'élasticité de l'éther. 

Coefficientsthermoélectriques. 



I Composantes de la directrice d'Ampère. 



Fonctions des cosinus des axes. 

Id. 
Coefficients d'élasticité. 
Coeflieients piézoélectriques. 
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g 3. — Fonction orthogonale. — Formules de Rodriguea. 
L'étude de ces types de fonction n'est autre que celle bien 
connue des équations reliant les coordonnées d'un point 
relatives à deux systèmes d'axes rectangulaires. Si on pose 

«'.M = ! Pi. Pi. Pu h'iM 

j P,t Pu Pï3 

ou 

w = [p]v, 

on a, par définition, les douze relations 
Spî — Spl =Sps' = 1, Sp,Pj = S/>3p, =ôp,p,= o, 
Sp^= — 1, Sp;pa =^0. 

11 en résulte que le déterminant D, = 1, que ses mineurs 
sont égaux aux termes qui leur correspondent, que l'inverse 
de [p] est le même que son transverse, c'est-à-dire que 

Les p sont des fonctions de trois variables pour lesquelles 
on peut adopter, soit les angles d'Euler, soit ceux de Rodri- 
gues. En employant ces derniers, on est amené à considérer 
w et ui non plus comme un même vecteur rapporté à deux 
systèmes d'axes, mais à regarder le vecteur w comme résul- 
tant de la rotation du vecteur v d'un angle donné autour d'un 
axe donné. 

Comme j'utiliserai par la suite cette manière de voir, je déve- 
loppe les formules de Rodrigues. Elles résultent de ce qu'une 
fonction vectorielle orthogonale peut se mettre sous la forme 
w^\a\v 
I (ï,, = /f (1 — cose) + cosO 

«it = fifi (* — cos6) + /isinS 

t ff„ = fJ^ (1 — cos 6) — /.sine 

I "il — 

(5) 



«.. = t.f, (I - 


- cos 8)— /.sine 


",.= A'(i- 


- cos 6) + cos 6 


»,.= ^./,(i- 


-cos 6) + /jsino 


».. = r.f. (1 - 


- cose) + /, sin 


«.. = f,r, (1 - 


-cosô) — ^,sine 


"..= flU- 


- C0S6) + cose 
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/,„ étant les composantes d'un vecteur unitaire. Ses coefficients 
satisfont, en effet, aux relations qui définissent une ronction 
orthogonale. Le vecteur / est Taxe de rotation, car les coor- 
données d'un de ses points sont les mômes dans les deux 
systèmes d'axes v et w>; autrement dit, ces équations sont 
satisfaites pour des valeurs de w,„ et v,„ proportionnelles à /,„. 
Ënfm, si cet axe se confond avec celui des v„ auquel cas 
f^:^f^^^ (I, /", = 1, on voit que 6 est l'angle dont il faut faire 
tourner Ow pour l'amener sur Ov,. Pour des valeurs cons- 
tantes de /■„, et 6, la fonction exprime donc que les vecteurs ic 
ont tourné d'un angle 6 autour de l'axe f, pour venir se 
confondre avec les vecteurs v. 

On obtient la fonction inverse en changeant le signe de H. 
Cette dernière exprime que l'on peut amener les axes v sur les 
axes w, à l'aide de la rotation autour de l'axe f. 

§ 4. — Quadriques et ellipsoïdes. — Surface des ondes. 

Si les vecteurs v et w représentent, au point de vue concret, 
respectivement un agent et l'effet qu'il produit, j'appellerai 
énergie du phénomène le produit de ces vecteurs par le cosinus 
de leur angle, produit qui peut se représenter par 

ou, explicitement, par 

b„vl + btiVl -i- fcja!'^ 

Lorsque l'on considère y„, comme les coordonnées d'un 
point, cette quantité, égalée à une constante que je prendrai 
égale à l'unité, représente une surface qui n'est pas nécessai- 
rement fermée et que j'appellerai quadrique inverse. Les carrés 
de ses rayons vecteurs sont les inverses de la projection de w 
sur V. 

Le lieu des projections des extrémités de v sur to s'obtiendra 
en remplaçant dans l'équation précédente les b par leurs 
inverses b~', et sera la quadrique directe Qj. 
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La considération de la fonction transverse donne les mêmes 
quadriques. 
Remarquons de plus que 

0d|6|d = 1 est la quadrique inverse de |6| et \b\ 
et la quadriqne directe de |6|-'; 

Sv\b\ti>= i est la quadrique inverse de \b[-' et \h \-' 
et la quadrique directe de \b\. 

On obtient une autre surface en cherchant le lieu de l'extré- 
mité du vecteur v, lorsque l'extrémité de w décrit une sphère 
de rayon unité, c'est-à-dire lorsque l'on a 

ou 

tfîSBi'+DÎfiBi' + DîSB;* 

+ 2r,p,SB;B; + 2p,r,5B;B;-t-2f,c,9B;B,= 0p|B>=l, 

en posant, d'après les conventions des § 1 et § 2, 

( S b;' = b\, + bl + bU, ÔB, b; = fc„&., + b„b„ + ft„6„, 

(6) j Sb;' = 6Î, -h 6Î. + b\^, S B;Bi = 6„fc„ + fc„6,. + fc„fc,„ 

( 3 Bi' = h\, + b\, + bl, ÔB.B, = fc.,&„ + b„b,,+ b„b,,. 

Cette surface est un ellipsoïde que j'appellerai ellipsoïde 
inverse. Comme d'après les formules (S) ses coefficients sont 
les mêmes que ceux de l'ellipsoïde inverse delà fonction |fe|\fc|, 
son équation peut encore s'écrire 

E(, âr|6|\6/p = l. 

On appellera ellipsoïde direct de \b\ le lieu de l'extrémité du 
vecteur w, lorsque l'extrémité du vecteur v décrit une sphère 
de rayon unité. Il suffira de changer b en b'' dans les équations 
précédentes, ce qui donne, pour équation de cet ellipsoïde, 

E^, &v\b\-'\b\~'v=2l. 

Les ellipsoïdes de la fonction transverse de \b\ s'obtiendront 
en permutant les indices des b dans les équations précédentes. 
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ce qui transforme les B' en B et conduit, pour l'ellipsoïde 

inverse, à 

Sv\^\v=Sv\b\\b\v = i, 
et pour l'ellipsoïde direct, à 

Sv\b[-'\b\-'v = \. 
Remarquons encore que 

0ti|{)l\ti|f — 1 est Tellipsoïde inverse de \b\ 

el rellipsoïde direct de \b\-', 
3ii|6|"'\6|-' = 1 esl l'ellipsoïde direct de |&1 

et l'ellipsoïde inverse de 1&|~'. 
De même 

6v\b\ \b\ est l'eilipsoïde inverse de \6| • 
el l'ellipsoïde direct de \b\-'. 

Nous verrons plus tard que les ellipsoïdes inverses de cpb 
fonctions ne diffèrent que par leur orientation, et qu'il en est 
de même des ellipsoïdes directs. 

Lorsque la fonction est symétrique, c'est-à-dire lorsque 
6j, ^ b„, ..., on peut, par un changement d'axes, faire dispa- 
raître les doubles produits de la quadrique inverse et mettre 
la fonction vectorielle sous la forme canonique 

(7) Wi = b,v„ Wt = btVi, i", — bjPa ou w=\h!v. 
Les quadriques inverses et directes deviennent alors 

(8) 0.-. Mî+MÎ + 6,t-î-l, ^-.f^ + ^^1, 0„ 

et les ellipsoïdes, 

(8') E„ blvt + blvi-hblvl = i, ^H-^-H^=i E^. 
Oî bl oj 

Ces derniers sont polaires réciproques par rapport à la 
sphère de rayon unité; il en est de même des quadriques, et 
les axes des uns sont les carrés de ceux des autres, il en 
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iiUe entre les rayons vecteurs des relations simples que je 
i indiquer. 




1 la direction dont les cosinus directeurs sont ni„„ corres- 
id dans l'ellipsoïde E^ la direction conjuguée (b*»î),„=jn[,3, 
i celte dernière correspond dans l'ellipsoïde Qî la direc- 

■(f),..=- 

oient : 0) =»îOm' l'angle de ces directions; 

)i, Qrfi K.-, Ed leurs points de rencontre avec les quadrlques 

Bs ellipsoïdes; 

. et K' leurs intersections avee les plans tangents en Oi et Oj, 

is pei^endiculaires au plan de la figure; 

lX,„, rx,u, pÇ,„ les payons vecteurs de Ej, Q^ et K et leurs 

rdonnées pour une dii-ection m ; 

l'X',„, r'x\„, p'ï'„„ ceux de E;, O^et K' sur la direction m'. 

In a symboliquement pour les cooi'données des divers 

its : 

K, ; = »((:, Qd, x = mr, E^, X = m R = frm', 
K', ^'=m'p', Q;, x'=m'r', E^, X'=m'R'=|, 

secondes valeurs des X se déduisant des équations (7) 
que le vecteur v ou lo est égal à l'unité. 
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Les définitions précédentes fournissent de plus : 

coaw — m,tn; + ni,m', + mjWij = âmm', 

X, X, X, /m' mî mî\„ R 

eos to — m, Y^ + m, -r^ + ni, -r^ = ( -j-i 4- T^ + -r^ I R = -=1 

0, p, 0, \o, 0, 0, / r* 

cosh; r=6,XJffli + btX',m'i + frjXJiMJ 

Les triangles semblables nous donnent 

p' p , 
coswi: ^ = £;; 

et enfin, la réciprocité de Q, et Od, 

rp = l, r'p' =1. 

L'élimination de p et p' entre ces relations conduit à 

r' =rR, r = r'R, RR' = 1, rr'cosio— j. 

Ces relations permettent d'arriver rapidement à une surface 
(les ondes si on leur adjoint les hypothèses suivantes, corres- 
pondant aux hypothèses admises par Fresnel dans la théorie de 
la lumière : 

1" Un déplacement du vecteur v égal à l'unité dans la direc- 
tion m' met en jeu une force fonction vectorielle du déplace- 
ment et dont les composantes sont, par conséquent, 
X,„ = (bm\„. 



â" La projection Kfc,m;' + (*,»»,' + 6,™;* ou ôXm' de cette 
force sur m' est seule efficace pour propager une onde plane, 
dont la normale Oq' est perpendiculaire à OK'. 

3" De même que dans la théorie de l'élasticité des corps 
pondérables, la vitesse de propagation est la racine carrée de 
la force etficace dont la valeur est 

ÔXm' =\/^ = l^R^^^ = OYi'. 



5x»' = »/^ 



Si Og'=OK', le lieu des extrémités de q', lorsque la 
direction m' varie, s'appelle surface des vitesses. Soient q\,. 
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les cosinus directeurs de Og', f' son rayon vecteur; Og' étant 
perpendiculaire à OK' et parallèle à K'Q^, nous avons 

g[m', •+■ g^m', + qim', ou Sm' g' = 0, 
Ç; — a;, : ï; — X, : Çj -a?, = vi : «i : îi- 
D'autre part, en supprimant les indices 1, 2, 3, 

R P 

L'élimination de m' entre ces équations donne pour équation 
de la surface des vitesses en coordonnées polaires : 

(9) -T/Tr + ^ + ^=»- 

P — 1*1 P — *i P 0) 

La surface des ondes est Tonveloppe des plans des ondes, 
c'est-à-dire du plan normal à la figure et passant par qq' . On 
sait que le point de contact du plan qq' avec son enveloppe 
est le sommet q d'un triangle qOq' semblable au triangle OdOK'. 

Soient g,„ les cosinus directeurs de Og, r ou — son rayon 
vecteur égal à OQ^; nous avons, comme précédemment, 

I, — *; : Ç, — ic, : Ç) — ^j = Vi : 9i : ?), 

d'où 

r' fc, r" fc' r" fc, 
équation de la surface des ondes. 

Dans la seconde partie, nous aurons à étudier une surface 
des ondes plus générale déduite de la théorie électromagnétique 
de la lumière. 

§ 5. — Décomposition d'une fonction vectorielle en une somme 
de deux autres, l'une symétrique, l'autre gauche. — Défo7-- 

mations. 

La décomposition en somme résulte simplement de ce 
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que l'on peut mettre la fonction vectorielle (1) sous la forme 
suivante : 

2 2 " 

l)„+b„ . b„-i-b., b„--b„ 

m 




«;- 2 + î ■ 

La première partie de cette somme est une fonction symé- 
trique et représente un vecteur dont la direction est conjuguée 
du vecteur v relativement à la quadrique inverse : Qv\b\v=^i. 
Quant à sa grandeur, c'est l'inverse de la normale abaissée de 
l'origine sur le plan tangent au point V à une quadrique 
semblable à la précédente et passant par l'extrémité V du 
vecteur v. A l'aide d'un changement d'axes, annulons les 
coefficients b„ + b^^ des doubles produits de la quadrique; la 
fonction vectorielle se ramène à sa forme canonique 

Ces équations expriment que tous les points d'une sphère 
de rayon unité sont venus se placer sur l'ellipsoïde 



Les points du milieu ont donc été soumis à une triple 
extension de valeur &„„ suivant les trois axes de la quadrique 
inverse. Telle est la signification géomélrique de toute fonction 
symétrique. 

La deuxième partie de la somme est une fonction gauche, 
qui représente un vecteur perpendiculaire à v et au vecteur r, 
dont les composantes sont 
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Fig. !. 

Sa grandeur est, d'après la forme des deux derniers termes 
des équations (10), 

®v.'3trsinOV.OB, 

c'est-à-dire le segment VW sur la figure 2, dans laquelle les 
trois vecteurs r, v elw sont représentés par OR, OV et OW. 
Soient en effetTP la perpendiculaire abaissée de l'extrémité de V 
sur OR, 9 l'angle VPW ; on a 

VW — yp langO — OT. sinÔYTÇB. lang», 

c'est-à-dire l'expression précédente, si 

. '3[r= lange. 

Mais VW résulte d'une rotation amenant V en V, et d'une 
dilatation VW perpendiculaire à OR et telle que : 

"FV:^PWcose. 

Si donc le vecteur w est une fonction vectorielle générale 
d'un vecteur v, le premier se déduit du second en supposant 
que les points de l'espace lieu des extrémités des vecteurs v, 
qui ont la même origine, aient subi : 

l" Une triple extension suivant les trois axes de la quadrique 
inverse ; 

2° Une rotation 6 autour de l'axe r, telle quo 
%.r— langO; 

iC Une extension peqiendiculaire à cet axe (radiale) et égale à 
i 



VW ou PV 1 



cosO/ 
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Ce modo de décomposition appliqué à une fonction ortho- 
jronale conduit à une simple rotation. En effet, l'ellipsoïde 
d'extension dont on évaluera les coefficients à l'aide des angles 
de Rodrigues est alors de révolution autour de l'axe de 
rotation. Ses axes sont égaux l'un à l'unité et les deux autres 
à eos 6, c'est-à-dire qu'il se produit une contraction qui 
compense l'extension radiale. 

Pour des déplacements finis, ce mode de décomposition en 
somme n'est pas avantageux, par suite de la présence de 
l'extension radiale. Nous aurons à signaler plus loin une 
décomposition en produit qui ne fait intervenir qu'une triple 
extension et une rotation, quels que soient les déplacements. 

Mais lorsque le vecteur r est un infiniment petit de premier 
ordre, l'extension radiale rapportée ù l'unité de distance, c'est- 
à-dire — ~rx — étant infiniment petite du deuxième ordre, est 
négligeable. On peut écrire alors 

(ange = 6 = '3tr, 
, et la fonction vectorielle représente une triple extension ac- 
compagnée simplement d'une rotation. 

Cette simplification se présente lorsque l'on considère les 
déformations d'un milieu élastique. Soient, en effet, xelx + dx 
les coordonnées de deux points voisins, x+u et x+dx+u+du 
ce que deviennent ces coordonnées quand le milieu a subi 
une déformation. Le vecteur qui relie ces deux points dans 
leur nouvelle position est une fonction du vecteur dx, qui 
peut s'écrire, en ne conservant que les premiers termes du 
développement. 



d{x + u) = 



(10') 



du, 

àx. 



dXy 



dx. 



dXt 



du, 
dXg 
du. 
àx, 
, du^ 



dx. 



Il suffira, pour appliquer à ce cas particulier les théorèmes 
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ppécédents, d'identifier les coefficients des fonctions vectorielles 
(10) et (10'). A cause du fréquent emploi que vont recevoir 
dans ce qui va suivre certaines expressions formées avec ces 
coefficients, il est bon de rappeler le sens de ces expressions. 

Les neuf coefficients de la fonction (10') sont les cosinus 
directeurs des trois éléments de ligne dx„ dx„ dx,. Ils 
permettent de calculer les cosinus de ces éléments entre eux, 
et, en remarquant que les -p sont des infiniment petits du 
premier ordre, on trouve ainsi pour leurs valeurs : 

''i\âx,'^ dxj' ^\dio,^' dx,r i\àx,'*' dxj' 
ou bien 

b„ + b,, &„+<>., 6,1 + fc.. 

i ' 2 ' 2 ■ 

Je représenterai par 3,, g„ 3, ces quantités que l'on a 
appelées glissements (de Saint-Venant). 
Ses dérivées 

du, <ïw, au, , , , , . , ■ 

ff^i dXj aie, 

représentent les allongements par unité de longueur des seg- 
ments dx,, dx,, dx,; on les appelle les dilatations et je les 
représenterai par d„ d„ d,. Enfin, d'après ce qui précède, les 
différences 



dxj' i\dx, dxj' i\d. 






ou bien 



ne sont autres que les composantes de la rotation déjà repré- 
sentée par r„ r„ r,. 

Les six quantités d,„, g,„ sont les coefficients d'une fonction 
symétrique que je représenterai dans certains cas par Xlj. 
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§ 6. — Fonction vectorielle rapportée à de nouveaux axes. 

La question à résoudre est la suivante : Tpansformer la 
fonction 

(11) 2^ = 161.1; en la fonction |' = |0[5 

à l'aide des substitutions orthogonales 

Il y a avantage à considérer la question dans le cas de substitu- 
tions quelconques, la solution devant nous être utile par la suite. 
Soient donc posées les fonctions générales 

x- = m-'^' et x=m. 
Elles transforment x' = \b\x en 

égalité qui peut s'écrire 

dU * V=\g\\b\\s\l 

L'élimination de % et %' entre les équations (11) et les for- 
mules de substitution aurait conduit à l'équation 
(11,) ■i''=lî|-'IPI|sl-'a'. 

La comparaison des relations (H), (11,) et (lis) donne alors 
(Hc) I31==lyllft|l«l et \b\ = \q\-'\?,\\s\-'. 

Ces égalités sont des fonctions vectorielles à neuf termes. 
Le développement de la première conduit aux neuf équations li- 
néaires dont les coefficients sont inscrits dans le tableau qui suit : 
b., 6.! f>„ b„ 6„ 6„ &„ b„ b„ 
Pu ?ii*ii 9ii*ii 5ii*ji Vii^ii ?ii*ii ?ii*M 9ijSii Vij*ii în^M 
Piï 9u*ii 5ii*i! ?ii*3i îii'u 5it*ii ?ii*ji îu^ii 5ii*ii 5i»**i 
3t» ?ii*is 9ii*i» Sii'aa Vii*i3 9ii*M 9ti*3J 9ij'u SiJ*n Sia'j» 

/JO^ Pli î.l'll îït»ïi îïl'll ?îi«il ?»»» ?.i«M ?M*M ¥«««1 9>a*« 

P« îtlSlt 3,1*11 ?11«»1 îll«l, ?M*»1 îll*31 «««11 ?I.«1. 9ll*S« 

&1» îîlSn Vll**l îll^JS îll^tJ «I1*U 5lî*S3 ïis'll Î,)*M Sll's» 

?St ÎM^U îjl'll 9si*ll 9>i*ll îjl^tl 9)Î*JI Vj)*ii ?J»*J1 ?Sl*il 

&»i Î3l*ll ?ai*tl 9st*3i V3i*ti ?Ji*ii 9ji*si 533*11 ?ii*ii îsa *a( 

^H 9)1*13 Vll^IS 93t*S] 931*19 ?>l'l3 $31*33 933*13 9ll*ll 9l3*3t 
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Dans le cas particulier actuel de substitutions orthogonales, 
nous devons poser : 

I»I = [P] et [q\'' = [p] ou \g\ = [p]-'=\p], 
c'est-à-dire pemplacep dans le tableau précédent s par p et g 
par cette même lettre /j, mais en permutant les indices daus 
ce deuxième cas. 

Désignons papp, o, y et X ce que deviennent respectivement 
dans la fonction transformée les quantités r, d, g et l définies 
plus haut. Le tableau écrit précédemment montre que les 
rotations se transforment par la substitution inverse, c'est-à- 
dire que 
(13) P^[p]-'r ou r=[,,]f. 

Les dilatations et les glissements se transforment par une 
substitution sextindme dont la forme explicite est donnée par 
le tableau suivant : 



i, i. 


II. 


3, 




S, 




9. 


Pi'f P'n 


Pi, 


^P„P„ 




iPnPn 




2P„P„ 


PU Pi, 


P'„ 


2P„P. 




iP„P„ 




2p„P,. 


Pu Pi, 


P'n 


ip„P„ 




3P,.P.. 




2P,.P,. 


i.Pi. P„P„ PuP, 


P,.P..+P. 


P, 


P.,l>,.+P, 


P. 


P„Pi.+P,.P 



(13.) î. 
Te 

ï, PnP,i P«Psi P»Pii P,.P«+P.iP). ftiPH+PoP» P,iPi^+PiiPt^ 
Ts PiiPi, P,iP» PmP» P,iPa+P»P>i P.iPi.+P.tP» PiiP«+PnP,i 
Se représenterai cette fonction par le symbole 
(130 X — [ [p]-'| ( on encore X = | \p] | (, 

afin d'indiquer qu'elle est dérivée de la fonction orthogonale [p]. 
Cette dernière égalité résolue par rapport à ï devient, avec 
les conventions établies, 

i=l[p]-l-'X. 

Mais on peut lui donner une autre forme. En effet, la deuxième 
égalité (Ht) diffère de la première par le changement de g et s 
en leurs inverses. L'inverse de la fonction sextinôme (iS,) 
s'obtiendra, par suite, en changeant dans cette fonction les 
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quantités p en leurs inverses ou en leurs transverses (puisque 
[p] est une fonction orthogonal?), c'est-à-dire qu'il suffira de 
permuter les indices, et, par suite, on a 

l = |[p]|X. 

La comparaison de ces valeurs de î permet d'écrire 

i[p]-'i = i[^]r' 

et établit un théorème relatif à ce genre de fonctions. On peut 
faire porter le signe d'inversion — 1 soit sur la l'onction 
orthogonale [p], soit sur la fonction sextinôme dérivée, sans 
altérer cette fonction. 

Comme application de ces formules, on peut citer, en pre- 
mier lieu, les relations entre les coefficients d'une même 
quadrique rapportée à deux systèmes d'axes, rectangulaires de 
même origine; X et l sont alors les coefficients de la quadrique 
dans les deux systèmes. 

En second lieu, 6 et ^ peuvent représenter les dérivées 
■p et -j- des déplacements y et u d'un milieu déformé rap- 
portés à deux systèmes d'axes différents 5 et x. Dans ce cas, 
il est facile de vérifier, par un calcul direct, les formules (12) 
qui lient ces dérivées. Posons, en effet 

a~\q\u et l'^lsl-'x, 

et le calcul des dérivées conduit au résultat suivant, que 
j'écris sous la forme symbolique habituelle 

^' = (q,,dti^-i-q,,dv, + g,,âu,){^ + ^ + p.), ' = 1,2,3, 

c'est-à-dire aux formules {M). Par suite, les rotations, etc., 
seront liées par les formules (13). 

% 7. — Invarhiiitn. 

J'appellerai invariant d'une fonction vectorielle toute l'onc- 
tion des coefficients h dont la forme reste inaltérée lorsqu'on 
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opère un changement d'axes rectangulaires. Ils méritent d'être 
étudiés spécialement, car c'est avec eux que l'on formera les 
expressions de l'énergie. 

Les invariants résultent de la combinaison d'un nombre 
limité d'entre eux que j'appellerai indépendants et que je vais 
énumérer. 
Un vecteur n'a qu'un invariant qui est son tenseur. 
La rotation se transformant comme un vecteur pour une 
substitution orthogonale, on a donc le premier invariant 
(15) B' = r; -Hrî + rj. 

On en obtient immédiatement trois autres en formant les 
coeflicients de l'équation cubique 

X'-I,X' + I,X-I, = 0, 
dont les racines sont les carrés des axes de la quadrique 
inverse Q^. D'après les notations du § 5, on peut écrire 
I, = d. + (/, + rf„ 
1 1, = (/jrf, + rfjrf, + (/,rfi — gl — gl — Oi, 
(Ï6) \à, g, g,\ 

lh = \g. à,_ gA 
\ 9i 9i '^3 I- 
Nous verrons que la déformation d'un milieu est complète- 
ment défmie, abstraction de l'orientation, par trois dilatations 
suivant trois axes rectangles et une rotation de direction 
donnée relativement à ces axes, ce qui implique la connais- 
sance de six quantités. Nous devons donc trouver deux autres 
invariants outre les quatre précédents. Les coefficients de la 
cubique relative à l'ellipsoïde inverse E,- les fournissent en 
effet. 

Leurs formes étant cortipliquées, je les simplifierai en choi- 
sissant les axes de façon à annuler les glissements. Nous 
aurons donc : 

b„ + fc„=0, fr„ -ft„ = 2ri, 
d'où 

b,, = - t„ = r. .... 
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Les invariants précédents deviennent 

|^ = (f, + rf, + d„ I, = rf,d, + d^d, + (/,rf„ I, = d,ei,d,. 

Les coefficients de la cubique relative à l'ellipsoïde E„ 
calculés à l'aide des coefTicients de cet ellipsoïde donnés 
au § 3, sont : 

J, -i;-2I, + 2R', 

J, = lî - 21,1, -H R* + 2K„ + 2K„ 

ils contiennent les trois invariants nouveaux 

/ K, = rld, + r,M, + rjrfj, 
(17) K„= r;d,rf, + r;d,d, + rjd.rf,, 

liés par l'équation 

K„- RM, -I,K, + K,. 

Comme vérification, on peut supposer symétrique la fonction 

vectorielle |h|, c'est-à-dire annuler r„„ et la cubique devient 

X' - (i; - 21.) X* + (lî - 21, 1,) X - lî = 0; 
c'est réquation aux carrés des racines de la cubique relative à 
la quadrique Qo ce qui implique qu'en ce cas, les axes de 
l'ellipsoïde E.- ri.it les eari-és de ceux de la quadrique Qi, 
comme on l'a déjà vu. 

Les coefficients !'„, de la cubique relative à la quadrique 
directe Qa sont : 

DÎI,' = I, + R', oui =UI, + 2K„, ntV, =IÎ + I,K. 

Quant à ceux de l'ellipsoïde direct E^, on les déduira des 
précédents par un procédé identique à celui employé pour 
déduire les invariants relatifs à Q; de ceux de E,. Ils n'intro- 
duisent aucun invariant nouveau indépendant. 

Ces invariants permettent d'en obtenir d'autres; car on 
pourra, avec les dérivées par rapport à x,„ du vecteur r, 
former de nouvelles rotations, dilatations et glissements sem- 
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blables aux rotations, dilatations et glissements fopmûs avec 
les dérivées du vecteur «, et construire avec ces quantités des 
expressions invariantes telles que celles trouvées précédem- 
ment. Elles interviennent dans l'explication des phénomènes 
physiques pour lesquels il est nécessaire de prendre en consi- 
dération les dérivées du second ordre dans le développement 
en série du vecteur u. 

Soient r', d' et g' ou /' ces nouvelles quantités; leur combi- 
naison avec r, d el g ou l fournit un nouvel invariant que 
j'adjoins aux précédents, parce qu'il a servi <ï M. Jiiveii à 
expliquer la polarisation rotatoire des cristaux. C'est 

(18) G = /„ /;, + /„ (;, + /„ (;, +2/,, („ + 2 (,. ;;, + 2 1„ (;, o. 

Cet invariant provient de la considération des invariants 
simultanés de deux formes quadratiques ternaires, qui, au 
nombre de deux en général, se réduisent, en ce cas, à un seul, 
en raison de la symétrie des quantités l et I'. On peut se 
rendre compte de son caractère d'invariance par le raisonne- 
ment suivant : 

Les quantités / et V se transforment, d'une part, par la 
fonction (13). On sait, d'autre part, que si on opère un chan- 
gement d'axes dans le premier membre 

û„irî + ... 4- 2ir„ic,a;, + ... 

de l'équation d'un ellipsoïde, l'expression obtenue reproduira 
cette même forme, à la condition d'appliquer aux nouvelles 

quantités a\, ... et x\' ... des substitutions de la forme (13). 

§ 8. — Décomposition d'une fonction générale en un produit de 
deux facteurs, l'un fonction symétrique, l'antre fonctmi ortho- 
gonale. — Cas particulier. 

. On a vu, au § 5, qu'une fonction symétrique w = \/i/v 
pouvait, par un changement d'axes, être amenée à la forme 

{') Voir Salmoii-Fiedier, Kegelschnitte. 
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canonique f ,„ ^^Ctu),,,. Les points, ppimitivement sup'tinc 
sphère, viennent donc se placer sur l'ellipsoïde direct 



9 



(î)'-' 



et tes points de l'ellipsoïde inverse 

S {hvy — i 

viennent se placer sur une spliope. 

Cette l'onction représente une triple extension du milieu, dont 
les composantes sont mesurées en grandeur et direction par 
les axes h„ h„ h, de l'ellipsoïde direct, auquel nous réser- 
verons le nom d'ellipsoïde de déformation. On a vu aussi an 
g 3 qu'une fonction orthogonale [p] représentait une rotation 
des points du milieu autour d'un axe. 

Nous nous proposons, dès lors, de chercher les fonctions \hj 
et [p] satisfaisant à l'équation 

16| = [p]\ft/. 

Le problème n'a, en général, qu'un nombre limité de solutions, 
puisque le nombre des paramètres est le môme (neuf) dans 
chacun des membres, la fonction orthogonale [p] étant définie 
par trois paramètres et la foneîîon symétrique \hj par six. 

A. — Déterminalio» de la lonclioo \'i/. 

Remarquons que les opérateurs orthogonaux tels que [p], 
ne produisant qu'une rotation, ne peuvent que modifier les 
orientations des ellipsoïdes directs et inverses, mais non leurs 
formes, c'est-à-dire les grandeurs de leurs axes. Bien plus, 
les ellipsoïdes inverses des fonctions \hl et ji>|, liées par 
|i>| = [p]\hl, sont orientés de la même façon, car ils repré- 
sentent la position initiale des points du milieu qu'une 
extension \fe/ amène sur une sphère et qu'une rotation 
subséquente [p] amène à leurs positions finales. 

Par la résolution de la cubique relative à l'ellipsoïde inverse 
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E( dont les coefficients sont explicitement donnés au paragra- 
phe précédent, on obtiendra ses axes et, par conséquent, 
ceux de l'ellipsoïde de la fonction \h! . 

8,„ étant les racines de la cubique, les directions de ces 
axes seront données par des équations connues du premier 
degré. 

D'autre part, comme ta fonction \/t/ est symétrique, les axes 
de son ellipsoïde inverse sont les carrés de ceux de la quadrique 
correspondante. Le calcul des coefficients de \ft/ se fera donc 
il l'aide des formules (13), dans lesquelles les lettres p repré- 
senteront les cosinus directeurs déjà calculés de ces axes, et 
les quantités f,„ l',„ l'si seront remplacées par s„ s„ s,, les 
quantités î\„ ..., V,,, par zéro. 

Si l'on avait identifié les fonctions [6j et \ft/[p], il aurait 
fallu considérer l'ellipsoïde direct, car c'est lui qui est le 
même comme orientation et comme grandeur pour les fonc- 
tions |bj et \hl, puisqu'il représente la position finale des 
points du milieu après une rotation. 

B, — Caleol des «tsinas dircctenrs /',„ de l'aie de rolatisn. 

Il suffit d'identifier Jes fonctions 

(19.) \b\v et [p]Wr, 

en remplaçant les quantités p par leurs expressions en /" et 6 
(équations 5 du § 3). 

Désignons par [0]*' ce que devient [p]*', lorsqu'on emploie 
les variables de Rodrigues; par B„j et H,,, les trois vecteurs 
dont les constituants sont les termes des rangées horizontales 
de |ft|, \ft/. L'identification donne trois systèmes de trois 
équations que l'on peut écrire 

(I9i) B, = [e]-'H„ B, — [e]-'H„ B, = [«]-'H„ 

ou explicitement pour le premier système, 
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b, = \a\ A. 
;„ = f\{i — C0s6) + cos6 
Il ^ /iA(l — cosO) + /jsine 
!„ =^ /;/;{l — cos8) — ^.sino 
!„ — /;/■,{! ~ cosO) —fjSinfl 

(„ = /"î (i — coso) + C03 e 

i„ — /;/, (I — cosfl) + /isinO 

i„ = /,/,(! — cos9) + /;sine 

t„ =^ fJ^I^X — cos6) — /;sin6 

; a„ — /^ (1 — cos8) + cos8 . 

En multipliant ces trois équations respectivement par f„ f„ /"„ 

et ajoutant, on a 

*iiA + *,./■, + hff^ = Kif, + KJt + KU 

L'adjonction de deux autres équations tirées des deux autres 
systèmes donne les trois suivantes : 

I fc.i K 6m \h»=\K, K Aji |/i«' 

(i9„) 6i, K K A.. K K 

I &ii i-ti ft,. 1 I K K K I 

ou 

\h\f=\hlf, 

et par suite les valeurs de f. 



Pour obtenir 6, opérons avec [p]~' sur la gauche des deux 
membres de l'équation (19„), ce qui donne 

[p]-'\b\v = \hlv, 

et, en identifiant les coefficients, trois groupes d'équations 
semblables î\ {\%) ou (19,). En les retranchant de ces dernières, 
on peut résoudre chacune de celles obtenues par rapport à 
tang 6 et les mettre sous la forme suivante, dans laquelle on a 
écrit & pour b — h: 
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b„ + *„ 


i,, + »„ t„ + *„ 


r,K-r.i. 


f.K-fA, f,K-f,6„ 


s.,-»., 


6,, + »., !>.. + «„ 


r,K-f.\ 


n»„-/',s., 7A.-A5.. 


K + »„ 


6„ + ft„ 6„ + Aj, 



^ lange = 



fc„ + fc,. + />,3 -(Al, + fe„ + A.) 



: cosinus Pi, ... se déduiront de /"et 6 par les formules de 

:lrigucs. 

appliquons cette dénn m position à la fonction vectorielle 

résentée par les équations (0'), en négligeant 'es 1 j- I 

ilipsoïde et la quadrique inverses ont respectivement pour 
lations 

{! +2d,)cî + ■.. +23',P,tf, + ... = i, 

(1 + rf,)p; + ... + j/,f, «,+ ... = 1, 

ir axes 1 + 2î,„ et -1 + e„„ en désignant par s des infini- 
[it petits du premier ordre, et ont la même orientation. 
.es équations (19^) vl (10/) se réduisent dans ce cas à l'ap- 
ximation fixée à 

^* = y* = ^ =: -atr = tang li = e, 

iltats déjà obtenus au g 5. 

In autre cas particulier est celui d'une fonction vectorielle 
che ne contenant comme coefficients que les six quantités 
r,„. Le vecteur w est alors perpendiculaire aux vecteurs r 
', et son tenseur est 

%w = '^r.%v.siQvr. 

lOrsque l'extrémité de v décrit une sphère, celle de w décrit 
deux faces de la surface d'un cercle auquel se réduit l'ellip- 
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solde inverse. On passe du vecteur v au vecteur w,' en projetant 
le premier sur ie cercle, lequel est perpendiculaire à la direc- 
tion r, et en faisant tourner cette projection de 90° autour 
de r. Les formules donnent 

Enfin, si la fonction \b\ est symétrique, la déformation du 
milieu consiste en une triple extension transformant une sphère 
en un ellipsoïde. Cette extension peut ôtre envisagée comme 
s'opérant suivant les axes de coordonnées, à la condition de la 
faire suivre d'une rotation. Ceci revient à décomposer la fonc- 
tion en un produit dont l'un des facteurs [a] produit la rotation, 
et dont l'autre, relatif à l'extension, est une fonction symétrique 
canonique. On les détermine en identifiant les deux membres 
de l'égalité 

\bi = H \ki, 

ou encore 

\bi = \l; [•]-■. 

§ 9, — Théorèmes généraux sur le calcul des fonctions 
vectorielles. 

La fonction vectorielle linéaire est le symbole de certaines 
opérations que l'on peut interpréter soit en géométrie, soit en 
physique. 

Une interprétation géométrique consiste à la considérer 
comme une déformation résultant d'une triple extension, sui- 
vant trois axes rectangles, accompagnée d'une rotation. Les 
fonctions symétriques et orthogonales pourront, dans un pro- 
duit, être appelées respectivement facteurs d'extension et de 
rotation. Le premier facteur transporte sur des droites et des 
plans différents les points placés primitivement sur des droites 
et des plans, et sur des ellipsoïdes les points placés sur des 
sphères. Le deuxième fait tourner les figures autour d'un axe 
invariable relativement aux axes d'extension. 
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ous savons que l'addition de plusieurs fonctions vectorielles 
■espond à l'addition géométrique de vecteurs de même 
ire, et que cette opération est commutative. 
ans les produits, les vecteurs peuvent être de natures diffé- 
es. Un tel produit de fonctions vectorielles est associatif: 

{\a\\b\)\c\ = \a\{\b\\c\), 

s il n'est pas commutatif. 

y a cependant avantage à intervertir les facteurs pour 
'udre certaines questions. C'est pourquoi je vais indiquer les 
lifications que cette opération introduit dans le résultat, 

HÉORÈHE I. — Étant donnée une égalité entre des fonctions 

orielles 

1A|=|B|, 

leut en multiplier les deux membres par une fonction quel- 
que |P| sans altérer l'égalité 

|P1|A| = |P||B|, 

'Ucope opérer avec ces fonctions' sur une troisième, 

!*11P| = |B1|P|. 

l'aide de cette propriété, on pourra isoler un facteur quel- 
jue d'un produit. Si on a, en effet, 

|a||(.||^||c||d[=A, 



I \x\ = \b\~^\a\~'\M\dr\c\-K 

BÉORÈHE II. — Si on intervertit les facteurs d'un produit de 
i facteurs en changeant les fonctions en leurs transverses, 
isultat est la fonction transverse du premier produit. Ainsi, 
a équivalence entre les égalités suivantes : 
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Ceci résulte des valeurs, donoées par les formules (3), des 
coefficients d'un produit de fonctions. 

TiiÉoitËHG III. — Si on ïnterveptit les facteurs d'un produit en 
changeant leurs fonctions en leurs inverses, le résultat est 
l'inverse du premier produit. Il y a donc équivalence entre les 

égalités 

(2J) |1>I = I«M<I| el \b\-' = \dl-'\c\~'. 

Ceci résulte soit du théorème [, soit des identités 
\b\\bl-- = \c]ldlld]-'le\-' = l. 

Si l'on particularise les fonctions, ces propositions condui- 
sent à d'autres propositions nouvelles ou déjà signalées. 

Ainsi, par le théorème II, on déduit de 

in = \»; [pi 

régaliti'- 

puisque 

[Pl-' = \p]i 
d'où, par multiplicatio ii, 

(23) \b\\b\ = \hr- 

Or, d'après les formules (3), le carré d'une fonction symé- 
trique étant symétrique, il en résulte que l'application it un 
milieu d'une déformation générale que l'on fait suivre de sa 
transverse le déforme sans rotation. 

Soit encore 

|l'l = W[pl. 

et par suite 

ll'i-' = W-'W~' ei \b'r' = \hi-'[p]; 

d'où, par multiplication, 

\bi-< ibf = \hj-; 

résultat qu'on aurait pu tirer de l'égalité (23). 
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La permutation de deux fonctions symétriques donne le 
transverse du produit primitif; ainsi on a 

|h| = \A/\K/, \6| = \K/\A/. 

D'après les formules (â), deux déformations successives 
résultant de pures extensions produisent une déformation 
générale. Cependant, si les fonctions pont symétriques et du 
type (10'), le déterminant des formules (3) devient symétrique 
en négligeant les infiniment petits du second ordre (^p) , et 
par suite le produit de ces fonctions correspond à une défor- 
mation sans rotation. 

Lorsqu'on a décomposé nne fonction générale en un produit 
de deux fonctions, l'une symétrique, l'autre orthogonale, les 
ellipsoïdes direct et inverse de cette fonction ne peuvent diffé- 
rer, d'après les formules précédentes, des ellipsoïdes corres- 
pondants de la fonction symétrique que par l'orientation. Bien 
plus, les ellipsoïdes directs seront orientés de même, si le 
facteur d'extension précède le facteur orthogonal, c'est-à-dire 
si l'extension est consécutive à la rotation. Dans le cas con- 
traire, ce seront les ellipsoïdes inverses qui coïncideront. 

Le produit de plusieurs fonctions est une fonction générale 
dont les éléments extension et rotation ne se déduisent pas 
d'une façon simple de ceux de ses facteurs. 

Dans le cas particulier où les facteurs sont tous orthogonaux, 
les rotations successives qu'ils représentent se composeront en 
une rotation unique représentée par leur produit. 

Dans le cas oij tous les facteurs sont de pure extension, 
chacun d'eux peut's'écrire 

\ki=-[^]\i! ou \^!Wr'. 

On pourra donc, par des permutations convenables, associer 
toutes les fonctions canoniques et, 'd'autre part, toutes les 
fonctions ortliogonales, et par suite ramener leur ensemble au 
produit d'une fonction canonique par une fonction orthogonale. 
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Les axes de l'ellipsoïde résultant s'obtiendront à l'aide 'de 
cette remarque que si trois fonctions canoniques satisfont à 
l'équation 

\ll = \Si \(i, 

on a entre les axes des ellipsoïdes correspondants des relations 
de ta forme 

C'=:A.B. 



lauyGoOgle 



CHAPITRE I! 
Fonctions rapportées à. des axes obliques. 



§ 1. — Relations entre les coefficients d'une fonction oblique. 
Contraordonnées. 

Lorsque les fonctions vectorielles représentent des phéno- 
mènes relatifs à des milieux anisotropes de symétrie oblique, 
il est évident qu'il y a avantage k rapporter les fonctions aux 
axes de symétrie oblique du milieu. La fonction orthogonale 
précédente est alors remplacée par une autre que j'appellerai 
oblique, et dont les éléments sont assujettis à des relations que 
je vais indiquer, en ayant recours aux symboles vectoriels. On 
trouvera la forme explicite de quelques-unes de ces relations 
dans un travail intitulé : Des Constantes d'élasticité dans les 
milieux anisotropes, et publié dans les Mémoires de la Société 
scientifique de Bordeaux, 1886. Pour faciliter la correspondance 
entre les formules des deux textes, j'ai conservé la même 
signification aux lettres p, o, m, ^ et a. 

Les coordonnées d'un même point x„„ ï,,, par rapport à 
deux systèmes d'axes obliques, sont liées par les fonctions 
vectorielles 

(24) x—\p\i et ^ = \rs\x ou \p[~'x, 

que Ton obtient en appliquant le théorème relatif à la projection 
d'un contour. 

Les coefficients p et o sont proportionnels aux cosinus que 
font entre eux les axes de coordonnées d'un système avec les 
normales aux plans de coordonnées de l'autre système. 
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L'iMimination de x ou ? cuire ce» groupes d'équations 
conduit aux dix-huit relations 



6P,-t,^1, 6P,x; = i, 

6p;t:j =0, ôp.icj-o, 



i, j = 123, 



dans iesf];ielles P, P', x, t.' sont les vecteurs ayant pour consti- 
tuants les termes des lignes horizontales ou verticales de \p\ 
et |n|, d'après la convention du § 2 (ch.I). Ces relations, dont 
neuf seulement sont distinctes, montrent que les dix-huit coel- 
ficients p et o dépendent de neuf paramètres, qui sont les six 
angles des deux trièdres de coordonnées et les trois quantités 
définissant leurs positions relatives. 

Rankine a appelé contraordonnées d'un vecteur |„„ ^,„ ses 
projections orthogonales sur trois axes obliques. Par l'applica- 
tion du théorème des projections, on trouve que ces contraor- 
données sont liées aux ordonnées par la l'onction vectorielle 

1=1 1 cos 5, 5, cos 5,5, I 5 
(25) cos 5, 5, 1 cos 5, 5, 

I eos^a^i cos5)5, 1 11 
que je représenterai par le symbole 

On aura de même, en substituant des x aux |, 

X = \ml X. 

Dans ces fonctions, les lettres m„ ..., affectées d'indices égaux 
entre eux, représentent l'unité. Celles affectées d'indices iné- 
gaux, m„ ..., représentent les cosinus des angles des axes de 
coordonnées, que je désignerai aussi par m„ m„ m,. 

A ces fonctions, j'en adjoins une autre dépendant des angles 
que les arêtes d'un trièdre font avec ceux de l'autre. Elle est 

I cos 3*1 Ç, C0S3T, îi COSÎCjîj I 

cosiB,;, cosa;,5, cosic,5j 
I cosa;,5, cosaîjï, cosa;,?, j. 
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Je la désignerai par \a\ en posant 

Les déterminants de ces diverses fonctions satisfont aux 
relations suivantes : 

D^Df, = Dî, DoDu = D,D„ = D„, 

d'où on déduit 

D,Dw = i. 

D'autre part, si A^ et A, sont les vecteurs formés avec les 
lignes verticales ou liorizontales de rang i de la fonction \a\, 
l'application du théorème des projections donne 
(25j) a; = \ml P.-, Ai = \n/ ir;. 

La substitution de ces valeurs dans la fonction \a\ la trans- 
forme en 

(25,) |«|==Wip|,. |a| = W(W, 

d'oii l'on déduit 

(25,) |p] = M-|fll:z=\«l-^V/, H=W-'\o| = |onin/- 

Ces formules permettent d'exprimer les contraordonnées 1 
et X en fonctions les unes des autres ; on a, en effet, 

(25,) r=\m/|pi5=\m/\m/-»|(.U=|a|;=|«l\ix/-' = = \or;, 
(25,) ; = \|*/ Hx=\i>.l \i>.j-'\a\x=\a\x-\a\\ml-'x=\p\x. 

Elles différent des relations entre a; et S en ce que l'on 
remplace les fonctions \p\ et |n| par les transverses de leurs 
inverses, et pourraient, d'ailleurs, se déduire directement du 
théorème des projections. 

Les vecteurs sextinômes p, o, m et ja ont entre eux des 
relations indépendantes du vecteur a et que l'on peut écrire 
de bien des façons. Éliminons, par exemple, \a[ ou \a\ entre 
deux des quatre équations (25J, (25^), et nous trouvons, 
suivant le mode d'opérer, 

\m/ = Wl \î*/ |o|, \v.! = \p\\fnl\p\, 
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OU encore 

W-' = i"! w-wi, \»;-' = ipi w-'\pi 

Or, si on développe les calculs relatifs, par exemple, à la 
première de ces égalités, on trouve que les m sont exprimés 
en fonction des jj. par une fonction du type (13,); on est donc 
ramené à écrire simplement 

m =iW|lii, V. ==\\p\\m, 

(j.-'— lia 11 m-', m-» = ||p]|i*~'. 
On vérifie d'ailleurs facilement les deux premières formes, si 
dans l'égalité 

»w„a;î+ ... +2m,ja?,a!s+ . .. =vi„?î -h .■■ 2|i„ÇjÇ, + ... 

on élimine les x ou les 5 à l'aide des relations (24) et si on 
identifie les deux membres. 

§ 2, — Chatigement d'axes. — Invai'iants. 
Le problème consiste, comme au § 6, à transformer 

x' = \b\x 

en la fonction 

par les substitutions 

a;"=|p|r, x= IpI;. 

Mais ta fonction |pj n'étant plus orthogonale, on ne trouve 
plus entre les quantités d, g et r et leurs transformées les 
mômes relations que celles énoncées dans le § 6 (13). 

En conservant les notations déjà employées, cherchons donc 
tes conditions que doivent remplir les substitutions générales 

x' — \q]-'^' et X — |«|; 
pour que l'on retrouve entre les quantités b et ^ les reltt- 
tions (13). Or, pour que la rotation p exprimée en fonction 
des btj à l'aide du tableau (1â) ne soit fonction que de la seule 
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p=\,|-'r=|.|-'r, 


'.= llïll' =i\»ll'- 


r=ls|p =\î|f. 


l =||,l-'|),= |\s|-'|X 



4Ï B. ÉLIK. 

l'otalioii Vf- il l'aut ([ue les termes des colonnes du détermi- 
nant de \q\ soient égaux aux termes des lignes correspondantes 
du déterminant |s|, autrement dit, que les fonctions \s\ et j^j 
soient transverses l'une de l'autre; les équations {^^3„) se 
trouvent alors vérifiées. 

Si donc \q\ = \s|, on en déduit, par les mêmes procédés 
qu'au § Q, 

(26) |p|^\sl|6||«i, 

et par suite, 

(26J 

Tel n'est pas le cas actuel, puisqu'on opère avec une seule 
fonction jp|. Pourtant, si on transforme simultanément les 
variables x par la fonction |p|, et les coefRcients (t par certaines 
fonctions, je vais montrer que les nouveaux coflicients 3 se 
transforment d'après les fonctions (â6„). 

Posons, en effet, les relations 

î; = \m/B;..., i= 1,2,3, 

dont l'une, {t = 1), s'écrit explicitement 

^ b„ = b„ + m,bn -h Wjb,,, 
(27} \hi = ^hb,i + à., -^W(,fc,„ 

' lu, = m, &i, + m, *„ -I- 631, 

et d'autres semblables pour les coefficients 3, mais où la 
fonction \yi/ remplace \m/, 

j}' = \p\^', ir=:|p[=. 
On en déduit 

\p\~'\b\\p\ = \^\, 
ou, puisque alors 

I6| = \i»;-'|f|, ip| = w-'ill, 
]p\-'\mr'\t\\f] = vt'm, 



lauyGoogle 



FONCTION VBCTOHIELLE. 4& 

et, d'après les relations (25,), 

■ |a|-M5|lpl = W-'|fl ou \p\\b\\p\ = m. 

It en résulte les équations (â6„), dans lesquelles on fait s = î>. 
On arriverait à la même réduction, mais sous une Torme un 
peu différente, si on posait 

B,= \m/B.., 1 — 1,2,3, 

ou, explicitement pour i ^ \, 

!h„ = fc,, + «a bi, -4- Wj b,„ 
6,, = in,fc,, + m^b^i + (",,; 
de ces relations on déduit, en effet, 

\b\ = \b\\ml, |?l-=lè|V/- 
Par suite, 

\P\-'\b['\P\ = \^\ 
devient 

Mlè!\m/[pI-Wl|611al=|élW ou kl|ft|W[=iPJ; 

d'où, les relations (26a), si on y fait {s| = \cj|. 

Ce qui vient d'tUre dit relativement à la transformation de 
la fonction a.'' ^=\b\x S'applique, d'après les développements 

du S 6, i\ la transformation des dérivées -;;— et -r| . Si Ton veut 
" ' âx d^ 

que les rotations ... qui leur sont relatives se transforment 

par les formules (-13) et (13,), il faudra opérer simultanément 

sur X avec la fonction |p], et sur tt et s avec les fonctions 

(29) n = \m!u et ? = \i./?. 

Ceci revient évidemment à faire un changement d'axer do 
coordonnées et à substituer aux projections obliques u et o du 
déplacement ses projections orthogonales sur les axes auxquels 
on le rapporte et que nous appelons les contraordonnées. 

On peut résumer ainsi les résultats obtenus : 
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Étant données les égalités 

x' = \b\x a 5' = isiç, 

si les variables x et x' , 5 et Ç' sont liées entre elles comme le 
sont les coordonnées d'un même point rapporté à deux sys- 
tèmes d'axes obliques, c'est^-dire vérifient les équations 

si, de plus, on substitue aux coefficients b et p d'autres 
coefficients 6 ou & et g ou g liés aux précédents par les égalités 

B,' = \ro/BJ ou B,= WB,, ,-^^ g 3 
b; = \ml b; ou Bt=zi \ii/ B;, 

Xii-I et \m/ étant des fonctions vectorielles du type (25), les 
coefficients & et g vérifieront une fonction vectorielle à neuf 
termes, qui sera, dans le premier cas, 

l = \p\ft>\\PU 
et dans le deuxième, 

^ = H|fc]\o|, 
en posant 

N = |p|-'. 
Formons avec les demi-différences et les demi-sommes des 
coefTicients 6 ou p les quantités que j'ai appelées rotations, ... 
et que j'ai représentées par les lettres r, d, g et l affectées 
d'une barre ou d'un point lorsqu'il y aura lieu de les dis- 
tinguer; on aura 

p = |a|f, l^Wnill 
ou, par inversion, 

f=|p|p, Ï::=|WM>-. 

et de même 

p = \p\r, X=||n|/. 
De même que pour les coordonnées rectangulaires, il sera 
, dans le cas des coordonnées obliques, de former des 
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invariants. J'appellerai ainsi des fonctions composées avec les 
quantités r, d, g et m, et telles que leurs valeurs et leurs 
formes ne chaînent pas, lorsqu'on change d'axes de coor- 
données. 

Ces invariants se déduisent des équations cubiques dont les 
racines sont les diamètres conjugués des surfaces relatives à 
la fonction jh|. Je me bornerai îi écrire les seuls que j'utiliserai 
par la suite. Ce sont les analogues de ceux qui ont été désignés 
dans le cas des coordonnées rectangulaires par R, !„ I„ I„ 
et leurs formes sont bien connues dans le cas de la substitu- 
tion (27); ils sont : 

iK = rî + rî+rî+2iB,rjrs + 2m,r,r, + 2mjr,ri, 

3, = «iTi'rf, +«J7,'rf, + mf,'rf, +2m7,'sr,+2m7,'j, + 2»i7,'g'i, 

3» = ?î — à^â, + g', — rfj»/, + gl — rfj rf, 

+ îm» (ffiffa— 9 A) + 2m„ (>,?,— fl',rf,) + 2m„(ff ,5-, — 3,^,), 

m ' s\ l g\ 

I 9i 9i ^3 I 
Si l'on opérait avec la substitution (28), il faudrait changer 
les coefficients de m en leurs inverses. 

S 3. — Forums canoniques. 

J'ai appelé forme canonique d'une fonction générale en 
coordonnées rectangulaires celle où les coefficients d'indices 
différents sont égaux et de signes contraires. Je vais montrer 
qu'on peut simplifier encore cette forme et annuler ces coetïi- 
cients lorsque la fonction est rapportée à des axes obliques. 

Reprenons la fonction 

x' = \b\x, 
transformée en celle-ci : 

par les substitutions 
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Il en résulte : 

C'est l'équation {ii„) dans laquelle on pose 
|.| = lpl, « = o, 
DU bien encore, l'ensemble d'équations représentées dans le 
tableau (12). On a, sous forme explicite, pour l'un des coeffi- 
cients, 3„ par exemple, 

, K = {b^lP^^ + 6|,Pji + *itP,>)ci„ 
;32) \ ■+ (fi„p,. + 6„Pi, + &«p,,)t3„ 

+ (fraiPu + *jiPii + ''aiPM)f'»s; 
en général. 

Or, il est facile de voir que l'on peut, par un choix d'axes 
convenables, annuler tous les coefficients dont les deux indices 
5ont ditîérents. Multiplions, en effet, respectivement parc,,, c„ 

:es équations 

p„ — et &„ T= 0, 

ît retranchons les résultats. Opérons de même sur les autres 
;[uantités p, et nous aurons 

P» ~ P» 

~ Pi. 

3U 

6b,p._6b.p;_ 0b.p, _ 
p^, ~ p» ~ p» ^ ' 

jt deux autres groupes d'équations semblables. On tirera de 
les équations les neuf valeurs de p,, par la résolution d'équa- 
;ions du troisième degré en S. 

L'élimination des quantités p des équations (.32) par la 
nême méthode conduirait aux valeurs des ts à l'aide d'équations 
léduites des précédentes par la permutation des indices des 
ettres g. L'équation cubique resterait la même; ces valeurs 
le es sont d'ailleurs les inverses de celles de p. 
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Les trois directions ainsi définies et donnant à la fonction 
vectorielle sa forme la plus simple, sont les trois diamètres 
conjugués communs à la quadrique et à l'ellipsoïde de la fonc- 
tion |6]. 

Supposons, en effet, l'ellipsoïde de déformation rapporté à 
trois diamètres conjugués et ayant pour é<|uation 

Miî wî wl a. /'"V . 



})' 



C'est le lieu des points primitivement situés sur l'ellip- 
soïde vî 4- vî + vî =; 1, lorsque les coordonnées de ces points 
deviennent 

Formons l'équation de la quadrique directe définie pap la 
condition que le scalaire des vecteurs v et lu soit constant : 
Sviu = i. Si nous éliminons v entre les équations précédentes 
et ce scalaire, nous obtenons 



'(7^) = 



J. -" -'--'• 



Ces surfaces définissent une fonction vectorielle générale, 
dont les neuf coefficients seraient donnés par les équations (32), 
dans lesquelles on remplacerait les b ù indices égaux par 
P.i ?i. Pj. 6t les autres par zéro. Elles sont d'ailleurs rapportées 
à un même système de diamètres conjugués. 

On peut donc substituer à une extension suivant trois axes 
rectangulaires suivie d'une rotation une triple extension 
suivant trois axes obliques. 

Il y a cependant cette remarque à faire que la sphère, lieu 
primitif des points qui viennent, dans le cas d'axes rectangles, 
se placer sur un ellipsoïde, doit être remplacée par un ellip- 
soïde wî + vî -1- vj = i, possédant trois axes conjugués égaux 
dans la direction des axes de coordonnées! 
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DEUXIEME PARTIE 

APPLICATIONS 



CHAPITRE J 
Énergie de déformation. — Ëlasticllé. 



g 1. — État actuel de la théorie. 

Le calcul symbolique exposé dans la première partie va 
nous faciliter l'étude de quelques théopies de ta physique. Je 
commencerai par celle de l'élasticité, ou plutôt, des phéno- 
mènes statiques qui se présentent lorsqu'un milieu est déformé. 
Les théorèmes qui en résulteront trouveront en effet à chaque 
moment leur application dans les chapitres suivants. Mais 
auparavant, il est nécessaire de préciser quelques définitions, 
et il y aura avantage à exposer rapidement les hypothèses sur 
lesquelles s'appuie actuellement la théorie. 

Un milieu est dit hétérogène ou homogène relativement à 
une propriété déterminée suivant que cette propriété varie ou 
non avec la position du point oij on la considère. Dire qu'un 
milieu est homogène, c'est dire que la fonction vectorielle qui 
définit cette propriété en fonction d'un vecteur considéré 
comme variable a ses coefficients constants en tout point du 
milieu. Je ne m'occuperai, par la suite, que des milieux homo- 
gènes et des propriétés de nature vectorielle. 

Un milieu est dit isotrope relativement à une propriété 
lorsque l'ellipsoïde de la fonction vectorielle se réduit à une 
sphère. Sinon, il est anisotrope. 

Dans ce paragraphe et le suivant nous supposerons les axes 
rectangulaires. 

Les déformations l et r rapportées à deux systèmes d'axes 
différents sont liées par les équations (13) de la première partie. 
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Il en est de même des pressions résultant de ces déforma- 
tions, pressions que je désignerai par X ou S suivant qu'elles 
sont rapportées aux axes x ou ï. Soient, en effet, P,„ et r,,, les 
composantes suivant les axes a; et 5 de la pression s'exerçant 
sur un plan dont la normale a pour vecteur n ou i, suivant 
qu'on la rapporte aux x ou aux |. 

On sait, en effet, par le théorème du tétraèdre, que cette pres- 
sion est une fonction vectorielle de la direction de la normale 
dont les coefficients sont les neuf composantes des pressions 
s'exerçant sur les plans de coordonnées, et qu'on a par suite : 

P=|X|«, , = |S|,. 
D'autre part, les vecteurs P et w, n et v sont liés par les mêmes 
relations que celles qui existent entre les coordonnées d'un 
même point. On a donc 



d'oii 



= !2|v 



'-[p]-'m[ph 

c'est-à-dire les équations (-12). 

Op, on admet dans la théorie de l'élasticité que les forces 
extérieures qui équilibrent les pressions n'ont pas le caractère 
rotatoire, c'est-à-dire qu'elles ne peuvent équilibrer les couples 
5^.1 — X,.. Ceux-ci devant être nuls, les relations précédentes 
conduisent bien aux relations (13„). 

Les relations entre les pressions et les déformations ont été 
l'objet des recherches des physiciens depuis Poisson. Le pre- 
mier il admit entre ce» vecteurs une relation vectorielle 
sextinôme générale à trente-six coefficients, qu'il réduisit à 
quinze à l'aide d'hypothèses sur les forces moléculaires. Cette 
simplification revient à supposer symétrique la fonction 

X- ._ _ 

c,, 

(1) 
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adjoignant 

( C,( = C,j. C„ = C,„ C„ = C4, (') . 

îen admit simplement que la fonction était symétrique, 
la condition pour qu'après un cycle fermé de déformations 
vail élastique effectué soit nul. L'énergie étant le scalaire 
imi-produit des vecteurs X et l, 

2E = 0X( 

IV suite une fonction quadratique homogène de sis varia- 
à vingt et un coefficients réductibles à dix-huit par un 

d'axes convenables (î). Les pressions sont les dérivées 
tte fonction par rapport aux variables. 

nombre de ces coefficients est moindre pour les milieux 
i de plans de symétrie ou d'axes d'isotropie. 

des plans de coordonnées, a;,a;, par exemple, est dit de 
:rie si l'expression de l'énergie n'est pas altérée par le 
;ement de x^ en — x^ et de m, en — m„ c'est-à-dire par 
ingement de signe de g„ g^, r„ r,. 

plan bissecteur des plans de coordonnées, celui pas- 
par 0.-C,, par exemple, est un plan de symétrie si 
gie n'est pas altérée par la permutation de x, en a;, et «, 
, c'est-à-dire par la permutation de d, en d„ g, en g, et 

r„ r, en — r„ — r„ — r,. 

axe est d'isotropie d'ordre n si l'énergie n'est pas altérée 

ne rotation de — autour de cet axe. 

R 

conditions réduisent à neuf, six et trois les constantes 
ticité des systèmes cristallins à axes rectangles, et à 
six et cinq celles des systèmes à axes obliques. 



.. l, représentent ici les trois dilatations d„, et les doubles de deux glis- 

>g„.. 

iigton. Traniactioni of Irlsli Academy, lHi6. — Rankine: Tranaactions 
I Society, lt»5. — De Saint-Venant, DUlribtilion dea élatlicilés dans 
le anitotrope (Joumaj de mathématiquei pures cl appliquées, 1863). 
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Dans le cas des milieux isotropes, l'expression de V 
élastique devant être indépendante des axes devient 

AIÎ -t-4BI„ 

I, I, étant les invariants du | 7, première partie, A et B les 
carrés des vitesses de ppopagation des ondes longitudinales et 
transversales. 

Ces simplifications ne furent plus trouvées suffisantes, au 
moins pour d'autres milieux que la matière pondérable, lors- 
qu'on voulut rattacher les théories de la lumière et de 
l'électricité à celle de l'élasticité. 

Partant de cette idée que l'anisotrope de l'éther doit être 
distinguée de celle de la matière pondérable, et qu'elle doit 
être attribuée aux modifications que produisent, sur un milieu 
primitivement isotrope, des tractions opérées suivant trois 
directions rectangles, de Saint-Venant arriva aux six relations 

c.i + 2c., = Vc., c,, = = c,,, 
/ 
1 c„ + 2c„ = VCii c„ = ■= c„. 

Si l'on accepte ces conclusions, l'expression de l'énergie est 
. la somme des quantités 

a,a, (yî - rf.rf,) + a, a, {g\ - d.rf,) + a. oc, {gt ~ rf,d.) 

multipliées respectivement par les constantes A et 4B. 

On remarquera que ces quantités contiennent des constantes 
qui ne dépendent que de la symétrie du milieu, et qu'elles 
deviennent des invariants si le milieu est isotrope. Quant aux 
facteurs A et B, ils ne dépendent pas de la symétrie, mais de 
la nature du phénomène. 

Dans les phénomènes statiques d'élasticité concernant les 
milieux isotropes, ces constantes A et B ne sont pas indépen- 
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danfes et doivent satisfaire à l'inégalité 2A >B, afin que 
l'énergie ne puisse changer de signe pour des valeurs quel- 
conques des déformations. Cette condition, d'après M. Beltrami 
(Il nuouo cimento, 1886) empêche d'expliquer les actions 
électrostatiques par les déformations d'un milieu élastique 
conçu comme d'ordinaire. 

D'après Sir W. Thomson (Philos. Mag., 1888), lorsqu'on 
s'adresse à d'autres faits que ceux de l'élasticité statique, ces 
constantes A et B pourraient être indépendantes. 

Il y aurait évidemment intérêt à savoir si l'expression de 
certains phénomènes élastiques lumineux ou électriques est 
aussi, pour les miHeus anisotropes,ia somme de fonctions ne 
contenant que des constantes dépendant de la symétrie du 
milieu, multipliées par d'autres constantes qui ne dépendraient 
que de la nature du phénomène. La séparation des constantes 
en deux groupes ayant chacun une signification spéciale amè- 
nerait une réduction dans leur nombre. C'est cette hypothèse 
que je me propose de développer afin de déterminer à quels 
cas elle est applicable. 



§ 2. — Syiiiétne rectangle. 

Avant de former, à l'aide des fonctions vectorielles, les 
énergies spéciales E,, E„ ... dont la somme constituera l'éner- 
gie totale de déformation d'un milieu anisotrope à symétrie 
rectangle ou oblique, je vais d'abord, en me limitant au cas 
de la symétrie rectangle, les déduire de cette seule considéra- 
tion qu'elles doivent se réduire aux invariants If, I, et R*, 
lorsque le milieu devient isotrope. 

Par hypothèse, une quelconque de ces énergies E est une 
fonction homogène quadratique des neuf variables -p ou b, ou, 
ce qui revient au même, des neuf quantités 
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Les forces élastiques développées lorsque le milieu est 
déformé, sont les coefficients de 5u,„ dans la variation de 



'///"" 



étendue à tous les éléments dn du milieu. Une intégration par 
parties la décompose en une intégrale étendue au volume, dans 
laquelle les coefficients de 8«,„ sont, en désignant par X„ les 
neuf dérivées de E relativement aux variables 6, 



ôx. 


<)X„ 


àx^' 




àx. 




àX., 


âX„ 


à^.. 



. dXi àx^ dXj 

et en une autre étendue à la surface, dans laquelle les coeffi- 
cients de îu,„ sont, en désignant par n le vecteur unité normal 
à la surface, 

JX„II, 4-X„B, + Xi,«, 

(4) j X„n. + X„n, + X„n, ou IX»]. 

'x,.«. +X„», + X„n, 

Les quantités (3) sont les forces élastiques qui font équilibre 
aux forces d'inertie ou extérieures, forces qui peuvent être 
nulles. Les forces (4) équilibrent les forces appliquées aux 
limites du milieu ou aux surfaces de discontinuité. 

Les forces d'inertie ou provenant d'agents extérieurs peuvent 
avoir un caractère rotatoire ou non et, suivant ces cas, sont 
capables ou non d'équilibrer les couples dus aux pressions 
tangentielles et dont les composantes sont X„ — X,,, ... (^s 
couples seront nuls ou non et, par suite, la fonction X sera 
symétrique ou non suivant qu'on se trouvera dans un cas ou 
dans l'autre. 

Soit Es ou E, l'énergie* exprimée en fonction des variables 6 
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ou d, g, 


r;oii a, à 


après 


la délintlion de X 


pour ime 


quelconque 


de ces quantités 


x„ 


db,, 


àEg 
ai, '■ 


., 




(3) 


X. 


i)E6 
àEb 


àEg àg, 

àg. àb,, * 


àEg àr, 
àr, àb„ 


1 àEg 

~ï àg, 

1 i)E!) 
~2 àg. 


làEg 

i àr,' 

^MEj. 

a <(r. ■ 



A»i -r Aé, — ■ ^— i A>» -^ An — — ^T * 

àg, àr. 

Si donc la fonction est syinétpique, il faut que E, soit indé- 
pendant des rotations. 11 ne subsiste alors que sis tensions 
liées aux déformations l par une fonction sextinôme X = \c\cl. 
Si elle n'est pas symétrique, on sait qu'elle se décompose 
eu la somme d'une fonction symétrique et d'une fonction gauche 

^^ M. + Vl. t^ M.-\oL ,, 

Admettons cette dernière hypothèse, mais supposons le 
milieu symétrique relativement aux trois plans coordonnés. 
L'énergie totale peut s'écrire 

Cherchons les valeurs que doivent prendre les coefficients, 
pour que dans le cas d'isotropie, c'est-à-dire si les bissecteurs 
des plans de coordonnées sont des plans de symétrie, cette 
expression devienne l'un des invariants lî, I, ou R'. 

La condition de symétrie donne à 2E la forme 

/7i aF - \ ''"^''' "^ '*"' "^ "^'^ '^ 2''"(''.''. + Mi + rf.rf.) 

^' i + f..( 12?, 1' + |2!/,1' + |2y,[') + c,,.{rl + r\ + rj). 

Si on ajoute la condition que cette expression doit contenir 
I au carré, on voit que les trois fonctions cherchées seront 
nécessairement 

/2E' — («;rf.+a;rf.+a;(/js 

( 2E"= atatrl + a^aïT', + «ï'aï'rî. 
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Par suite, le double de l'énergie totale est 
2E = AE' +4BE' + CE", 
A, D, C étant des constantes indépendantes de la symétrie, 
a , a', a™ étant d'autres constantes liées à la symétrie du 
milieu. Dans le cas d'un phénomène dont l'énergie est la 
somme de plusieurs de ces fonctions, il se peut que les a 
dépendent les uns des autres; c'est ce qu'a admis de Saint- 
Venant en posant « = a'. 

Lorsque l'on suppose « = »'::= «', les équations qui expri- 
ment les tensions en fonction des déformations prennent une 
signification remarquable"; ceci revient en effet à remplacer 
u ei X respectivement par Va.u et — r. c'est-à-dire qu'il faut 

modifier les trois projections du déplacement suivant trois 
rapports différents Ka,, Ka„ Ka,, et altérer les unités de 
mesure suivant ces projections dans trois rapports inverses 



Dans la théorie de l'élasticité statique, la constante C est 
égale à zéro. Désignons par P' et P' les produits «Isiai, ^'c^okl, 
et formons la fonction vectorielle des tensions. On l'obtient 
en prenant les dérivées relatives à (Z et 2j de la quantité 
-(AE' + 4BE'). Ses coefficients sont inscrits dans le tableau 



qui suit : 



(9) 



; — 2B-, Aa,' 

' P' P' 

7-2B-, A^-21 









2?. 2?3 
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Ces coefficients, que nous désignerons par C;^, satisfont aux 
équations 

I c„ + 2c„=Kc„c„, 

et, si on pose a' = et', aux équations (2) proposées par de 
Saint- Venant. 

Ils peuvent alors se partager en trois groupes, dont deux ont 
leurs termes ppoportionnels entre eux et inversement propor- 
tionnels aux racines carrées des termes du troisième, comme 
t'indique la suite des rapports 

(H) c,, : c.. : c,. = c., : e,, : c,, = — = : : 

Kc.. |/^. K^ 

Si a' diffère de a', mais si on supprime les indices de ces 
lettres, ce qui revient à supposer que les bissecteurs des 
dièdres coordonnés sont des plans de symétrie, l'énergie ne 
contient plus que deux invariants et deux constantes, — et - 1 
et le milieu doit, par définition, être appelé isotrope. Un cristal 
cubique ne différerait pas d'une substance non cristallisée au 
point de vue élastique, bien que la propriété du clivage les 
rende si distincts. Il se pourrait cependant que le clivage fût 
complètement indépendant des phénomènes d'élasticité, abso- 
lument comme la ténacité d'un corps non cristallisé est indé- 
pendante de ses constantes d'élasticité. 

Les équations (10) établissent à elles seules la différence 
entre la théorie de Green et celle fondée sur les invariants, au 
moins pour les milieux à trois plans de symétrie. La valeur (6) 
de E peut en effet s'écrire, lorsque les rotations disparaissent, 

{c„d, + c„d, + c„rf,)' + 4Cm(?Î — <^srfi) + iCisC^Î — <*.<it) 
+ ic„(5Î — d,rf,) 
+ 2(c„ + 2c„ — l^c,",c«) "^t*^) + 2 {c„ + 2c„ — Vc^^) d,rf, 
+ 2 (c„ + 2c„ — yc,iC„) d.rfj, 
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quantité formée des fonctions E', E" et d'une troisième qui 
s'annule lorsqu'on admet les équations (iO). 

§ 3. — Déformations et pressions dans un milieu 
à symétrie oblique. 

Je vais maintenant examiner sî l'emploi des fonctions vecto- 
rielles conduit à des résultats plus généraux que ceux trouvés 
précédemment. Il y a lieu de diviser la question en deux 
parties. Dans ce premier paragraphe, je chercherai, en partant 
de l'expression de l'énergie adoptée par Green, si l'emploi des 
coordonnées obliques, associé à la notion de symétrie, peut 
conduire à des simplifications dans le nombre des constantes 
d'élasticité, et à préciser celles que l'on doit adopter pour 
les milieux cristallins obliques. 

Dans le paragraphe suivant, j'étendrai à des milieux obliques 
l'hypothèse indiquée dans le paragraphe précédent d'énergies 
distinctes déduites des invariants que l'on rencontre dans les 
milieux isotropes. 

Soient, comme dans la première partie, u et 9 les vecteurs 
représentant le déplacement d'un même point rapporté à deux 
systèmes d'axes obliques x et Ç. On a vu que ces quantités 
sont liées par la même fonction oblique [p|, c'est-à-dire que 
l'on a 

x=\p\% et u=\p\r, 

que les dérivées - - 

(H) lî'|-'i6||Pl = ifli 

que, si l'on substitue aux projections obliques u et ■ 

déplacement leurs projections orthogonales définies par 

û = \mlH et ç =; \\>-h, 
l'égalité précédente devient 

(Hj \p\\s\\p\ = m 
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d'où Ton déduit, entre les déformations î »■ et X p, les relations 

(H.) ï=l\p|-'|Â. f=\p\ê- 

Désignons par E ou e l'énergie de déformation, suivant 
qu'elle est exprimée en fonction des h ou des ^, par T ou 3 
les dérivées de ces expressions en b ou p, c'est-à-dire par 
définition, les tensions. 11 reste à voir d'après quel type de 
fonctions vectorielles se transforment ces tensions et quelle 
est leur signification mécanique. 

Supposons que la fonction E ne contienne pas les rotations 
comme en élasticité statique, et ne soit fonction que des six 
déformations /. On a 

^Xi~^l^^\,'^ dl, âh'^ dl,d\,'^ ... (_l ... (i, 

ou 

S;=X,CT(1 + X,n,î + X,na+ ... = ||o||X, 
puisque les dérivées des énergies ne sont autres que les 
tensions et que l et À vérifient l'équation (iU). La résolution 
de ces équations par rapport à X conduit à l'équation 

(18) X = I iPl IS, 

c'est-à-dire à une fonction à la fois transverse et inverse de 

celle liant les déformations letX 

Dans le cas où l'énergie contient les rotations, il existe neuf 
tensions distinctes qui sont les dérivées partielles de l'énergie 
par rapport aux neuf quantités b. 

Par la même méthode que précédemment, on trouve que 3 
etX sont liés par une fonction vectorielle à neuf termes dont 
les coefficients -^ s'obtiennent à l'aide de l'équation (11) : 
(lï.) ' Ë-[pllXl\p|. 

Mais cette dernière conduit à l'équation (12) à la condition 
d'y remplacer les tensions par les demi-sommes — ' ■ ■- " des 
tensions symétriques par rapport à la diagonale de leur déter- 
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minant. Si, de plus, on désigne par M, ^= X„ — T„ îes couples 
formés avec les différencea des tensions, on déduit de (12,) la 
fonction vectorielle à trois termes 

(12.) M = \p|-'n;. 

La relation (12) peut être obtenue à l'aide de considérations 
purement tirées de la mécanique. Cette méthode, outre qu'elle 
servira de vérification à la précédente, mettra en évidence le 
sens qu'il faut attribuer aux tensions X. 

Désignons par S et X les projections obliques sur les axes 
des tensions s'exerçant sur les plans de coordonnées; 

Par Xç, la projection oblique sur Ç, de la tension s'exerçant 
sur le plan x, x,; 

Par Al une aire parallèle à l'un des plans x interceptant sur 
les plans des Ç des aires Aç,„; 

Par x' et 5' les normales aux plans des axes a: et Ç; 

Par le théorème du tétraèdre appliqué i\ des coordonnées 
obliques, on a les relations 

K,Xu - Ae,S„ 4- Ae,S„ + Aî,H,i, i = 1, 2, 3. 
De plus, les projections des aires A,, et Aj, sur 5', étant égales, 
on a aussi 

Aj., cos 5(3;! = k\f cos^i %'i . 



La définition dp- 


,uanlités p 


^ donne 




Pijcosa;, 


3?; = cos 5(3;;. 


et par suite 








A,, cos a;, 3^; 


p„=iA5,cos=;! 


Posons alors 






2; 

COS^ 







égalités qui définissent les quantités X et S au point de vue 
géométrique, et éliminons les aires A des équations précédentes ; 
le n'fsultat est 

Xïi — cosa;,a;;(pi,¥,i + p„3,j + p.,2„), i — l, 2, 3. 
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D'autre part, la tension exercée sur un élément de surface se 
comporte comme un vecteur, et on peut lui appliquer les 
formules de transformation relative aux vecteurs, ce qui 
permet d'écrire 

X,i = p..Xï, +p..Xt, -HprtXE,, 1 — 1,2, 3. 

Il suffit d'éliminer les Xe entre les deux groupes d'équations 
qui précèdent en introduisant X pour obtenir la relation 
cherchée entre X et 3. 

Si on admet le théorème de réciprocité relativement aux 
tensions X, c'est-à-dire si 

ou 

X,j=X,o 
la relation entre X et S prend la forme (12) déjà établie 
T=\\p\\Z. 
On arrive à des conséquences analogues si l'on considère 
'énergie comme une fonction non p 

du 

mais d autres que j écrirai —r = b. 
dx 
Posons, \a\ étant la fonction définie au § 1, chapitre II, pre- 
mière partie, 

«— |p|ç, a:— IpIS, x=\a\x, i=\a\\. 
Le calcul des dérivées conduit entre h et la quantité— *(>*© 
j'écrirai g à la relation 

qui ne diffère de {H„) que par le changement de \ip\ en \cj|. 
On trouvera de même 

i^|[p]|x. f = wip. 

Cette autre forme de l'énergie pourra aussi être considérée 
comme le scalaire des variables l et de quantités X qui seront 
les tensions. Ces tensions que l'on désignera par X ou S, 
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suivant qu'elles seront rapportées aux axes x et £, sont reliées 
par l'équation vectorielle 

X = l\a||É. 

Elles ont une signitication analogue i^ celle des tensions X, 
Pour former ces dernières, il a fallu considérer le trièdre des 
x' conjugué du trièdre des x, c'est-à-dire formé avec les 
normales aux plans des coordonnées x, et projeter orthogona- 
lement sur chaque normale les composantes obliques de la ten- 
sion s'exerçant sur le plan correspondant des coordonnées x' . 
Au contraire, les tensions X senties projections orthogonales 
sur les axes x des composantes obliques de la tension s'exer- 
çant sur les faces du trièdre x' . 

Les déformations et les tensions rapportées à des axes 
obliques se transforment donc, lorsqu'on passe d'un système 
d'axes à un autre, absolument comme lorsqu'on opère avec des 
axes rectangles. Nous pouvons donc adopter pour expression 
de l'énergie une forme quadratique des déformations à vingt-un 
coeificients, pareille à celle de Green, mais où les déformations 
sont prises dans le sens indiqué, c'est-à-dire fournies avec les 
projections orthogonales du déplacement sur les axes obliques. 
Cette expression une fois déterminée par des raisons de symé- 
trie pour un milieu donné, on devra, pour les recherches 
expérimentalej à l'aide de la relation (Ha), revenir ù des 
coordonnées rectangles. On trouvera à la fin du paragraphe 
suivant un exemple de cette transformation. 

Pour le moment, je fais remarquer que ce procédé, non 
seulement diminue dans bien des cas le nombre des constantes 
d'élasticité, mais fournit un choix bien plus considérable pour 
les systèmes de constantes à adopter. Le classement des 
systèmes cristallins peut se faire en eifetjde la façon suivante 
quant au nombre des constantes d'élasticité. 

Les milieux ayant trois plans de symétrie rectangle ou 
oblique dépendent de neuf constantes. Dans le cas des milieux 
obliques, lorsqu'on sera revenu à des coordonnées rectangles. 
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un se trouvera en présence de vingt et un coefficients, mais 
dépendant de neuf papamètres et des angles des axes obliques 
antre eus et avec les nouveaux axes rectangulaires adoptés. 

Ces neuf constantes se réduisent îi six ou à trois si un ou 
deux des bissecteurs des plans coordonnés sont des plans de 
symétrie. La possibilité de cette réduction est basée consé- 
quemment sur l'application à des coordonnées obliques de la 
définition du plan de symétrie donnée en coordonnées reo- 



II n'y a lieu de considérer un nombre de constantes supérieur 
\ neuf que si le milieu est hémiédrique, et il suffit de chercher 
combien de plans de symétrie sont supprimés par les facettes 
liémiédriques. C'est ainsi que certains cristaux cubiques hémié- 
iriques exigeront vingt et une constantes en axes rectangulaires 
[tour leur détermination ; pour d'autres, pour les tétraèdres par 
exemple, il suffira de trois constantes, à la condition de rap- 
Dorter ces tétraèdres à des axes obliques. 

§ 4. — Énergie élastique dans les milieux à symétrie oblique. 
Discussion, 

Revenons maintenant à l'hypothèse du § 2, consistant îi 
'egarder l'énergie des milieux isotropes comme la somme 
i'invariants. Si, dans ces milieux, les forces dépendent des 
Jérivées -p, la généralisation la plus simple de cette idée sera 
le substituer il « et à a? des fonctions vectorielles générales de 
îes quantités, puis de remplacer dans les invariants I,, I,, R', 
es valeurs que prennent les déformations j- après cette subs- 
titution. On obtiendra ainsi des expressions qui devront repré- 
ienter dans les milieux anisotropes les énergies spéciales à 
îertains phénomènes attribuables aux déformations du milieu. 

Posons 

o=|?]tt et |— jsl-'iF, 
q\ et [s] étant des fonctions générales, et cherchons ce que 
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4s 
deviennent les quantités -t| après cette double substitution, 

puis les rotations... formées avec ces dérivées. Le calcul se fera 

à l'aide du tableau (12) de la première partie. Si nous voulons 

que ces relations, glissements... se transforment d'après les 

équations (lia), nous savons que les fonctions \q\ et |3| doivent 

être transverses l'une de l'autre. Si cette condition est réalisée, 

on 3, en désignant par ^ et 6 les dérivées ^ et j-^ 

I&I-UIK'II^I, 

dont on déduit, entre les relations et les glissements, 
(13) (=|H|-'|)„ .— I.Ip, nec \}1 = 1.1. 

il est entendu qu'on devra employer pour cbacun des inva- 
riants une paire distincte de fonctions \q\ et [s|. Mais, pour 
simplifier l'écriture, je conserverai les mêmes lettres pour les 
trois invariants, en les affectant de un, deux ou trois accents 
lorsqu'il sera nécessaire de les distinguer. 

Par la suhstitntion (1 1„), l'invariant linéaire 

I, = (/. + rfj + d, 

devient, en dési;3;nant par ©(Qi)"^', (Oj)~'i -.., les scalaires for- 
més avec les colonnes verticales du déterminant de \q\~\ 

A e ^ 5.â(0i)-'+ î,5(o,)-' ^- î.s (o;)-' 

' +2ï.S(Q;)-'(O.}-'+2ViS(O;)-'(Qi)-'+2T,0(Q;)-'(Q;)-'0. 
Pour le calcul de la transformée de 

Ij — jî — rfjrfj + gl — d^d, -\- gl — d,d,, 
on a d'abord, pour g\ — dtd, seul, 

(YÎ-3i3>)9ri + (ïî-B,3,)îî.+ ---+2(î,Y,--,v,09„V„ + -, 
et, en lui ajoutant deux expressions semblables, 

[ [(v;-3.î,)â'Oî-i-(Yî-B3î,)9'o:+(ïî-3,s,)Soî 

' ' +2(Î3Ï.-Y.V.)90,0,- 
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]{*= 1^ + rl + rl devient 
'SOî + pjSoi + fîSQi-H 

)tale sera, d'apn^-s notre liy[iolli(''se, 

z = \e +4Bs' + Cî'". 

is et les couples qu'elles produisent seront donnas 
es en Syp et dépendront des vinj^l-sept coefficients 
s les fonctions q', q' et g". 
us immédiatement dans l'hypothèse admise en 
ique, où les tensions sont symétriques, c'est-à-dire 
nsidère que les énergies e' et t, la première carac- 
mpressibilité du milieu, la deuxième sa rigidité. 
on de cette énergie totale revient à celle des fonc- 
, lesquelles contiennent dix-huit paramètres. Les 
emarquables qui se présentent sont, en commen- 
)Ius simples : 

ïs. — Les fonctions l^'j et \q'\ sont symétriques 

(première partie) qu'une fonction symétrique est 
me fonction canonique symétrique par une fonc- 
lale. On posera donc 



\ir 



K, K, K, 



^Mllp], 



constantes, pt il en pésultepa 

.il = K|l>,. K,P,. K.P., j = 
K,P„ K,p„ K.p„ 
K,p„ KjPi, K,p„ 
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Les valeurs depiL, ....s'obtiendront en annulant les scalaires 
fiQ,Oa, ■.-■ Les scalaires 0(01)'"', {Oi)"'> ■-, s'annuleront par 
cela même, résultat que l'on pouvait obtenir par Tapplication 
des formules (4) (première partie), ou encore en remarquant 
que les ellipsoïdes direct et inverse de \ql ont la même orien- 
tation. Or, annuler ces scalaires revient à rapporter ces ellip- 
soïdes à leurs axes principaux. 

Les énergies t et s' prennent alors les formes (8) du § 2, 
dans lesquelles on admet a' = a'. On retombe ainsi dans l'hy- 
pothèse de de Saint-Venant, considérée antérieuFCment, d'un 
système de coefficients d'élasticité ne dépendant que de trois 
paramètres et satisfaisant aux équations (2). 

Cette réduction, à la forme normale par un changement 
d'axes, revient à remplacer 

Deuxième cas. — Les fonctions 1^' | et \q'\ sont symétriques ■ 
distinctes, mais telles que leurs ellipsoïdes ont la même orien- 
tation. 

Les égalités précédentes sont applicables; la fonction ortho- 
gonale est la même dans les produits qu'on identifie aux 
fonctions q' et q', mais les facteurs canoniques diffèrent. Les 
coefficients d'élasticité sont ceux donnés par le tableau (9) et 
dépendent de six paramètres. 

Réduites à leurs formes normales, les déformations sont, 
suivant qu'elles mesurent la compressibilité ou la rigidité, 

(16) ^ ^ .... 



Tfoisième cas. — Les fonctions \q'\ et |}'| sont f 
mais ou égales ou inverses ou transverses l'une de l'autre. 

La discussion de ce cas va se trouver complètement contenue 
dans celle du cas suivant absolument général. 
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Quatrième cas. — Les fonctions jg'j et \q'\ sont générales 
et différentes. 

Formons avec les dérivé(>s de l'énergie relative à d et '2g les 
fonctions sextinômes représentant les tensions. Pour simpli- 
fier récriture, je représenterai par i„, ..., et Tju, ..., les coeffi- 
cients qui entrent, sous forme de scalaire, dans les expressions 
e' et e' (14). Ces quantités e et t^ sont les coefficients des ellip- 
soïdes inverses des fonctions \q'\~' et |g'|. 

Cette fonction sextinônie, qui lie les tensions aux variables d 
et 2y, est représentée par le tableau suivant, dont je n'écris 
que les termes situés d'un côté de la diagonale, puisqu'il est 
symétrique : 

: *'. <'. ft, 2Yi 2ïi ' 

■ S„ Aef, As„Ê„— 2Bi;„ Ae„E„— 2Bt;„ As„e„+Br;„ Aï„î„ A€„s, 

12„ . As», Ae„e„— 2Btî„ As,,E„ As„e„ + BY;„ Ae^Sj 

(19)j3„ . AeÎ, A£„s,3 Ae„s„ Aê,jE. 

Is.j . . . Ae|j +Bt;,, As„s„— Bt;„ Ae„£, 

'S„ . . . . AeÎ, +Bï;„ Ae„E, 

U„ . . - . ■ . Aeî, 

La réduction de ce système à sa forme normale, c'est-à-dire 
au moindre nombre de coefficients, est différente suivant qu'on 
utilise des axes rectangulaires ou obliques. • 

Dans le premier cas, le seul procédé de réduction est de 
rapporter à ses axes un des ellipsoïdes que je désignerai par (t) 
ou (ï;), suivant que ses coefficients sont les quantités t ou ij. 

Si c'est l'ellipsoïde de oonipressibilité (e) qui est rapporté à 
ses axes, on a 

Les coefficients C du tableau précédent restent au nombre de 
quinze, mais ils vérifient les équations (10), et de plus 

C, — — c„, c„ — — c„, c„ = — c,;. 

Si, au contraire, on rapporte l'ellipsoïde de rigidité (t;) à ses 
axes, on a 
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Les vingt et un coefficients diffèrent tous de zéro, mais ne 
dépendent cependant que de neuf paramètres. 

On aura le cas où le milieu a un plan de symétrie en écri- 
vant la condition pour que les ellipsoïdes aient un axe dont la 
direction soit la même. 

Pour réduire les coefficients ù six, il faudra que les ellip- 
soïdes (e) et (»i) aient même direction d'axes ; ceci arrivera si les 
fonctions \q'\ et |^'] étant différentes, l'ellipsoïde inverse de 
l'une est orienté comme l'ellipsoïde direct de la fransverse de 
l'autre. La discussion du troisième cas est par cela môme 
achevée. 

11 servirait peu de faire l'énumération des cas où l'un des 
ellipsoïdes est direct ou a des axes égaux à ceux de l'autre, etc. 
Je crois plus utile d'appeler l'attention sur les formes normales 
que peuvent prendre les énergies lorsqu'on les rapporte à des 
coordonnées obliques. 

Les deux fonctions \q'\ et \q'\ étant quelconques, on peut 
toujours, par un choix convenable d'axes obliques, ramener le 
système des coefficients d'élasticité au type (!)). Il suffit de 
rapporter les ellipsoïdes (t) et (tî) au système de diamètres 
conjugués qui leur sont communs. 

On a vu (première partie) que toute fonction vectorielle 
générale pouvait être considérée comme le produit d'une fonc- 
tion canonique symétrique par une fonction oblique, f! n'y a 
qu'une solution si l'on veut que l'ellipsoïde et la quadrique de 
la fonction aient leurs équations ramenées à la forme canoni- 
que, mais il y en a une infinité si l'ellipsoïde seul doit être 
ramené à cette forme. 

\p\ étant une fonction vectorielle reliant les axes rectangu- 
laires, auxquels sont rapportés les invariants primitifs, aux axes 
obliques, auxquels sont rapportés les nouveaux invariants, je 
poserai 

[g'\~'= K;o„ K,a„ K>„ —\jz\\i'l 



k;o„ 


K,a„ 


k;iii 


k;o„ 


K,o„ 


k;o 


Kjr^,, 


K>,. 


k;o 
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lî'i= K>.. k;p,, k;/.,, =\k'/\p1, 
k;p„ kip,, Kip,, 
KjPi. k;p„ k;p„ 

K' etK' étant des constantes, les fonctions |g'| et \q'\ seront 
ramenées à des fonctions symétriques et canoniques. 

Les six équations obtenues en annulant les scalaires des 
produits des termes de deux colonnes de |5'|~' etdedeuxlignes 
de \q'\ détermineront les angles des axes obliques entre eux 
et avec les axes primitifs. 

Il résulte de là que le système le plus général des coefficients 
d'élasticité en coordonnées obliques et avec l'hypothèse d'éner- 
gies déduites de formes invariantes peut toujours prendre la 
forme (9), ce qui implique d'ailleurs que le milieu a toujours 
trois plans de symétrie oblique. Ses neuf coefficients vérifient 
la formule (10) et ne sont fonctions que de six paramètres, 
plus des trois angles des axes obliques entre eux. 

Lorsqu'on réduit le système des coefficients îi sa forme nor- 
male, les substitutions \q'\ et |^'|, qui ont transformé les 
dérivées^, reviennent simplement à remplacer ces dérivées 



tt étant la projection orthogonale du déplacement sur les axes 
obliques x, et stmam six constantes, dont trois dépendent de la 
compressibilité et les trois autres de la rigidité du milieu. La 
substitution directe de ces quantités dans les invariants I,I,R' 
donnerait immédiatement les fonctions t' s' e'". 

On peut résumer le contenu de ces g 3 et 4 en les quelques 
lignes suivantes : 

i° Lorsque des cristaux ont des axes obliques, on peut 
rapporter les déformations, les tensions et les coefficients 
d'élasticité à des axes obliques, et ces quantités sont liées entre 
elles par des équations absolument semblables à cellesadmises 
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lorsqu'on opère avec dos axes rectangles. On obtient une 
réduction dans les nombres des coeHicients en appliquant aux 
cas d'axes obliques la définition de la symétrie adoptée dans le 
cas d'axes rectangles. Ce nombre est de neuf s'il y a trois plans 
de symétrie rectangle ou oblique. 

2" Un deuxième mode de réduction consiste à admettre que, 
dans les milieux anisotropes, l'énergie d'une déformation est 
la somme de diverses fonctions qui déterminent respectivement 
la compressibilité, la rigidité et la torsion du milieu. Ces fonc- 
tions sont des invariants dans le cas d'un milieu isotrope, et 
on les déduit de ces invariants, pour les milieux anisotropes, 
en rempla^îanl la coordonnée x et le déplacement u par les 
fonctions vectorielles générales 

D'après les résultats fournis par le calcul, cette hypothèse 
implique que le milieu a trois plans de symétrie rectangle ou 
oblique, et que les coefficients d'élasticité vérifient les trois 
relations (10) du § 2, dont la raison d'être se trouve ainsi expli- 
quée. Comme dans le premier mode de réduction, les vingt et 
un coefficients qui subsistent lorsque les axes des coordonnées 
sont quelconques, ne sont fonctions que de neuf paramètres. 
Mais ceux-ci ne dépendent que des six longueurs d'axes des 
ellipsoïdes (e) et {■[) et de trois angles ; d'où résulte la réduction 
à six paramètres lorsqu'on rapporte le milieu au système de 
diamètres conjugués conmiune aux deux ellipsoïdes. 

Quel que soit le mode de réduction auquel on s'adresse, il 
faudra toujours, pour la vérification expérimentale, connaître 
le système des coefficients rapporté à des coordonnées rectan- 
gles. La marche du calcul est la suivante : 

Soient \a/, et \aje les fonctions sextinômes formées avec ces 
coefficients lorsque le milieu est rapporté à des axes obliques 
ou rectangles, [p| la fonction qui lie ces systèmes de coordon- 
nées, X := I \p| 1 1 celle qui lie les déformations. Si on substitue, 
dans l'expression de l'énergie exprimée en fonction des a, les 
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leurs (le À, les coefficients des variables i dans l'expression 
insfopmée seront les coefficients d'élasticité en axes rec- 
igulaires. • 

Appliquons ce mode de calcul au cas d'un milieu ayant trois 
es obliques formant entre eux le même angle de cosinus égal 
t. et avec la diagonale principale un angle 9 vériflant la rela- 
m (cas du rliomboèdre) : 



Supposons, de plus, que les faces du rhomboèdre et les 
ssecteurs de ces faces soient des plans de symétrie ; on trouve 
le les six constantes d'élasticité doivent vérifier, outre la 
lation 2fflOT = o,i — a„ obtenue dans la théorie ordinaire, qui 
admet que la symétrie rectangle, les trois suivantes, déduites 
! l'hypothèse d'une symétrie oblique ; 

tang*f ^ 2 — ' — —^! 

(Tu tang'? — a,i lang ç — 2ff,j — 0, 
fl„tang*<?+ 2fli, langç — 4(7,, = 0{'). 

A ce premier mode de réduction, adjoignons le deuxième, qui 
insiste, comme on l'a vu, à admettre que les équations (10) 
tnt vérifiées quels que soient les axes rectangulaires ou obli- 
les, et nous obtenons 



(£^ 



2 ' 



Voyons maintenant si les faits expérimentaux confirment ou 
firnient les résultats théoriques obtenus précédemment. En 
! qui concerne les milieux isotropes, j'ai été conduit à laisser 
déterminé le rapport des vecteurs de propagation des ondes 
ngitudinales et transversales A et B, ce qui est d'accord avec 

s, Élie, Coiislaitles 
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les expériences faites jusqu'ici. Ces expériences montrent, en 
effet, que la valeur de ce rapport peut varier depuis zéro, 
comme dans le cas des liquides parfaits, jusqu'à environ j, 
comme dans le cas de certains verres. 

Les expériences sur les cristaux sont rares, et je ne connais 
que celles exécutées dans ces dernières années par M. Voigt, à 
Gœttingen, et par M. Pockels. 

Elles montrent que les coefficients d'élasticité relatifs aux 
corps anisotropes ne vérifient ni les équations de Poisson, ni 
celles de M. de Saint-Venant, ni les équations plu8généraies(10) 
du § â. 

L'inexactitude des équations de Poisson s'explique par le peu 
de généralité des hypothèses qui ont servi à les établir. Pour 
que les équations de de Saint-Venant soient vériiiées, il faut, 
dans notre théorie, admettre que les ellipsoïdes de compressi- 
bilité et de rigidité, non seulement aient la môme orientation, 
mais encore que les longueurs de leurs axes soient inversement 
proportionnelles, hypothèse qui semble par trop particulière. 

Quant aux équations (10), qui caractérisent notre deuxième 
mode de réduction, c'est-à-dire l'hypothèse de deux ellipsoïdes 
distincts définissant l'un la compressibilité, l'autre la rigidité 
du milieu, elles devraient toujours être vérifiées; et cela quelle 
que soit l'orientation des axes de coordonnées. 

L'allure entre les résultats expérimentaux et celle indiquée 
par les équations (10) est, il est vrai, ta même. Pourtant le 
désaccord est assez grand, dans certains cas, pour m'empêcher 
de reproduire aucune donnée numérique. 

De ce que le désaccord existe entre l'expérience et la théorie, 
fauHl considérer cette dernière comme inexacte? ou bien faut-il 
admettre des incorrections dues aux difficultés de l'expérience? 
Je ne le pense pas. Sans contester l'exactitude des expériences, 
on peut, en effet, contester le point de départ dont on est parti 
pour les exécuter. 

Hemarquons qu'un coefficient quelconque ne peut être obtenu 
indépeiidaiiiuient des autres. Il faut le déduire des relations 
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qu'Us oiil entre eux et avec les dé l'ornia tiens et les actions 
mécaniques exercées. Le chois de ces relations a été fait d'après 
ce principe, admis jusqu'ici sans discussion, que la symétrie 
cristallographique d'un cristal commandait sa symétrie élasti- 
que. Or, c'est là que je veux en venir, on n'est pas en droit 
d'affirmer a priori qu'il y a coïncidence ou même des rela- 
tions entre les axes cristal lograpliiques et ceux qui définissent 
sa rigidité et surtout sa compressibilité. La propagation des 
ondes lumineuses nous offre un exemple d'un cas où la symétrie 
d'un phénomène ne coïncide pas toujours avec la symétrie du 
cristal; de plus, la détermination des axes cristallographiques 
s'appuie principalement sur le mode de cliva^ du corps. Or, 
le clivage semble dépendre principalement des propriétés mal 
connues que l'on désigne sous les noms de ténacité, dureté,.,, 
et jusqu'ici aucune relation n'a été établie entre ces propriétés 
et l'élasticité du corps. 

On est conduit ainsi à supposer qu'il n'y a pas de dépendance 
immédiate entre les données de la cristallographie et de l'élas- 
ticité. Ces propriétés élastiques se rattachent, au contraire, 
intimement aux phénomènes vibratoires de l'acoustique, et il 
est probable que si l'on pouvait déterminer par l'expérience les 
surfaces d'ondes relatives aux vibrations longitudinales et 
transversales, on aurait en mains le moyen le plus sûr de fixer 
les valeurs des coefficients d'élasticité. 

Ces considérations me semblent suffisantes pour nous empè- 
cTier de rejeter d'une fa^on définitive une théorie qui apporte- 
rait avec elle une si notable simplification dans la recherche 
des coefficients d'élasticité. 
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CHAPITRE II 
Piézoéleclrlcilé et pyroélectrlclté. 



§ 1. — Hypothèse relative à la distribulion de l'électricité sur les 
molécules d'un cristal. Classification des milieux piézoélec- 
triques résultant de cette hypothèse. 

Les expériences les plus précises relatives à la piézoélectricité 
ont été faites par MM. J. et P. Curie. En 1880, ils ont montré 
qu'un cristal de tourmaline comprimé parallèlement à l'axe 
ternaire dégage des électricités égales et de signes contraires 
aux extrémités de cet axe; que les signes de cette électricité 
sont renversés si on le décomprime, et que la quantité dégagée 
est proportionnelle à la pression. La compression d'un cristal 
de quartz suivant son axe ternaire ne produit aucun dégage- 
ment d'électricité; mais soit qu'on le comprime suivant l'un 
de ses trois axes binaires, soit qu'on le décomprime perpen- 
diculairement à ses axes, on obtient un dégagement d'électricité 
de même signe, dans les deux cas, sur les faces perpendicu- 
laires à l'axe binaire. Si s est la surface qui dégage l'électricité, 
p la pression par unité de surface, la quantité dégagée est égale 
à Ks.p. K est une quantité constante dans les cas précités, 
mais qui en général dépend de la direction de la pression et 
do la direction de ia surface pressée lorsqu'on fait varier cette 
surface. Je renverrai au traité de cristaltographie de M. Mallard 
pour plus de détails sur cet ordre de faits et pour l'exposé des 
principes de cristallographie dont j'aurai à l'aire usage. 

Dans ce chapitre, je me propose la résolution du problème 
suivant : Étant donné un système dépressions quelconques 
agissant sur un milieu piézoélectrique, détei'miner en 
grandeur et direction le dégagement d'électricité qui en 
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uUe. Je m'appiiiepai, pour la solution du problème, sur 
te notion qu'un corps est composé de molécules dont la 
ne peut être précisée dans certains cas. Bien que ce genre 
considération puisse être écarté, je le conserverai, parce que 
tt le procédé qui m'a conduit à la solution de la question, 
jue, de plus, il a l'avantage de fixer a priori les valeurs des 
fBcients piézoélectriques, valeurs que l'expérience devra 
ifier. Le mot molécule ayant été employé dans des sens très 
srs, il est indispensable de déterminer le sens que je lut 
(ne ici. 

,es molécules seront pour nous des surfaces tracées autour 
n point, et dont les rayons vecteurs représenteront comment 
! propriété vectoriplle d'un corps est distribuée autour du 
nt. S'il s'agit de symétrie, par exemple, ces surfaces devront 
séder la symétrie du milieu; s'il s'agit d'une propriété 
laire, elles devront être les surfaces équipotentielles relatives 
iette propriété. Elles ne sont assujetties, dans un milieu 
nogène, qu'à être toutes égales et pareillement orientées. 11 
a à s'occuper ni de leur nombre, ni de leur grandeur, ni 
leur pénétration mutuelle; les centres de chacune d'elles 
lient un milieu continu, condition nécessaire pour que le 
^ul différentiel soit applicable à l'étude des phénomènes, 
'admettrai avec Sir \V. Thomson que tout cristal piézo- 
îtrique est, à l'état neutre, dans un état de polarisation; en 
Bndant par là que chaque molécule matérielle à l'état neutre 
associée à une certaine quantité de substance électrique ou 
Iber, que la distribution de cet éther à la surface est abso- 
lent la même que celle distribuée sur un corps conducteur 
dimensions finies semblables à la molécule, qu'enfin cette 
iribution est indépendante de celle des molécules voisines. 
mt à la forme de la molécule, elle n'est assujettie qu'à avoir 
ymétrie du milieu cristallin. 

les hypothèses conduisent aux conclusions suivantes : 
Supposons le milieu et, par suite, ses molécules soumis à 
! déformation générale consistant, comme on l'a vu, en trois 
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dilatations transformant une sphère tracée autour d'un point 
en un ellipsoïde, et considérons : 1" le centre de gravité du 
polyèdre supposé homogène ; 2" le centre électrique, c'est-à-dire 
celui obtenu en assimilant à de la matière pondérable les 
quantités d'électricité distribuées à la surface de la molécule» 
et en cherchant leur centre de gravité. 

Deux cas se présentent alors : i" ou la surface moléculaire a 
un centre géométrique de figure, et alors le centre de gravité 
et le centre électrique coïncident avant comme après la défor- 
mation eUipsoïdale, car les masses matérielles ou électriques 
situées aux extrémités d'un diamètre se sont modifiées dans le 
même rapport et sont restées symétriques; 2° ou la molécule 
n'est pas centrée, et alors tandis que le centre de gravité n'a 
pas changé, le centre électrique s'est déplacé relativement 
à lui. 

Pour rendre ces idées sensibles, il suffit de considérer un 
tétraèdre régulier recouvert d'une couche d'électricité en équi- 
hbre sur sa surface. Toute élongation du tétraèdre suivant une 
hauteur aura pour résultat de déplacer l'clectrlcité vers le 
sommet correspondant; c'est le déplacement relatif des deux 
centres que j'identifierai avec la quantité que Maxwell a appelé 
le déplacement électrique, et dont il a désigné les compo- 
santes par f, g, h. 

Au fond, c'est formuler d'une façon plus précise la règle 
posée par MM. Curie que l'électricité positive dégagée par la 
compression d'un cristal se porte vers le ciité le plus pointu. 

C'est, de plus, une application de ce théorème dû à Liou- 
ville (*) que le centre de gravité d'un ensemble de masses 
électriques est le même que le centre de gravité de la couche 
équipotentieile qu'on peut leur substituer. 

Avant d'entrer dans le détail des calculs, une dernière remar- 
que est nécessaire : c'est que la dyssymétrie cristallographique 

(') Reprint ofpapers on Elevti'oUatics. de Sir W. Tliomsoii, 
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n'implique pas nécessairement la dyssymétrie (électrique. Si, 
par exemple, les extrémités d'un axe d'un cristal sont pourvues 
de facettes constituant des angles polyèdres qui ne sont pas 
symétriques l'un de l'autre relativement à un plan perpendi- 
culaire à cet axe, mais si ces angles peuvent être rendus tels 
par une rotation de l'un d'eux autour de l'axe, nous admettrons 
que ces axes sont symétriques au point de vue électrique, 
puisque la quantité d'électricité distribuée sur les deux sera la 
même. 

Pour éviter des redites, je vais procéder immédiatement à 
rénumération des types de cristaux pouvant donner naissance 
à la piézoélectricité. Les cristaux centrés, et par suite les subs- 
tances isotropes, étant exclus, certaines formes hémiédriques 
seront seules ii considérer. 

Dans le système cubique, il ne reste alor-s que les formes 
dont la symétrie est représentée par les symboles cristal logra- 
pliiques qui suivent : 

3L* 4L' 6L* OC Or. OP (non observée), hémiéJrie holoaxe, 
(a) 3L* 4L' OL* OC Ot: 6 P (boracile), anlihémiédrie. 
{(7,) 3L' 4L' OL' OC Oie OP (clilorate de soude), lélarloédrie. 

La première est à rejeter, car d'après la remarque faite pré- 
cédemment, les sommets opposés, dyssymétriques au point de 
vue géométrique, peuvent le devenir par une rotation de l'un 
d'eux autour de la ligne qui les joint. Les deux autres formes 
doivent présenter quatre axes piézoélectriques suivant la hau- 
teur d'un tétraèdre régulier, c'est-ii-dire quatre axes dont les 
extrémités sont différentes au point de vue électrique. 

Dans le système sénaire, des trois hémiédries holoaxes 

A' 3L' 3L" OC OP OP', 
i* OL" OC 3P 3P', 
4' OL' OC OP, 

la première a ses sommets symétriques au point de vue élec- 
trique; la deuxième a six axes piézoélectriques distribués 
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autour d'un septième, Taxe du cristal, et la troisième les mêmes 
six axes alternat ivemeiit de sif;nos contraires. 
Les liémiédries 

A' 3L' OC :: 3P', 

i' OL' OC T., 

pp.^sentent trois nxes piézoélectriques espacés de 120" dans un 
même plan. Mais ces cinq formes n'ont pas été rencontrées. 
Dans le système phoniboédriqiie, on a 

(b) A» 3L' OC OP (quarl?.), liémiédrie lioloaxe, 

qui a trois axes à 120° l'un de l'autre dans im môme plan; 

(6,) i' OL* OC 3P (loiirmaline), anliliémiéJrie, 

et 

4' OL' 00 OP (non observée), lélarloédrie, 

qui ont trois axes rangés symétriquement autour d'un qua- 
trième : l'axe du cristal. 
Dans le système quadriquo, en rejetant 

4' 4L' OC Ot. OP{héiniédrielioloaxe), 

qui n'a pas de sommets électriques distincts, il reste les liémié- 
dries 

4^ OL' OC 4P ou OP (non observée), 

qui offre huit ou quatre axes piézoélectriques rangés autour 
d'un cinquième : l'axe du cristal, et 

(c) 4' 2L' OC 0- 2P ou OP (cuivre pyrileux), anlibémiédrie, 

qui a quatre axes dirigt'-s suivant les hauteurs d'un tétraèdre 
isocèle. 
Dans le système orthorhombique, on a les formes 

((/) 3L' OC OP (sel de seignelte), hémiédrie lioloaxe, 

qui a quatre axes dirigés comme les hauteur.* d'un tétraèdre 

isocèle, et 

(rfi) L' OC OP (calamine), antihémiédrie, 
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qui, outre quatre axes dirigt-s vers une même région, en pré- 
sente un cinquième : l'axe binaire du c?tstal. 
Dans le système binaire, les hémiédries 

, . ( V OC OP (sucre), hémiédrie tioloaxe, 

' ( OA' OC P (non observée), anlihémiédiie, 

ne peuvent avoir que deux axes piézoélectriques. 
Enfin, le système asymétrique 

CL OC OP (non observé) 

ne pourrait avoir qu'un seul axe. 

Ces résultats se déduisent de l'hypothèse admise et montrent 
qu'en général les axes piézoélectriques sont des directions 
rangées symétriquement autour d'une autre qui peut être elle- 
même un axe piézoélectrique. 

§ 2. — Fonctions vectorielles liant le déplacement électrique 
aux défornifitions ou aux pressions. 

Pour passer de ces vues générales ù l'établissement des 
équations qui déterminent le dégagement de l'électricité en 
fonction des pressions exercées dans le milieu, il ne faudrait 
rien moins que connaître, non seulement la distribution élec- 
trique à la surface de la molécule, mais encore la nouvelle 
distribution lorsque la molécule est déformée; mais on peut 
tourner cette dilïiculté par l'artifice suivant : 

Ne laissons subsister du cristal constdf^pé que les directions 
des axes piézométriques, dont l'intersection commune sera 
prise pour origine. Portons sur ces droites, à partir de l'origine, 
des longueurs égales, et imaginons aux extrémités de ces lon- 
gueurs des masses électriques que l'on déterminera par des 
raisons de symétrie. Ces masses devront remplacer la distri- 
bution électrique sur la surface de la molécule non déformée. 
J'appellerai pôles les extrémités de ces droites, et pflle positif 
de l'une d'elles celle des deux extrémités qui aboutit au plus 
aigu des deux sommets qu'elle rencontre. 
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Supposons qu'il se produise une défopinatioo générale trans- 
formant la splière décrite autour d'un point en un ellipsoïde. 
J'admettrai que la nouvelle distribution électrique sur la molé- 
cule peut être remplacée par un accroissement des masses - 
électriques des pûles, et que cet accroissement est mesuré par 
la distance infinitésimale de sa position primitive sur la sphère 
au point où son rayon vecteur perce l'ellipsoïde transformé de 
cette sphère. 

Soient alors les coordonnées et la masse d'un pôle représen- 
tées par les lettres ;»,„ et q, et 2 le symbole d'une sommation 
étendue aux n pôles, d celui d'un accroissement infinitésimal, 
j/i„ les coordonnées du centre de gravité des masses q, Q la 
masse totale. On a 

Qy, — '2qx, Qy^=Sgx,, Qs^—^gx^ 

avant la déformation, et 

(i/i + ày,)Q = Z{x, + dxt) (g + dq), ... 

après la déformation; or on a 

= 29 ^2(9 + dq); 
d'où 

Qàp, = I,{x,dq + g à Xi), .... 



Admettons que Iqàx, soitnégti^able relativement à ^x,âq, 
ce qui revient ù dire qu'à des déformations très petites toPres- 
pondent des dégagements finis d'électricité, ou encore qu'une 
couche électrique rigide sur la molécule conserve le même 
centre de gravité lorsqu'on déforme infiniment peu cette sur- 
face. 

Notre hypothèse fondamentale consiste à assimiler ày au 
déplacement électrique et à écrire 

, _ ^Xidg 
A- g ,- 

Pour un milieu homogène et un système de pressions le 
même en tout point, ce déplacement est constant. Par suite, 
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la densité électrique en volume, r'est-à-dire la convergence 
de f, est nulle partout : 

La densité superficielle aux surraces-limites du corps ou de 
iliscontinuité est, par délinition, la projection du déplacement 
sur la normale dn à l'élément da de ces surfaces. 

La quantité totale d'électricité dégagée sur cette surface sera 
donc, ^ et n étant les vecteurs du déplacement et de la normale, 



//« 



ou explicitement, 

'/(An, +A», + A».)ii', 



//'' 



et si la surface est plane et indéfmie, 
Q.f.ZQsnf. 

Il faut encore obtenir les valeurs dq des accroissements des 
masses électriques de chaque pôle. Considérons pour cela un 
deuxième système d'axes rectangles Ç,m lié au premier par la 
fonction [p]: 

? = [?]», 
et réqiiation 

%\ (1 + ;,)• + :|(1 + î,)' ^■ s; (i 4- s,)* - R', 
de l'ellipsoïde sur lequel sont venus se placer les points de la 
sphère |t + Sî 4- Sî ^ Ri, lorsqu'une déformation a réduit les 
coordonnées de $ à ; - - -t, les S étant des quantités de signes 
quelconques assez petites pour qu'on puisse négliger leurs 
carrés. 
Rapportée aux axes primitifs, cette équation devient 

l a-? [i+2(pî,î,-Hpî.3, +'V\M + ... 

( + 43;,a!,(p„p„î, +p„p„î, +Pj,P„3,) -+-...— U*. 
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On aura le rayon vecteur s de cet ellipsoïde suivant une direc- 
tion piézoélectrique quelconque, de cosinus directeurs am, en 
y faisant x,,, = fa,„, ce qui donne 

,J = * + 2(p,i ï, +Pi,ît, + Pijai)*S, + 2(p„ «, +Ph«, +p„a,)'5, 
+ 2(p„a,+p„a,+p„a,)'3v 

ou 

^=i + 2(0P.A)*î. + 2(SP,A)'S, + 2(0P,A)3,. 
f 

Prenons l'inverse de cette quantité, effectuons la division, 
puis Textraction de la racine, en négligeant S', et nous obte- 
nons, en soustrayant le résultat de R, l'accroissement du rayon 
vecteur, et par définition Taccroissement d'électricité du pôle 
considéré : 

^ = ^^ = {âP.A)*Ë, + (0P.A)'£, + (3P,A)'£3, 

où A et P|,j sont les vecteurs ayant pour composante a„i et 
les termes des diverses lignes de [p] d'après les notations du 
premier chapitre. 

Si l'on ordonne le développement de cette expression par 
rapport i'i a, on trouve pour coefficients de aî ,.. a,ai ... les 
dilatations et les glissements définis au §5 delà première partie. 
Car on a, en vertu des équations (13) du § 6, puisque S,,, 
représentent les dilatations principales suivant les axes ^, 

Formons ensuite les valeurs de /■,„ 



en posant 

et pour 



R=l, 0=1, 
Ix? — A,', lin} = A/ A,-, 2a;3jai — A^AjAi-. 
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Les valeurs définitives de /■„, sont 

( fi ^ AX +A*A,rfj-*-AÎA,</,+A,A,A,2fl',+ AÎA,2ff,+A'A,2tf,.. 
(3) j f, = A'A,rf,4-Aîd, +AÎA,rf,-i- AîAjîfl-. + AiA, Ajîj.+A; A, «fl-, , 
(/", ^AÎAX+AîAjrfjH-AJrf, + A;A,2?,+ AîA,2v,-i-A,AiA,2Sa- 

Le vecteur f est donc une fonction vectorielle sestlnôme, 
mais incomplète et non symétrique des dérormations. En dési- 
{(nant par /, ..., les six variantes d et g, je l'écrirai symboli- 
quement 

/■=|A|.l, 
et l'appellerai une fonction vectorielle fractionnaire. Nous ad- 
mettrons, d'autre part, que Ton a entre les déformations et les 
tensions l et X les relations 

x = \c/,i, d'où / — \c/r'x, 

et, par suite, 

(3„) /•=|A|;\c/7'X, soit f=\gU^X, 

c'est-à-dire que le déplacement électrique est une fonction 
sextinOme incomplète des tensions. Remarquons que l'inversion 
des fonctions vectorielles telles que (3) ne donne la valeur des 
six quantités i ou X en fonction des trois variables f que si les 
quantités sont, en outre, liées par trois équations. Dans le cas 
présent, les quantités peuvent s'exprimer par les équations (2) 
à l'aide des trois tensions ou dilatations principales. Ces der- 
nières sont donc simplement liées au déplacement électrique 
par une fonction linéaire à trois termes. 

Ce résultat aurait pu être posé a priori; il suffit, en effet, 
d'admettre que les composantes du déplacement électrique 
sont des fonctions linéaires des défonnatious et, par suite, des 
tensions. Mais les coefficients dépendent des angles que font 
entre eux les axes piézoéJectriques et les tensions, et doivent 
être déterminés par l'expérience. Notre méthode a l'avantage 
de déterminer les valeurs de ces coefficients. 



GooqIc 



^^ 
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§ 3. — Application de la méthode. 

Pour appliquer la méthode, je considère, en premier lieu, les 
cristaux dont la molécule a la symétrie de l'octaèdre ou du 
sphénoèdre, comme le chlorate de soude, le cuivre pyriteux, 
le sei de seignette; cas : (Oi), (c) et (d) du § i. 

Sojent 



les cosinus directeurs des rayons vecteurs joignant le centre 
pris pour origine aux quatre sommets d'un télraèdre isocèle. 
Je supposerai ces sommets chargés de masses égales à l'unité 
avant la déformation. On vérifie que les coefficients A sont tous 
nuls, sauf -ViAjAj, qui vaut «,0,0,, produit que je désignerai 
par a. Par suite, on a 

Les cristaux considérés n'ayant pas trois plans de symétrie 
rectangle, il est douteux qu'on puisse leur appliquer le système 
des coefficients (9), § 2 du chapitre précédent. Je l'adopterai 
néanmoiiis, puisqu'il s'agit surtout d'éclairer la méthode. On a 
donc Xi, .,., X, désignant les' tensions ; 

X,^c,.(., X.-c„(„ X, = c„(,. 

Introduisons les pressions principales 2,2, E, en l'onction 
desquelles ces tensions X peuvent s'exprimer à l'aide des for- 
mules (13) du § 0, cliapitre I, et nous obtenons 

/; — ^ (PiaPiiH, +p„p„2, +Ps,P„ï,), 
f> = ^iPnP„^,+P„P„~> + P,iPi,^z)' 
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Considépons im espace anisotrope, limiLé par une splière, 
soumis suivant un dianiètpc quelconque de direction n,nà une 
pression Si, les autres pressions SjSa étant nulles. Les équa- 
tions précédentes deviennent 

(4) A^f w.«,s„ /;-.... 

Le dégagement d'électricité est donc nul si ce diamètre est 
un des axes binaires du tétraèdre, et maximum si ses cosinus 
satisfont à 



Il se confond, dans ce cas, avec les hauteurs du tétraèdre si 
celui-ci est régulier, et il est facile d'en déduire le plan de 
nulle électricité. 

J'appelle surface piézoélectrique le lieu des extrémités du 
déplacement électrique /", lorsque l'extrémité du vecteur repré- 
sentant la pression unique S, décrit une sphère. Elle est le 
résultat de l'élimination du vecteur n entre les équations (i) 
et 7i't -i- n\ + n^ = i, et s'écrit 

K,» étant des constantes. 

Si ces constantes sont égales, comme dans le cas d'une 
anisotropie dépendant du système cubique, c'est une surface 
qui se déduit simplement de la surface de Steiner. 

Je ferai remarquer que si les quantités a. sont les cosinus 
directeurs des rayons vecteurs relatifs aux sommets d'un 
polyèdre régulier, l'application des formules précédentes con- 
duira à des surfaces ayant les mêmes symétries que ce polyèdre. 

J'appelle surface projetante piézoélectrique le lieu des extré- 
mités du déplacement électrique f projeté sur la direction de 
ta pression principale unique £, supposée constante et prenant 
toutes les directions possibles. 

Son rayon vecteur r mesure la quantité d'électricité dégagée 
sur les faces d'un cristal taillé perpendiculairement à ce rayon 
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vecteur, lorsqu'on le comprime dans la direction de ce rayon. 
On obtient cette surface en coordonnées polaires en éliminant f,^ 
entre les équations (4) et 

r — /■,!), + /"jWj + /"jWj ou Sp/b. 
En coordonnées cartésiennes, elle s'écrit 

in + n + nr = ur. x const. 

Un autre cas mieux étudié au point de vue expérimental est 
celui du quartz, c'est-à-dire des hémiédries des systèmes sénaire 
et rhomboédrique présentant trois axes électriques à 120" dans 
un plan (cas (ft) de la classification). 

J'admettrai que les trois pôles situés dans un même plan 
sont chargés, avant la déformation, de masses électriques 
égales ît l'unité par raison de symétrie. Les directions de leurs 
rayons vecteurs sont 



^/3 



2 2 







Les coetRcients de |A| seuls différents de zéro sont 

' * '4 

ce qui donne pour le déplacement 



.=-^., 


f. 


=!<-. 


+ «■). 


f, = o. 


On obtient des résultats concordant 


avec les 


i expériences en 


adoptant les coefficients du 


système sénaire : 




d. d. 


t. 


2j, 2», 


ig. 




*^ii f^ij 


^n 








X, 


c„ 


Cu 








X, 




'n 
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«.. 





X. 






c,. 





X, 
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X, — Xj — c,.(rf, ■ 



/; = -2cx., /; - - c{x, - Xj), A = o, 
u 

ti exprimant les X en fonction des pressions principales S. 
Supposons que la pression principale Si existe seule, les 
utres étant nulles, et posons 

p,, — cos<i)sin6, Pi, ^ sinu>sin6, p,, ^ cosfl, 

étant l'angle de cette pression avec Oa;,, et &> celui de sa pro- 
jction sur le plan Oxjar, avec Oa:,. On a alors 

fi = — cSsinïusin'O, U— — c2cos2û)Sin'6, A — 0> 
{\ 4- fl =:c'S'sin'fl. 

La surface piézométrique est donc {lorsque 6 varie) celle d'un 
ercle. Autrement dit, une compression constante en grandeur, 
ariable en direction et dans le plan des axes électriques 
) =90°), dégage une quantité constante d'électricité (propor- 
lonnelle à cette pression). La direction du déplacement élec- 
rique a est donnée par tga;= tg^u, et fait un angle droit ou 
ul avec celle de la pression, suivant que cette dernière est 
erpendiculaire ou se confond avec l'un des axes piézoélectri- 
ues. Si celle-ci, confondue d'abord avec l'axe Ox„ tourne 
utour de l'origine d'un angle u de Ox, vers Oxu le déplacement 
lectrique, d'abord dirigé suivant la direction négative de Oj:„ 
[)urne vers ta direction négative de Oj^i d'un angle double 2a. 
les résultats sont bien d'accord avec les expériences effectuées 
lar MM. Curie sur des lames de quartz taillées dans diverses 
lireclions, et par M. Rontgen sur une sphère de quartz. 
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On obtiendra la quantité d'électricité dégagée sur ia face 
d'une lame taillée perpendiculairement à la direction de cosinus 
Phi Pm Pm «t Comprimée suivant cette direction, en cherchant 
le rayon vecteur r de la surface projetante piézoélectrique. Sa 
valeur est 

v — f^p^^ + ^,pi; — — a3(3cos'w — sin'u)siD(i>sin'8. 

L'examen des divers systèmes cristallins a montré (§ 1) que 
les axes piézométriques étaient symétriquement rangés autour 
d'un axe principal. Je me bornerai à énoncer les résultats 
relatifs aux cas de trois et six axes. 

i" Trois axes de masses égales à l'unité, symétriques autour 
d'un quatrième de masse m. 



X' + iJi'=l. 



COORDONNÉES DES 


DÔLES. 










1 


m 





li 


x 


1 


\/3 


2 


X 


1 




9 


>. 


,, 



DÉPLACEMENT ÉLECTBIQCE. 



/» - ^ n'X(2d, + d,) + (3X' + f»)d,. 

Ces formules reproduisent celles déjà établies relativement au 
tétraèdre régulier, lorsqu'on adopte pour Ox, une des hauteurs 

2^3 I 

du tétraèdre. Il sufiit d'y poserX=:—5—' 1* = ;;, et de ren- 
verser la direction de l'axe central, ce qui revient à admettre 
m = ~~ I . 
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Elles s'appliquent à la tourmaline (6,, § !). Employons le 
système des coetlîcients d'élasticité du type sénaire,et supposons 
que les pressions X|X,Xi, dirigées suivant les axes de coordon- 
nées, existent seules; alors f, = 0, f, et f, sont des fonctions 
linéaires de ces pressions, et le calcul montre qu'une même 
espèce d'électricité pourra apparaître à une extrémité de Taxe 
central du cristal, pour des pressions perpendiculaires ou 
parallèles à cet axe. C'est un résultat constaté par MM. Curie. 

2° Six axes symétriques autour d'un axe central. 
Il suffit d'ajouter aux actions des pôles du cas précédent 
celles des trois autres ayant pour coordonnées 



±\, 



— n ±x 

Le signe — convient lorsque les signes de ces pôles alternent 
avec ceux des précédents. Le déplacement f qui en dépend a 
pour composantes 

3 

= 4^ 

En ajoutant ce résultat au précédent, il ne reste, en prenant 
le signe —, que l'action de l'axe piézoélectrique central. 

Comme il ne s'agit que d'une question de méthode, il est 
inutile de développer les calculs d'ailleurs faciles pour d'autres 
cas particuliers, d'autant plus que, sauf pour le quartz, des 
expériences de mesure n'ont pasencoreété faites. Les systèmes 
binaire et terbinaire présenteraient deux ou quatre axes symé- 
triques deux à deux autour d'un axe central. La fonction lAI 
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se simplifie beaucoup dans ces cas, mais, à cause des nom- 
breux coefficients d'élasticité, le résultat peut être compliqué. 

Je terminerai par un cas singulier qui, non rencontré jus- 
qu'ici, il est vrai, présente cet exemple de complication due aux 
coefficients d'élasticité, bien que son axe piézoélectrique soit 
unique. 

C'est l'hémiédrie OL, OC, OP du système asymétrique. Elle 
ne peut présenter qu'un axe piézoélectrique. En le prenant pour 
axe Oxjet en adoptant les deux coefficients d'élasticité de ce 
système, on aura 

Supposons qu'il n'existe qu'une seule pression E de direc- 
tion n,a, on aura 

La surface piézoélectrique se réduit évidemment à tous les 
points d'une droite. La quantité dégagée sur une lame taillée 
perpendiculairement à la direction n sera la projection de f, 
sur cette direction ou ftti,. 

On obtiendra le lieu des extrémités du vecteur représentant 
les pressions capables, en agissant aux extrémités d'un dia- 
mètre d'une sphère, de produire un même dégagement d'élec- 
tricité suivant ce diamètre, en égalant ^Hb à une constante, 
puis en remplaçant n,„ par -|r et S.' par ai + xj + x',. L'équa- 
tion résultante est du cinquième degré. 

g 4. — Phénomènes de réciprocité. 

Lorsqu'il n'y a que deux variables en jeu, le phénomène de 
réciprocité résulte de l'interprétation du théorème suivant : 
Soit u une fonction de deux variables x et y, et par suite 

1 vj iTj ji ■ f^X rfY 

ilu = Xdx + Ydy, dou -3-=- — 

dy dx 

Donc, l'accroissement que subit la quantité X lorsque y croît 
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de dij, la variable x restant constante, est égal à l'accroissement 
subi par Y lorsque x croit de dx, la variable y restant constante. 

Considérons, comme exemple, une lame diélectrique déga- 
geant la quantité d'électricité Q lorsque la pression exercée sur 
ses faces croît de zéro à P. On a 
(8) Q = a?. 

D'autre part, soit L Paccroissement de longueur de la lame 
lorsque la différence de potentiel devient $ sous pression cons- 
tante. On doit avoir 
(8.) L ^ - a* O- 

En effet, l'énergie électrique Q* est créée tout entière aux 
dépens de l'énergie mécanique LP, l'énergie totale est nulle, 
et l'on a 
(6) Q* = -LP. 

Cette équation, combinée avec (5), fournit (5,). La différen- 
tielle exacte du est ici 

/QL\ 



d(aP*) ou -^ir^) 



suivant que l'on considère comme variables indépendantes P, * 
ou 0, L. 

Pour étendre ce théorème au cas d'un nombre quelconque 
de variables, c'est-à-dire à un système de pressions quelconques 
et à un champ électrique quelconque, je rapporterai la quantité 
d'électricité et la force de pression à l'unité d'aire, ce qui revient 
à leur substituer le déplacement électrique Z' et la pression X. 
La relation entre ces vecteurs, l'un trinôme, l'autre sextinôme, 
est, d'après (â„), une fonction sextinôme d'indice f qui rempla- 
cera la première équation (5), et que j'écrirai 

r,„ = | «M •■■«,. |x,--, 

(6) }„... ou /■=1«1.X. 

(<) Lîppmann, C. A. de VAcadÉmic des Sciences, L XCII, 1881. 
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De même, en pappoptant l'allongement L et le potentiel * à 
l'unité de longueur, c'est-à-dire en leur substituant les déforma- 
tions et les dérivées tp,,, de "Ii par rapport aux coordonnées, la 
deuxième équation (5) devient une fonction fractionnaire dont 
je désignerai l'indice par |, et que j'écrirai 

?,..., I S^^ S,i S,, I î,j, 

(7) ou /^|sN?. 



Elle exprime que les six déformations d et 23 sont des fonctions 
linéaires des trois dérivées de 4>, dont nous allons déterminer 
les coefficients. 

Les énei^ies électriques et mécaniques doublées et rapportées 
â l'unité de volume sont : 

/,X.-h ?,X, ou Six. 

Substituons, dans ces scalaires, Si f et l leurs valeurs (fi) 
et (7); leur somme doit être identiquement nulle, et par suite, 
les coefficients des produits X©, ce qui donne 



-«,1, *u = -?îi. 



C'est la relation de réciprocité généralisée. 

En éliminant X et l entre (5) et (7) et la relation 

X - VU, 
on a 

(8) f=\g\t\cU]s]if, soit |6|ç. 

Le développement de cette fonction trini3me montre qu'elle 
est symétrique en même temps que la fonction \cj. A ces rela- 
tions il faut encore adjoindre celle dont on a déterminé les 
coefficients à l'aide des éléments de symétrie du milieu : 

(9) r-|A||I. 

& et 5 sont les coefficients piézométriques à pression ou à 
potentiel constant; s et c~', les dilatations dans ces mêmes 
conditions. 



:flnvGoOgle 



9k B. ËLIE. 

On se rendra peut-être mieux compte de ces diverses rela- 
tions en remarquant que les variables i et X se transforment 
respectivement, lorsqu'on change d'axes de coordonnées, delà 
même façon que se transforment les coefficients d'un ellipsoïde. 
Appelons ellipsoïdes des dilatations et des tensions ceux qui 
ont respectivement pour coefficients i ou X et ont pour équa- 
tions 

La relation \t=\cj,1 exprime que l'un des ellipsoïdes est 
complètement déterminé par l'autre. 

La relation f^ \6/ç exprime que les vecteurs f et o sont 
aussi déterminés l'un par l'autre comme étant des directions 
conjuguées de l'ellipsoïde de \bj. 

Ënlin, par les relations 

/•=|,||X el i = -\«/.,, 

on voit qu'il existe la même correspondance entre l'ellipsoïde 
de tension et le déplacement qu'entre l'ellipsoïde des dilatations 
et la force électromotrice ç. 

Une application très simple de ces considérations peut se faire 
à l'aide du premier exemple traité au paragraphe précédent, 
pour lequel on a 

'' ~ * * ~ c„ * ~ c^^ ^'' 

g 5. — Dilatations calorifiques et phénomènes pyroélectriques. 

Les phénomènes électriques dus à la chaleur peuvent se 
classer en deux catégories que l'on désigne sous les noms de 
pyroélectriques et thermoélectriques, et que je définirai de la 
façon suivante. Soit T la température à l'époque i d'un point 
de coordonnées x,„. L'expérience apprend que si en ce point 

dT 
la dérivée -j- n'est pas nulle, il y a dégagement d'électricité. 

Ce phénomène a reçu le nom de thermoélectrique. Je le laisserai 
complètement do côté, car son étude exige d'autres considéra- 
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tioM que celles abordées dans ce travail. D'autre part, si eii 

(JT 
un point la dérivée -j- n'est pas nulle, autrement dit s'il y a 

échauffement en ce point, l'expérience montre qu'il y a déga- 
gement d'électricité. Ce phénomène rentre dans la catégorie de 
ceux appelés pyroélectriques, et va faire l'objet des remarques 
qui suivent. 
Limitons-nous au cas où réchauffement, c'est-à-dire où la 

rfT 
dérivée -^ est constante en tout point. La question à résoudre 

est de trouver la valeur du déplacement électrique. 

La seule différence entre ce problème et celui relatif à la 
piézoélectricité consiste en ce que l'ellipsoïde de déformation 
est remplacé, dans ce cas, par l'ellipsoïde des dilatations. Le 
premier dépend du système de pressions que l'on applique et . 
des coefficients d'élasticité du milieu. L'étude de la pyroélec- 
tricité dépend donc de celle de l'élasticité, à laquelle je l'ai 
rattachée. Le second ellipsoïde dépend des dilatations calorifi- 
ques, sur lesquelles il est nécessaire de dire quelques mots. Cet 
ellipsoïde est le lieu des points en lesquels se transforment les 
points d'une sphère lorsque la température varie. Dans le cas 
des milieux à trois plans de symétrie rectangle à toute tempé- 
rature, l'orientation de cette surface ne change pas, et les 
vecteurs d'un même point à deux températures sont liés par 
une fonction vectorielle symétrique, dont les coefficients dépen- 
dant de la température sont assujettis à la condition de main- 
tenir aux axes de l'ellipsoïde les mêmes directions. 

Il n'en est plus de même dans les milieux à axes cristallo- 
graphiques obliques; l'ellipsoïde des dilatations peut changer 
de forme et d'orientation lorsqu'on passe d'une température à 
une autre. Tout au plus y a-t^il lieu de supposer que les axes 
obliques du cristal, en se modifiant en grandeur et direction, 
forment à toute température un système de trois diamètres 
conjugués de l'ellipsoïde des dilatations. 

Mais la connaissance complète de ce dernier n'est pas néces- 
saire en pyroélectricité. L'examen des diverses hémiédries fait 
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dans le pfemiep paragraphe montre, en effet, que les axes piézo- 
électriques sont rangi^s symétriquetnent autour d'une direction, 
et que cette direetlon, d'après sa position dans le cristal, doit 
être considérée comme un axe principal de l'ellipsoïde des dila- 
tations. Par suite, la somme des composantes perpendiculaires 
à cette direction des déplacements électriques dus aux dilata- 
tions suivant ces axes piézoélectriques est nulle; les compo- 
santes suivant cette direction sont seules à considérer. Il en 
résulte qu'il n'existe qu'un axe unique de pyroélectrieité, 
comme Texpépience l'indique d'ailleurs. 
Le déplacement électrique est donné par la formule 

d étant la dilatation de l'axe de pyroélectricité du cristal. La 
densité électrique sur les faces d'une lame taillée perpendicu- 
lairement à une direction quelconque sera la projection de f 
sur cette direction. 

Les conséquences des théories précédentes sont d'accord 
avec les expériences faites jusqu'ici, qui se sont bornées ii la 
mesure du dégagement d'électricité dans quelques directions 
et seulement à observer le sens du phénomène dans tes autres. 
Les fonnules proposées, en assignant la quantité d'électricité 
dégagée suivant une direction quelconque, faciliteront le travail 
expérimental. Il est vrai qu'elles supposent connu le système 
de coeflicients c qu'il faut adopter dans la fonction liant les 
tensions aux déformations. 

Ces coefficients pourraient môme différer do ceux admis dans 
la théorie ordinaire de l'élasticité, puisque les déformations dont 
il a été parlé sont plutôt celles de la molécule elle-même que 
celles du milieu. Pour cette raison, peut-être ont^ils une dépen- 
dance plus grande que ces derniers de la symétrie cristallo- 
graphique. C'est du reste ce que nous a montré l'examen 
des divers cas particuliers que nous avons traités. 
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CHAPITRE III 



Snrfaces ëqulpolenllelles et lensioas électroslallqnes 
dans les milieux anisotropes. 



§ 1. — Généralisation des vecteurs solénoïdaitx et lamellaires. 

On sait combien en physique est fréquent le cas où les coni- 
pusantes d'une grandeur sont les dérivées partielles relative- 
ment à x,a d'une fonction potentielle, Ainsi, dans la science 
(le l'électricité, le déplacement électrique, le déplacement 
magnétique, l'intensité d'un courant jouissent de cette propriété, 
c'est-à-dire que l'on a, en désignant par ç,n les trois déi'ivées 
d'une fonction *, 



Sir \V. Thomson a appelé lamellaires de tels vecteurs. Il a 
appelé solénoïdaux les vecteurs dont les constituants vérifient 
la relation 

du, du. du, 

dXi ax, ax. 

Je désignerai le premier membre de cette équation par con- 
vergence de u et l'écrirai 

On peut voir facilement que t out vecteur est la sgjQjne^ d'un 
vecteur solénoïdal et d'un vecteur l amellaire. En effet, repré- 
sentons par A, le paramètre de deuxième ordre d'une fonction 

rf* d' d' 
ax\ ax\ aw\ 

et posons, « étant un vecteur quelconque, 
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ns que m' soit solénoïdal et u' lamellaire, c'est-à-dire 
a fonction potentielle *; nous en déduisons 

Cm =3 A,*, 

ni détermine* si Cu est donné. L'étude d'un vec- 
inque se ramène donc à celle des vecteurs lamellaires 
laux. Telle est la manière d'envisager la question 
ilieux isotropes. 

la modifier dans le cas des milieux anisotropcs. Or, 
; la plus simple à faire sur la relation entre les déri- 
)tentiel et les composantes du vecteur est de les 
ées par une fonction vectorielle linéaire, c'est-à-dire 
m désignant le vecteur par f, 

te hypothèse que nous adopterons pour l'étude de 
lénomènes qui se produisent dans les milieux aniso- 

, dans la première partie, que cett^ fonction peut 

. \b\+_\b\ _ . m~\b\ 

f— --i^ ç + , _t_ ç, 

dit qu'elle est la somme de deux fonctions, l'une 
s, l'autre gauche. Le vecteur relatif à la fonction 
isfait à la condition 

c soIéflDïd»!. Ss direction êst parpendifiulaira «us 
orées, c'e8t>à*dire aux normales aux surfaces équi-* 
s, 

u vecteur sym^rique, il ne peut être assimilé coin- 
à un vecteur lamellaire. Sa direction n'est pas 
Lix surfacsB ëqulpotentlelles, mais seulement conju' 
3cteur çiM relativement à l'ellipsoïde de la fonction 
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vectorielle \6/, Il en est de môme, par coniéquent, de» lîgpes 
da flux qui ROTit définies par l^s équations 

dXi dx, âxj 
Pourtant la détermination de la fonction ♦, loraqu'an con- 
naîtra les valeurs de f en certains espaces ou sur certaines 
surfaces, se fera comme dans la théorie ordinaire du potentiel, 
et nous admettrons cette généralisation. Ainsi, désignons, 
oomma^pn le foit en physique, par densité de volume du vec^ 
teup f, la convergence changée de signa, fit popons 
-, = Gr. 
Soit e la densité Buperflcielle, c'est-Mire la quantité ; 
^e=x fin, -h z;», + fiHt, on tfn, 
n étant la normale en un point d'une surface. 

La densité Cf, lorsque f= \ft/<p, devient, si l'on écrit çu, 
f,„ ... pour les dérivées secondes de *, 

ou 

Â,( étant le paramètre du second ordre généralisé. 

La fonction potentielle en tout point x^ s'exprime par la 
formule 



-Hm-im- 



Les intégrations s'étendent à tous les pointe de coordon- 
nées a,M en lesquels les densités i ou « ne sent pas nulles. La 
fonction rt_i est définie par la relation 

D»-irVi== 67î{a'i— «i)' + — +'i^ï{^%—<'i) (a?.— 0))+ —i 
ou, symboliquement, 

(» - a) \fe/-' (aï — 0), 
dans laquelle D^-i est le déterminant de la fonction \&/*'- 
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jette généralisation des théorèmes relatifs à la théorie du 
ientiel était nécessaire pour nous permettre de traiter les 
estions qui vont nous occuper dans les paragraphes suivants. 

§ 2. — Lignes de force dans les milieux anisotropes, 
[Jne application très naturelle de la fonction vectorielle 

présente lorsqu'on étudie soit la propagation de la chaleur 
de rélectricité dans les conducteurs, soit la distrihution des 
dacements électriques ou magnétiques dans les suhstances 
lectriques ou magnétiques, lorsque le milieu qui est le siège 
phénomène est anisotrope. Le vecteur f représente alors 
t le flux de chaleur ou d'électricité, soit le déplacement 
îtrique ou magnétique, et o,,, sont les dérivées relativemeni 
; coordonnées soit de la température, soit du potentiel élec- 
[ue ou magnétique. Les coefficients b caractérisent l'aniso- 
aie du milieu, et celle-ci peut être naturelle ou produite 
ficiellement par des forces extérieures telles, par exemple, 
! des actions mécaniques, magnétiques, etc. 
,e cas d'une anïsotropie provoquée artificiellement est par- 
dièrement intéressant, parce qu'il établit une relation entrer 
forces de natures différentes et qu'il offre l'exemple de 
jrèmes de réciprocité remarquables. Il a permis, en outre, 
expérimentateurs, de mettre en évidence comment les lignes 
flux étaient- distribuées relativement aux surfaces équipo- 
:ielles. C'est ce côté particulier de la question que je com- 
ice par étudier. Que le milieu soit naturellement isotrope 
non, les coefficients b peuvent être modifiés en appliquant 
milieu un système de pressions mécaniques X,,,,...,,!. Cette 
iification a été établie par une expérience de Sir W. Thom- 
{*) consistant à tordre un tube conducteur dont les géné- 
ices sont parcourues par des courants, et à observer la 

Voir Wiedetnitnn's Amiakn. 1878. 
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déviation d'un aimant placé dans Tintérienr du tube. L'hypo- 
thèse la plus simple est de supposer que les coefficients crois- 
sent des quantités a, fonctions linéaires des six tensions et, 
par suite, des six déformations Z,...» qu'elles produisent. II y a 
même lieu de penser que les quantités a forment un détermi- 
nant symétrique, et par suite de poser, en désignant par c des 
constantes, 

La symétrie de la fonction [6| ne serait pas altérée, et l'expé- 
rience seule peut répondre à ce sujet. 
Si dans la fonction 

f=\b-i-a\<i ou a — \cl,t, 
on fait les substitutions 

f=[p]r, 9 = [pw, t==\[p]\i-, 

on obtiendra, en l'identifiant à 

f = \f +a'\>f', 

les valeurs de 6' + a' en fonction de fc + a et de l'. On aura 
ainsi passé d'un système de coordonnées à un autre. 

Supposons le milieu primitivement isotrope et les axes dirigés 
suivant les trois dilatations principales que produit le système 
de tensions, ce qui revient à poser 

b,i=p + g {di + rf, + d,) -H »d^„ et &(, = 0; 

les substitutions indiquées donnent, pour les valeurs des b 
dans le nouveau système d'axes, 

b'H = p -*■ q{\ + ^,-h 8,) + <îi„ et b'n = «Ym» 

oii p est la conductibilité primitive du milieu q, et s d'autres 
constantes de conductibilité introduites par la déformation du 
milieu. 

Les modifications de conductibilité dues aux actions méca- 
niques ont été peu étudiées. Au contraire, les modiflcattons 
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lues au iiidgnéli&mg ont faitlé sujet des travaux expérimentaux 
le MM. Hall, Leduc, Rlghî, etc., et elles sont d'autant plus 
Dtéressantes qu'elles nous donnent l'exemple d'une fonction 
reutopielle dyssymétrîque et, par suite, de la distribution sîn- 
{ulière des lignes de flux relativement aux surfaces équlpoten- 
ielles, distribution que je vais indiquer dans quelques cas 
jarticulieps. 

Supposons le milieu primitivement isotrope, indéfini et uni- 
'ûrmément magnétisé. Choisissons de* axes de coordonnées 
'ectangulaires tels que la fonction \èl soit ramenée à sa forme 
îanonique et puisse s'écrire explicitement 

*» — P. Pi, 



es coefTicients 6 et p étant des fonctions des composantes du 
;hamp magnétique. 
La densité, dont l'expression changée de signe est 

ioit être nulle partout, sauf aux électrodes par lesquelles 
entrent ou sortent les courants. S'il n'existe qu'une électrode 
\ l'origine, on satisfait â cette équation par 



(67'icf + bj'xl + b-i^xD^ 

Les surfaces équipotentielles sont donc d^ ellipsoïdes; elles 
sont indépendantes des coefRcients p qui caractérisent la dyssy- 
métrie, et ceci aura lieu quel que soit le nombre des électrodes. 

Les lignes de flux sont déflnles par les équations 

^3}^ àm^ 

^i~P»bî'ar^ + p,&7'3', ~" i»!— p,6j'ic, + pjfrr'a^s 
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Si les coefficients p sont nuls, elles sont des droites issues de 
l'origine et, par suite, ne rencontrant pas l'ellipsoïde à angle 
droit. Dans le cas contraire, elles sont données par les inter- 
sections des systèmes de surfaces ayant pour équations 

K (^ + y + 1^) = (Pi-T, -I- p,», + p,!r,)', 

6, 6) Pi ' 

K et c étant les constantes d'intégration et 9 un angle défini par 

(l'a —(Cl y ' >Vi^A l'ange 

Admettons que le magnétisme ne modifie pas la conducti- 
bilité principale, c'est-à-dire que &, = b, ^ &s = 6 ; ces équa- 
tions, en prenant pour axe Ox,, la direction du vecteur p,, ou 
en posant p, := p, := 0, se réduisent à 

K(x* + «;*♦• »*) = bpxî, 

d se, 

«1* + xl = Ce^- «t taDgOat'f- 

Les lignes de flux sont donc des spirales logarithmiques sur un 
cône de révolution, ce qui met en évidence le caractère rota- 
toire du coefficient p. 

Les expériences ont toujours été faites dans des espaces à 
deux dimensions, c'est-à-dire avec das couches minces de ma- 
tière! conductrices étendues sur un plan. On a dans ce cai 

f—kti + P?« 



La densité égale à 



V'éi 



+ *, 



à'\ 



3ÎÎ * »' fflj 
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doit être nulle partout sauf aux électrodes. Posons 



La fonction <b est égale h une somme de termes de la forme 
A, log r„. 

Dans le cas d'une source unique à l'origine, les lignes équi- 
potentielles sont des ellipses, et les lignes de flux les spirales 
logarithmiques définies par l'équation 

Vx*, -h s,'=Ke-9 — !_!. 
P 
Dans le cas de deux sources égales et de signes contraires, 
on sait et on s'assure, d'après les équations précédentes, que 
dans les milieux isotropes les lignes équipotentielles et les 
lignes de flux forment un double système de circonférences 
orthogonales. Si le milieu devient anisotrope, les lignes équipo- 
tentielles définies par t 'équation 

r,—rtX const 

sont des ellipses indépendantes de p et résultant d'une double 
extension suivant deux axes rectangulaires d'un des sys^'-mes 
de circonférences. 
Les lignes de flux ont pour équation -^ 

9,-9. .. rr- 
-^ = e ' X const. 

r, 

Ce sont des courbes joignant les deux électrodes et présentant 
une inflexion. 

Ces exemples nous paraissent suffisants pour mettre en évi- 
dence la relation qui existe entre les vecteurs a et f an point 
de vue de leur distribution dans l'espace. 

§ 3. — Définition des déplacements électriques. 

La fonction vectorielle f= \b\<f trouve une application, 
avons-nous dit, dans l'étude de la distribution des déplacements 
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ûleclriques ou magnétiques d'un milieu électrisé ou magné- 
tisé. Cette application va nous conduire à étudier cette fonction 
sous un autre point de vue. Des expériences dues à MM. Duter, 
Righi, Quincke ont montré, en effet, que les milieux diélectri- 
ques éprouvaient des modifications sous l'influence de l'élec- 
trisation; que, d'autre part, la capacité d'un condensateur était 
modifiée sous l'influence de pressions mécaniques; que, par 
suite, les coefficients diélectriques, qui sont précisément les 
coefficients b, devaient être fonctions de ces pressions, et par 
cela même des déformations qui en sont la conséquence. 

Déjà Maxwell, sans connaître les déformations d'un diélec- 
trique électrisé, avait montré que les actions électrostatiques 
pouvaient être représentées par un système particulier de ten- 
sions. MM. Helmholtz, KirchofF, Lorberg ('), etc., ont repris la 
question pour les milieux isotropes, en tenant compte des 
déformations. C'est le même sujet que nous allons traiter ici, 
mais en considérant des milieux anisotropes. 

Je rappelle d'abord la définition du déplacement électrique 
qui a été donnée au chapitre précédent, relatif à la piézoélec- 
tricité. Nous admettons qu'à l'état neutre le centre de gravité 
d'une molécule matérielle et de la masse d'éther qui lui est 
liée coïncident. Lorsqu'un corps non conducteur est électrisé, 
ces centres se séparent, et leur distance est ce que j'appelle 
déplacement électrique. Il correspond à ce qu'on désigne sous 
le nom d'électricité libre. Je le représente par f. Ce déplacement 
résulte de deux autres : 

i° Sous l'action des forces électriques extérieures et en sup- 
posant que la molécule pondérable ne puisse ni se déformer 
ni se déplacer relativement aux molécules voisines, l'éther 
subit un déplacement que je représenterai par f. 

2" D'autre part, il faut admettre que l'éther exerce des actions 
mécaniques sur la matière pondérable, que par suite la molé- 

(') Helmholti, Wiedemann's Annalen,iSS^ ; KîrcholT', Lorherg, même recueil. 
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cule se déforme; c'est cette déformation qui nous a servi à 
expliquer les phénomènes piézoélectriques. 11 en résulte un 
changement dans la position relative des centres de gravité de 
l'éther et de !a molécule. De plus, la molécule elle-même, 
abstraction ftiite de sa déformation, se déplace relativement ii 
ses Voisines, ce qui produit un nouveau changement dans la 
situation des centres de gravité. Nous désignerons, sauf à les 
distinguer plu» tard, l'ensemble de ces deux derniers déplace- 
ments par la lettre f. 
D'après notre hypothèse, nous aurons 

f=r-r- 

Enfin, si 9,^ est la force éleotromotrice dérivée d'une fonc- 
tion potentielle, nous admettrons que les déplacements sont 
liés aux forces électromotnces par les fonctions vectorielles 
suivantes : 

Il en résulte entre les constantes aa'a' la relation 
a^a' + 4ita'. 

Ces constantes, au nombre de vingt-sept, se réduisent à 
deux, a' et ft', dans les milieux isotropes, parce que a' = 1 
avec les unités adoptées. 

Les densités en volume relatives aux trois déplacements, 
qui sont par définition leurs convergences changées de signe, 
satisfont aux relations 

t — t'"- 1'. 
De même on a pour les densités superficielles, 

e == ^ &fn, c' =s — Sf'n, »'=:"" Bf'n, 

e=- e' — e', 

n étant la normale dirigée vers l'extérieur de la surface consi- 
dérée. 
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Les tensions électrostatiques qua nous cherchons se dédui- 
pont de la variation de l'énergie électrique du milieu. I^ est 
donc nécessaire de former l'exprasaion de cette énergie, ainsi 
que celle de la fonction potentielle. 

Pour la fonction potentielle, je considérerai comme connus 
les théorèmes qui lui sont relatlfo; sa valeur en un point quel- 
conque est 

en écrivant r pour r,, quantité qui a été définie plus haut, et 
en étendant l'intégrale à tous les points de l'espace. On peut 
encore écrire, en intégrant par parties, 

ou, en désignant par m les masses électriques tda et «du, 
(1.) * - 2 " 

Nous supposons nulle la densité de volume du diélectrique 
il l'intérieur duquel nous cherchons les tensions. En tous les 
autres points distribués en des volumes ou sur des surfaces, 
les masses électriques ont des densités données. Les a ont des 
valeurs assignées partout. On sait que, dans ces conditions, la 
fonction potentielle peut être déterminée partout^ ainsi que les 
déplacements. 

L'énergie doit être mise sous une forme proposée par 
M. Helmholta et commode pour le calcul. Par définition, 
l'énergie Ë, par unité de volume, est le demi'scalaire des vec- 
teurs ^et <f, c'est^ù-dire 



(2) ■= = 2///«ft^- 

onsidéré, ce qu'on pt 

. OU de / oeul, 



pour le volume total considéré, ce qu'on peut écrire, en l'expri- 
mant à l'aide de f seul ou de / seul, 
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lie intégration par parties permet d'écrire encore 

E = 5 II / e*(ici + 5 (/c*dii, ou encore â2]'~~' 

titroduisant l'une ou l'autre des valeurs * écrites plus haut, 
orme adoptée dans ce qui va suivre s'obtient en retran- 
it l'expression (2*) du double de l'expression (2J. On a 



— Tension» éîectrostatiquee dans un diélectrique anisotrope. 

i\xt obtenir la valeur des forces pondéromotrices, il faut 
itenant former la variation de cette expression relativement 
coordonnées. On sait, en effet, que, dans un état d'équi- 
t, le travail élémentaire, autrement dit le scalaire formé 
'- les trois composantes de la force et les déplacements 
lels SiCira, est nul. Les coefficients de oXi„, lorsqu'on aura 
vé la variation de l'énergie, seront les forces cherchées. 
expression (3) présente cet avantage que, dans la variation, 
leut supposer constante la valeur de la fonction potentielle ; 
I, en effet, pour la variation relative à *, 



-mi- 



]drs. 



r, si l'on intègre le dernier terme par parties, on trouve 
ectivement pour coefficients de S*, dans les intégrales 
tives au volume ou à la surface, 

(£+ Cf)dts et (C4- Sfn)d<i, 

itités qui sont nulles par définition dans l'état d'équilibre. 
î variation devra donc porter uniquement sur les coeffi- 
ts a. Comme l'énergie dont il est question ici est l'énergie 
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libre, nous savons que a est de la forme a' -{-aux. Les 
coefficients a' résultent de l'addition de deux coefficients : l'un, 
que nous continuerons d'appeler a, et qui dépend de la défoi-- 
ination de la molécule, n'aura d'existence que si cette défor- 
mation entraîne un changement dans la position relative des 
centres de gravité de l'éther et de la molécule; l'autre, que 
nous appellerons oT, dépend des déplacements relatifs des mo- 
lécules entre elles. L'hypothèse la plus simple est de supposer 
que l'ensemble de ces coetTicients a" est une fonction symé- 
trique sextinôme des aix déformations </,„, g,a. 

Pour former la variation de E, je ferai d'abord abstraction 
des coefRcjeats x" et ne considérerai que la somme a'+ à^a.', 
que je continuerai de désigner par a. Ces a sont des fonctions 
de Xin- On a, pour la variation du dernier terme de l'énergie (0), 
relative à Sic, sous forme explicite, 

+ l3-^Çi+'T-^?!+-r-^l?iSa-, dcj 



\àxi ^' 



ou symboliquement 



///' 



dxi\ 



Le coefficient de oxi, dans cette expression, peut s'écrire 
autrement à l'aide du lemme suivant, dont la vérification peut 
("tre facilitée en ramenant la fonction |a| à sa forme canonique, 
et cette forme est symétrique, puisque la fonction gauche, qui 
résulte de la décomposition de [a] en somme, n'entre pas dans 
la formation du scalaire qui représente l'énergie. Ce lemme 
consiste en ce que, si on pose 

(y»- Aï, y,,3=/;ç, y„=i(-/-.ï.-A?,+/;î,) 
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et 

da;, oï, àxj 
on a identiquement 

(7) -«f.-f.-K'i&.l'' 

La substitution de pe terme et de deux autres résultant des 
variations de Sa;,&r„ dans le dernier terme de l'expression (3), 
permettra d'écrire symboliquement 



(8) 



^ ! rrr? |a|çdn = ^ fffs&('ibx)dcs—fff6{nx)dr3. 



Cherchons ensuite quelles forces introduisent dans le diélec- 
trique les variations des densités e et c produites par la varia- 
tion Sx, en supposant le potentiel constant, dans l'expression (3). 

On a 

es = ~ SX, + -r — 0», + -r— ô 

àx. ax. ox^ 



- Sa;,, 



quantité qui, introduite dans le premier terme de (4), donne, 
après une iqlégratiop par parties, 

^^^ sfjJe*d'3=jJ*eS(Sa;.w)^5-j7T|8(ç5ai)dR. 

Les variations (8) et (9), substituée» dans l'expression (4), la 
réduisent à 

Or, la quantité ï{ni!3îi + »,8fl;,4- njSxj), qui entre dans le 
deuxième terme, n'est autre que la variation te de la densité 
è la surface, n étant la normale dirigée vers l'extérieur du 
milieu. Ce terme est égal et de signe contraire au dernier 
terme dans lequel le sens du déplacement de est dirigé à l'in- 
térieur. Il ne reste donc que 

(10) fff^^' ^'^^ "^ ^»'^* "'' ^''^^'^ '^''' 
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quantité qui, intégrée partieltement, donne l'expression 



im 



2J rrCY.in. + Y„n, + Y„B,)8a!i«îî, t = 1, 3, 3, 



dont les parenthèses représentent les forces agissant sur une 
surface de séparation. 

Le milieu se trouve donc soumis, pap unité de volume, aux 
forces Fi,3 (6), $t par unité de surface, aux presaîoni Y (5), 
absolument comme un corps élastique ioumis à des actions 
mécaniques. La seule différence est que le système de pres- 
sions (H) est dyssymétrique et que le milieu serait soumis, 
par unité de volume, aux couples Y,, — Y,„ .... Ces couples ne 
disparaissent que dans le seul cas où, conjme dans les milieux 
isotropes, les coefficients a se réduisent à un seul. Par consé- 
quent, si aucun phénomène physique correspondant ne peut 
équilibrer ces couples, il faut admettre que tout corps, quelle 
que soit sa constitution moléculflire, est isotrope au point de 
vue électrique, et qu'il ne devient tel que par les déformations 
qu'il subit sous l'iniluânce de l'électricité, , 

Dans ce cas, le système des tensions (5) se réduit à celui 
proposé par Maxwell duns le g 10S de squ traité d'éUctricité. 

Considérons maintenant la variation des coefficients que 
nous avons désignés par a" et que noua avous admis être liés 
par une fonction symétrique sextinôme aux six déformations l. 
Si ^ sont Us eeeffioients de cette fonotien, on a 

(11) «'■5::;\ft/,J. 

La variation lîoit porter, dans ce cas, «ur les déformations 
h ... !.. 

La valeur explicite du facteur de tîn, dans ta variation de 
l'intégrale 

est àomée par la lamme des termes du tableau suivant, dans 
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lequel on a inscrit les multiplicateurs des coeSlcients en tête 
des lignes et des colonnes : 



(«)( 



[2?.?) P.i 
Intégrons par parties et posons 



{«) 



„ =Zi,- =&;,?*+ ?,-,îî + ^,?.' + i^U<i..f, + 2&*«f,1 



(14) 



^ dZ» dZ;, dZ,, 

aXi àXj oXf 



= i, 2, 3. 



La variation cherchée de l'intégrale devient 
— ^ riT(G,3«, + G.Stt, + GjStiOdo 

Ces quantités doivent être respectivement retranchées des 
quantités (10) et {10^). Les fauteurs de Stt;, sous les intégrales 
triples et doubles, sont les forces pondéromotrices agissant sur 
l'unité de volume ou l'unité de surface. Les forces Z ne don- 
nent naissance à aucun couple et peuvent être équilibrées, 
comme en élasticité, par des forces mécaniques extérieures. 

Lorsque le milieu est isotrope, les quantités a'4-4sa' se 
réduisent à une seule, les constantes 3 à deux, et l'on retrouve 
alors les formules indiquées par MM. Kirchoffet Lorberg. 

Les équations (H) et (13) mettent en évidence le phénomène 
de réciprocité. Il consiste en ceci : lorsque les forces électro- 
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motrices varient et que les actions mécaniques extérieures sont 
maintenues constantes, comme en vertu de (13), les pressions 
électrostatiques internes sont fonctions des forces électromo- 
trices; le milieu doit se déformer. Au contraire, si, laissant 
constantes les forces électromotrices, on déforme le milieu, 
comme les coefficients diélectriques a", d'après les équations 
{H), sont fonctions des déformations, les déplacements élec- 
triques, et par suite la quantité d'électiicité libre, sont modifiés. 

L'analyse que nous venons de faire montre que les coeffi- 
cients diélectriques a peuvent se classer en trois catégories. 
Les coefficients de la première, qu'on a désignés par a.', ne 
dépendent que de l'anisotropie de l'éther intermoléculaire; ils 
se réduiraient à un seul si, en l'absence de toute déformation, 
le milieu qui transmet les actions électriques était isotrope. 
Les a' ne dépendent que de la déformation de la molécule 
matérielle; ils ne peuvent différer de zéro que si elle est dys- 
symétrique, comme dans les milieux piézoélectriques. Enfin, 
les a"' ne dépendent que des déformations mécaniques, c'esl^- 
dire des déplacements relatifs des molécules; ils sont nuls si 
ou empêche ces déformations à l'aide de pressions convenable- 
ment appliquées. Les a.' et a' pourraient alors recevoir le nom 
de coefficients diélectriques à volume constant. 

Dans cette étude des coefficients b d'une fonction vectorielle, 
j'ai toujours parlé de conductibilité ou de tensions électrosta- 
tiques, afin de fixer les idées et pour m'adresser à des expé- 
riences réalisées. Il est clair qu'on pourrait donner à / et o 
d'autres significations physiques et supposer les coefficients b 
fonctions de quantités autres que le magnétisme ou la défor- 
mation. 
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CHAPITRE IV 
>nde5 électromagnétiques dans les milieux anlsotropes. 



§ 1. — Énergies potentielle et cinétique. 

A l'aide des considérations développées dans la première 
lartie, nous allons maintenant étudier les mouvements vibra- 
oires auxquels on attribue certains phénomènes soit électri- 
[ues, soit lumineux. Cette étude, qui nous conduira à une 
lUpface des ondes très générale et comprenant comme cas 
larticulier celle de Frcsnel, est une simple affaire de calcul, 
me fois admises les formules qui servent de point de déparl. 
lais l'établissement de ces dernières n'est pas aussi facile. 
)oit-on les regarder comme la traduction la plus simple de 
âits expérimentaux isolés? Peut-on les déduire d'hypothèses 
dus générales? Quelle est la signification mécanique précise 
les grandeurs qu'elles renferment, ou du moins rencontrons- 
lous leurs analogues dans d'autres théories mieux connues? 
l'est sur ces questions que je me propose d'insister tout 
l'abord, en généralisant les notions acquises dans les problèmes 
irécédemment traités. 

Admettons, en premier lieu, que le phénomène à étudier 
ippartiennô à la catégorie de ceux qui dépendent d'un vecteur. 
Lvec les cofflpcsantSB de ce vecteur fondamental, leurs déri- 
'ées par rapport au temps et aux coordonnées et les coordon- 
lées elles-mêmes, on peut former une infinité d'invariants. 
ious appellerons énergie élémentaire d'un phénomène dans un 
nilieu isotrope le produit d'un invariant par une constante. Si 
e milieu est anisotrope, nous avons vu que son énergie se 
léduisait de celle relative à un milieu isotrope en substituant 
lU vecteur fondamental et à la coordonnée des fonctions vcc- 
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torielles inverses Tune de l'autre, et dont les coefficients 
doivent être regardés comme de simples nombres d'homogé- 
néité néro; 

On distingue deux espèces d'énergies, suivant que l'invariant 
dépend ou ne dépend pas des vitesses du vecteur fondamental. 
Dans le premier cas, elle est dite cinétique; dans le deuxième, 
potentielle. 

Soient T et P ces énergies relatives à un volume quelconque-, 
et formons la variation de la quantité 

relativement au vecteur fondamental u. On appelle force le 
vecteur dont les composantes sont les coefficients de Siii». 
Dans la mécanique des corps indéformables, le vecteur n'est 
autre que la coordonnée qui définit la position du corps. Dans 
le cas des milieux déformables, le vecteur fondamental u est 
le déplacement d'un point du milieu de sa position d'équilibre. 
Les lois de l'élasticité résultent de l'application du principe de 
moindre action que l'on vient d'énoncer aux fonctions 

P = r (Ai; + 4BI,)(ïcs, T^Tp^'iido. 

Les divers états des corps sont définis par les valeurs des 
coefficients de compressibilité et de rigidité A et B et du coeffi- 
cient d'inertie (densité) p. 

Lorsque la somme P + T est invariable, le milieu est dit 
parfaitement élastique. Si cette condition n'est pas remplie, il 
y a une dissipation de l'énergie qui se traduit par un dégage- 
ment de chaleur. Le milieu est dit alors visqueux. L'expérience 
a montré que la quantité d'énergie pendant le temps t pouvait 
être représentée par 



IX' 



(A'Ii' + iK\;}ilda. 
Bans cette expression, 11 1[ sont les invariants li 1, où l'on a 
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remplacé le vecteur u par sa vitesse u; A' et B' sont les coeffi- 
cients de viscosité-. Le quotient A'/ A de deux coefficients 
homologues est un temps appelé module de relaxation du 
milieu (en électricité, la résistance spécifique). 

Les coefficients de Sui^dt dans la variation de cette expres- 
sion sont des forces appelées forces de frottement que l'on 
devra ajouter aux forces trouvées précédemment. On aura ainsi 
déduit du principe de moindre action l'équilibre et le mouve- 
ment des milieux à la fois élastiques et visqueux. 

Ces propositions relatives aux milieux matériels étant ad- 
mises, il suffit, pour expliquer les faits électriques, d'admettre 
que l'espace puisse être le siège d'un deuxième vecteur fonda- 
mental V, auquel tout ce qui a été dit relativement au vecteuru 
est applicable. 

Ce deuxième vecteur caractérise'un deuxième milieu qui peut 
coexister avec le milieu matériel. Les lois de son équilibre et de 
son mouvement peuvent être diff'érentes, parce que les inva- 
riants avec lesquels on forme les énergies peuvent ne pas être 
les mêmes que ceux relatifs à l'élasticité des milieux matériels. 

Ainsi, la théorie des actions électrostatiques de Maxwell se 
réduit à adopter pour invariants le tenseur et la convergence 
du vecteur fondamental supposé lamellaire, tandis que l'énergie 
des phénomènes électromagnétiques est formée avec le tenseur 
de la version du vecteur supposé solénoïdal. A ces invariants, 
il faudra encore joindre le tenseur de la vitesse si le milieu est 
en mouvement, et enfin, si le milieu est conducteur, une fonc- 
tion dissipatrice formée avec les versions de la vitesse du vec- 
teur. C'est cette dernière qui conduit aux lois d'Ohm. 

Ces formes d'énergies sont encore insuffisantes pour expli- 
quer les actions mutuelles de la matière et du milieu électrique 
telles que les attractions des corps électrisés ou magnétisés, 
les phénomènes d'électrostriction, ceux de piézoélectricité.etc. 
Dans ces cas, il faut faire intervenir, comme on l'a vu dans les 
deux chapitres précédents, des invariants dépendant à la fois 
des vecteurs fondamentaux u et y. 



:flnvGoOgle 



PO.XCTION VECTORIELLE. 117 

Toutes ces énergies, qu'elles dépendent d'un seul vecteur ou 
de plusieurs vecteurs fondamentaux, sont de même nature, 
c'est-à-dire qu'elles peuvent s'associer par addition. Il n'en est 
pas de même des vecteurs et des forces que l'on en déduit. 
Par la variation d'une énergie relativement à u, on obtient les 
forces mécaniques; par la variation relative au vecteur électri- 
que V, on obtient les forces appelées électromotrices, et non 
associables par addition avec les premières. 

Il y aura lieu, dans ce qui suit, d'appliquer la théorie précé- 
dente d'une façon partielle ; mais grâce à cette vue d'ensemble, 
l'établissement des équations différentielles du mouvement 
électrique présentera une marche plus naturelle. 

§ 2. — Équations difféi-entieUes du mouveTnent loreque l'inertie 
et l'élasticité du milieu oii( les mimea plana de symétrie. 

On sait que les vibrations qui peuvent se propager dans un 
milieu élastique sont de deux sortes : les unes, longitudinales, 
dépendent d'un vecteur lamellaire ; les autres, transversales, 
ont un vecteur solénoïdal. C'est de ces dernières dont nous 
allons nous occuper. Le milieu dans lequel elles se propagent 
n'est pas la matière pondérable; nous le supposerons parfaite- 
ment élastique, et par suite nous n'aurons à considérer ni les 
énergies de dissipation, qui n'existent que pour les milieux 
conducteurs, ni celles dépendant du tenseur et de la conver- 
gence du vecteur, lesquelles conduiraient à des vibrations 
longitudinales d'ordre électrostatique. L'unique énergie poten- 
tielle que nous adoptons est celle formée avec l'invariant R, 
dans lequel on a fait les substitutions 



C'est donc, en écrivant C pour Ca,a,a„ 
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Nouç adopterons pour l'énergie cinétique l'expression 



(2) 



'-imm-^'m^-my 



Lee constantes 3p aont les coefficients d'inertie du milieu. Ce 
sont les analogues de la densité des milieux matériels isotropes. 
Mais comme nous considérons ici un milieu rendu anisotrope 
par son association à la matière, il est clair que ces coefficients 
d'ineptie doivent varier avec la direction. 
Fermons U variation 



5 r(P + T)dT, 



et après avoir intégré par partie les termes relatifs à l'énergie 
potentielle, annulons les coefficients de 8u sous l'intégrale 
triple. Eu égard à la notation 

(3) !H,(,)^p-.p. 

le résultat peut s'écrire 

Sous l'intégrale double, le coefficient de 5v, est 

X„»i + X„«, ■+■ X„«„ 

n étant le vecteur normal à un élément superficiel et les X 
ayant pour valeur 

X„ = 0, Xi, = - Xji = C<L^%{i>). 

Ces quantités représentent les tensions existant dans le milieu 
et permettent d'écrire le second nombre de l'équation {^) sous 
la forme 

<ÎX„ (»X„ dX„ 

QXi ça?j àx^ 

Ce qui vient d'être dit est relatif au cas où l'élasticité et 



:flnvGoogle 



FOIfCTtON VRCTOBIËLie. 110 

l'inertie du milieu, représentées respectivement par les cons- 
tantes a. et ^, ont des axes de symétrie communi; le cas plus 
général sera traité ultérieurement. Je remarquerai aussi que 
l'on aurait pu adopter l'invariant I, au lieu de l'invariant R, 
comme on le fait en élasticité. Les résultats eussent été les 
mêmes sous l'intégrale de volume; mais le» tensions sous l'in- 
tégrale de surface auraient été différentes et auraient vérifié 
Xy ^ X,j . Le choix de ces invariants doit donc devoir se faire 
par l'étude des phénomènes de réflexion et de réfraction. Ce 
cîlioix ne joue aucun rôle dans ce qui suit, puisqu'on ne consi- 
dérera qu'un milieu continu indéfini. 

Comparons maintenant ces équations à celles que l'on ren- 
contre en électromagnétisme. Ces dernières renferment quatre 
vecteurs: h et f, les déplacements magnétique et électrique; 
h et e, les forces magnétomotrice et électromotrice. 

Ils vérifient les équations vectorielles 

(5) 6 = V/ft. 

(6) f^\Tile. 

Dans les milieux isotropes, les constantes ji et K se réduisent 
à deux. Ce sont les coefficients magnétique et électrique qui 
entrent dans les lois de Coulomb relatives au magnétisme et 
à l'électricité ; 

(fi) force =AK^. force = -r-^^' 

La constante A introduite par M. Helmholtz peut être sup- 
posée égale à l'unité, en modifiant les homogénéités de ft et K. 

Les vecteurs doivent, en outre, vérifier les deux groupes 
d'équations diiférentielles 

(7) -Al.. ^ = ^.(6), 
(.8) AK,^ = rB,(ft). 

Les énergies magnétique et électrique sont formées respecti- 
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ent avec les scalaires de b.h et de f.e. Nous les écrirons, 
upprimant le facteur 1/4 «, 






Hiftj' + [jij,A;)dt3, 



+ K,e; 4- K,eî)(io. 

^Ue est la forme que les physiciens donnent actuellement 
problème de la propagation des ondes électromagnétiques. 
3ur trouver le sens mécanique des grandeurs électriques, 
tifions l'énergie magnétique (9) à l'énergie potentielle (1), 
'énergie électrique (10) à l'énergie cinétique (2). On en 
:lut 



=v/iLm .=Ya 






valeurs doivent, en outre, transformer respectivement les 
jpes (7) et (8) en les groupes (3) et (4). Ceci exige que 
ait 

& = ,, |; = 4, . = 1,2.3, 

Ji étant des constantes et, de plus, 

;ant une vitesse. 

ans les milieux isotropes, cette dernière condition devient 

AKK^ = i- 

systèmes d'unités électrostatiques ou électromagnétiques 

sistent à supposer 

soit A — 1, ]A=!, soit A = 1, K = 1. 
îmblerait peut-être plus naturel de supposer 
C = 1», K = p, d'où A = - > 
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OU encore 

C = -, K = -, d'où A=p. 
V. p 

Dans la première hypothèse, la version et la vitesse du vec- 
teur électrique représentent les forces ; dans la deuxième, les 
déplacements. On serait arrivé aux mêmes résultats en identi- 
fiant les énergies électriques et potentielles et les énergies 
magnétiques et cinétiques. 

Nous conclurons de cette discussion que les principes de 
mécanique adoptés dans l'étude de l'élasticité des milieux 
matériels conduisent aux équations que les physiciens, en 
s'appuyant sur des considérations toutes différentes, ont pro- 
posées pour expliquer les mouvements électromagnétiques. 
Mais l'étude de ces mouvements est insuffisante, et il faudra 
s'adresser à d'autres piiénomènes électriques pour préciser 
l'homogénéité des vecteurs et des constantes qui entrent en jeu. 

La symétrie des équations (7) et (8) conduit à quelques 
remarques intéressantes. Ainsi, l'élimination entre elles soit 
de e, soit de h, donne 

équation que l'on peut écrire 

(13') A'K,K,i.j -^ = i«(A.) ^, 



en généralisant les symboles A et C. 

Ces équations sont toutes semblables à l'équation (i). De 
même que celle-ci a été déduite des fonctions P et T, on 
pourrait déduire les équations précédentes de quatre fonctions 
formées avec les quantités e, h, K et [*, comme le sont les 
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fonctions PetTavecîf, a et p. En les désignant parT,TiP,Pi, 
on vérifie, à cause des groupes (7) et (8), que 

(15) T, = Pj, Tj = P,. 

On peut aussi arriver au résultat obtenu par M. Pointing dans 
le cas de milieux isotropes, que raccroissement par une unité 
de temps de l'énergie totale d'un espace s'exprime à l'aide d'unR 
intégrale étendue à la surface. Multiplions, en effet, les grou- 
pes (7) et (8) respectivement par les facteurs ^jft; !((«(, et ajou- 
tons les résultats. Si on intègre tes seconds membres par 
parties, on obtient 

c'est l'expression du théorème énoncé. 

§ 3. — Cas où l'élasticité et l'inertie du milieu ont des plans 
de syméti-ie différents. 

Arrivons maintenant à l'eKamen de l'hypothèse plus générale 
consistant à admettre que les énergies cinétique et potentielle 
ne puissent être ramenées simultanément à leurs formes cano- 
niques (10) et (9) par un choix d'axes convenables. Nous 
pourrons exprimer ce fait en disant que les axes de symétrie 
du milieu relatifs à l'inertiç diffèrent des axes de symétrie 
relatifs à l'élasticité, ou encore, pour parler le langage des 
électriciens, qu'il peut coexister dans le même espace deux 
milieux de symétries différentes, l'un électrique, l'autre magné- 
tique. 

Dans ce cas, nous déduirons encore les équations difFérçn- 
tielles du mouvement de la variation de l'action. Mais les 
énergies cinétique et potentielle dépendent des fonctions vec- 
torielles o--'- '-*'-- J- 
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symétriques, et en nous servant d'un ou de deux points pour 
représenter les dérivéoa première ou seconde par rapport au 
temps, nous les écrirons : 

(17) 2P = 33î(v) \a/-'^(D), 2T = 0é \&/ c. 
Il en résulte, pour l'une des équations d'équilibre ; 

( ^^^""^'^"^ -*■ "''^'^"^ ■*■ *^'-'^'^*'^ 

ou, symboliquement, 

(18) -Wv = ^\«/-'g!(o). 

On introduira les quantités ft, e, K et p. à l'aide des relations 
(H) et (12) du paragraphe précédent, en y supposant G et p 
égaux à l'unité et en mo'diflant l'homogénéité de a et p. 

L'équation précédente deviendra 

{18') -\K/e=^^W-'S(*). 

Posons alors 

(19) -h = \iLl-'^{e) on _\|^/A = S(e); 
elle se transforme en 

(20) \K/ë=^(ft). 

Enfin, l'élimination de e entre (19) et (20) donne 

(21) -\tt/A = S\K/-'3ï(A). 

Les égalités (19), (-20), (18') et (21) remplacent les égalités 
(7), (8), (14) et (13) du paragraphe précédent, mais sont 
d'une forme plus compliquée. C'est de ces égalités que 
M. Heaviside a déduit la surface des ondes électromagnétiques 
dans le cas général. 

Je vais montrer que, par l'emploi des coordonnées obliques, 
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il est possible de ramener les équations précédentes aux équa- 
tions canoniques du paragraphe précédent, et par suite d'arriver 
il la surface des ondes par un procédé de calcul qui reste le 
même, que les milieux électrique et magnétique coïncident 
ou non. 

Rappelons un résultat obtenu au chapitre I. Soient ;r et ^ les 
coordonnées obliques d'un même point rapporté à deux sys- 
tèmes d'axes obliques; 

« et V les projections obliques d'un même vecteur sur ces 
axes; 

X, I, « et tT les contraordonnées et contravecteurs corres- 
pondants, c'est-à-dire les projections orthogonales de ces vec- 
teurs sur les axes ; 

IpI et |d| une fonction oblique et son inverse. 

Entre ces quantités existent les relations 

(22) x=.\p\l, u=\p\v, 

(23) û = W|ë, ^(i) = |p|S(r). 

Pour transformer les équations (19) et (20), représentons 
par e et 1] ce que deviennent les vecteurs eeih dans le nouveau 
système d'axes obliques Ç, l'ancien système x étant rectangle, 
et opérons avec les substitutions (23). Le groupe (19) devient 

(24) -|pl^r= WWIÎ on -^!=lnlW\n|V. 

Le second membre de l'égalité (24) est une fonction vecto- 
rielle symétrique. Or, on a vu que la fonction \a\ peut être 
déterminée de façon à réduire ce deuxième membre à sa forme 
canonique. Il suffit d'annuler trois de ses coefficients, ce qui 
revient à écrire, d'après les notations du chapitre I, 

(25) SitA:*/ii,=0. t=i = i,2, 3. 
Posant, de plus, 

(26) S,t;V/^i = ^, 

nous ramenons le groupe (20) à la forme (7). On passera de 
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même des égalités (19) à la forme canonique (H) en posant 

(27) St.,\K.I-!Cj = Q, Si:,\K/T:i = Ki. 

Je remettrai à un prochain paragraphe la démonstration de 
ce théorème, que les axes obliques auxquels on rapporte ainsi 
les vecteurs ne sont autres que les diamètres conjugués com- 
muns aux ellipsoïdes inverses des fonctions \K/-' et \ji/~'. 

§ 4. — Surfaces des vitesses. 

Si nous considérons les positions qu'occupent, après l'unité 
de temps, les ondes planes parties simultanément d'un même 
point que nous prenons pour origine, on appelle surface des 
vitesses le lieu des pieds des perpendiculaires abaissées de 
l'origine sur ces ondes. Son rayon vecteur est donc l'inverse 
de ia vitesse Y de propagation. 

Partant des équations (7) et (8) et pour ne pas multiplier 
tes notations, représentant par eeth l'amplitude des vibrations 
correspondant aux déplacements électrique et magnétique, je 
substitue aux vecteurs e et ft tes quantités 



^"b r ) 

(t n,x, -4- n^x, -f- n,xA, 
T l )' 



At„cos2 

ce qui transforme ces équations en 

(28) Pj», - e,n, = [i,A,V, k,n, - li,it, = - K,e,V, 

d'où l'on déduit 

/ 1*1*1^1 + ^aKe^ + [ijAiSj ou Sbe = (i, 

\ K,e,ftt + K,eA + K^c,/*, ou Sfh = 0, 

I i/,fi,n, + ;j.jAjn, + [ji,A,nj ou $bn = 0, 

{ K,e,«, + K,p,n, + KjC,», ou Qfn =0. 

Si l'on suppose que les angles sont rectangles, les équations 
(7) et (8) correspondent au cas oij les plans de symétrie des 



^29) 
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milieux électrique et magnétique coïncident, n est ia normale 
au plan d'onde; V = ;^ est la vitesse de propagation de l'onde 
dont la distancée l'origine est L à l'époque T. Les équations (28) 
ou (29) expriment que le déplacement magnétique (induction 
de Maxwell) & ou rik est perpendiculaire à la force électrique e 
et à ta normale à l'onde, tandis que le déplacement électrique/" 
ou Ke est perpendiculaire à la force magnétique ft et à la nor- 
male à l'onde. 

Dans la figure ci-jointe, n représente la normale à Tonde 
perpendiculaire à 6 et à /"; m est la direction du rayon lumi- 
neux perpendiculaire à h. Les arcs indiquent les angles droits. 




Dans le cas où les axes sont obliques, les équations (7) 
et (8) supposent que les deux milieux sont rapportés à leurs 
l>ians de symétrie communs. Les quantités e, ft et n sont alors 
des contra vecteurs que je m'abstiendrai de représenter pare, h 
et n pour ne pas encombrer récriture. Au contraire, x„ Xs, xs 
sont les coordonnées obliques d'un point du plan d'onde dont 
l'équation est 

(29„) Sarii = a:,n, + x'n, + XjW, — \ = ~> 

L étant sa distance à l'origine ù l'époque T. 

On obtient, en effet, cette relation, si l'on projette sur la 
perpendiculaire abaissée de l'origine sur l'onde le contour 
formé par les coordonnées obliques x,a d'un point de l'onde. 
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La condition pour que deux vecteurs x et n soient pei'iiendi- 
ciilaires est alors : $xn = 0. 

Ces définitions posées, il n'y a plus lieu de faire une distinc- 
tion entre les cas de coordonnées obliques ou rectangles si. 
Ton veut déterminer la surface des vitesses. 

Éliminons e ou h entre les groupes (28), et nous les trans- 
formons en les groupes suivants, dont je n'écris qu'une équa- 
tion sur trois : 



(30) 



[ (K,n; + K.nî - ia,K,K,V')A, — K,n,«,A, - K,n,n,ft, : 



L'élimination de A ou de c de chacun de ces groupes conduit 
au même résultat, qui est : 

( ;x,ii.t.,K.K.K,V* - j:.,K,(ii,K, + i..K,)nî 

[ + (v-inî + !J-i»i' + l*inî}(K,«î+ K,«î + KjnJ) — 0. 

On a la surface des vitesses en posant nm = î/i«V. Si le 
système d'axes est oblique, n et y doivent être considérés, 
d'après ce qui a été dit précédemment, comme des contravec- 
teurs n et y. 

Lorsque les fonctions vectorielles symétriques \K/ et \ii/ sont 
ainsi ramenées à leurs formes canoniques, il est facile de 
rendre explicites les calculs symboliques à l'aide desquels 
M. Ileaviside, conservant aux fonctions leurs formes générales, 
est arrivé à cette surface. D'après la notation scalaire, on a 

1*1 1** t*i 
K,s! + K,}i + M; = 0»\6/!,, |i + |î+|i = Sj\K/-'i,, 

quantités que je désignerai encore, pour abréger, par M, M"' K, 
et K~*. Posons, de plus, 

li,lj.,;j.. = ;*, K,K,K. = /t, 



idByGoOgle 



lia B. EUE. 

et appelons et f deux vecteurs proportionnels aux déplace- 
ments magnétique et électrique 6 et /. 
Les égalités (30) peuvent s'écrire ainsi : 

V^ _ _ Jfi_ fi _ 

Chacun de ces groupes se transforme et, à l'aide des équa- 
tions (29), devient 

On peut encore multipllei' par tj— les différents numérateurs 
des égalités (32), puis ajouter les résultats. La somme des 
premiers rapports donne l'unité, et par suite ou a 

(M), (33) et (31) sont cinq formes équivalentes de la surface 
des vitesses. 

§ 5. — Surface des ondes. 

La surface des ondes est, par définition, l'enveloppe des 
ondes planes émanées de l'origine après l'unité de temps. Elle 
résulte donc de l'élimination du vecteur 9 entre les secondes 
équations (32) et leurs dilTérentielles, avec la condition que le 
vecteur y satisfasse à l'équation de la surface des vitesses (29). 
On serait conduit aux mêmes résultats en opérant sur le pre- 
mier groupe (32). 
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Or, les équations (33) écrites ainsi*: 
(35) ' », = „{^-^), y...... 

puis multipliées par fm et ajoutées, doonent 

,36) „„M(|i + |l-.|?)-.(i; + ^ + g, 

en tenant compte de la condition que le vecteur 9, propor- 
tionnel au déplacement Ke, est perpendiculaire au vecteur y 
d'après (29). On pourrait écrire ce dernier résultat en notation 
scalaire 

ây!H.ls-fîîlKl-'? - i*0f|i*|-*9 = 0; 

mais je préfère conserver les formes explicites pour plus de 
clarté. 

Prenons les différentielles de (35) relatives à y et 9, et mul- 
tiplions les égalités obtenues ainsi respectivement par ^,„. La 
somme des résultats, retranchée de la différentielle de (36), 
donne l'équation suivante, débarrassée des différentielles de ç : 



(37) 



l'4"'^''{^yi^r)\ = ' 



D'autre part, si Xi„ sont les coordonnées obliques d'un point 
de la surface d'onde correspondant au point de la surface des 
vitesses dont les contraordonnées sont ;/,„, on a, par la pro- 
priété des enveloppes, pour des valeurs quelconques de dy, 

(38) »iàSi 4- aî.^y, -t- «,(ïy, oa Sxày = 0. 

La comparaison de (37) et (38) doone 

Multiplions les deux termes de ces rapports paf fi:) et 
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ur valeur commune est alors 



kJ V K. K, K. y 



aî.| 



pès les valeurs (35) de ç, le deuxième facteur du 
n'est autre que le premier membre de (33.), et par 
. On a donc 

ou Bxe — d et 0a;^=:O ou 0aîfi=O. 
V- 

est donc perpendiculaire aux vecteurs e et h. 
riété va nous conduire à la surface des ondes. 
ir W la distance de l'origine au point Xm de cette 
m la direction de cette droite, ce qui donne 

3 qui précèdent deviennent 

0mA = O, âffle = 0, 

0m« = =P^=:cosVW, 

ecteur ot n un contravecteur. 

lions ou par de simples considérations géométri- 

ait les groupes suivants : 

£«> — -_^ 
ijk sont différents et égaux à 123. On peut les 

itions, toutes semblables à celles (^), on déduira 
nséquences, à la condition de substituer à V l'in- 
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verse de W, à i* et K leurs inverses, et à la coiitraordonnée y ht 
coordonnée x. 

L'équation explicite de la surface tirée du développement 
d'une des formes (34), (33) ou (31), est 

ij (1/1 1 \ , 

Cette surface, dans le cas d'axes rectangles, a été étudiée (') 
comme étant celle des singularités du complexe du second 
ordre. Ses sections par les plans de symétrie sont des ellipses, 
cap, pour x, = 0, on a 

Elle présente donc les quatre points coniques et les quatre 
plans tangents singuliers réels de la surface d'onde de Fresnel. 
On peut appliquer à son investigation le procédé de Lamé (*). 
En effet, cette surface est l'enveloppe du plan 

ttiXi + «,ir, + B,3;» = V, 

les paramètres n et V étant assujettis à satisfaire à l'équation 
ries vitesses (31). Or, si on pose 



elle est encore l'enveloppe du ^lan 

lia 
les y étant assujettis à satisfaire à l'équation i 
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d'ondç. Dii plua, le pUn considéré est le pUa poUire au 

point X relativement à l'ellipsoïde 

(41) 2;!'.''^sîXif;=<- 

ii3 

Les deux points x et y, appartenant à la surface d'onde, sont 
donc conjugués, en ce sens que chacun d'eux est le pâle, rela- 
tivement à l'ellipsoïde (41), du plan tangent à la surface au 
point correspondant. 

Je crois intéressant d'indiquer à quoi tient la simplidcation 
qui se présente pour l'équation de cette surface lorsqu'on la 
rapporte à des axes obliques. M. Heaviside, en la rapportant à 
des axes rectangles, écrit son équation 

1 ^0.-p\Ktt/-'iP -+- 0a!\K/-'*. !5aj\|i/-'ï = 0. 
où l'on a explicitement 

(42) Sa!\ii/~**=:it-'3;;+ ... + 2(*7,'rc,a!, + ..., 

(43) 0a!\K/-'ic=;E7i'aî;+ ... + 2K7;ic,aTs+ ..,, 
.(Di.D.Sj;\Kn/~'a:=(K„i*„-K„ii,.+K„;.,.-K„iJ.„)iFÎ+ .. , 
1 +2(K„n„-K„a„ + K„m,-K„i.,.Kï.+ .... 

D, et D(i sont Içs déterminants des fonctions \K/ et \ji/ ; liiî' ■ ■ . 
les inverses (iAnfij. — ni)- fr ■■■ des coefficients de ces fonc- 
tions; les termes supprimés se déduisit de ceux écrits par 
permutation des indices, 

La réduction à laquelle on est arrivé provient de ce que les 
trois ellipsoïdes correspondant à ces trois expressions admet- 
tent un môme sy&t^me de diamètres conjugués. C'est ce que 
je vais établir. 

Transformons, à cet effet, les expressions (42) et (43) à l'aide 
de la fonction vectoriellQ oblique 

»=ipi«, 

et annulons t«s doubles préduita des variaMes %. Les égalités 
qyi en résultent sont, pour la fonction (45), 

BPi\K;-' p; = t),i diJérenu et êgaiu a 1, J, 3). 
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Le résultat rendu explicite pour deux de ces équations est 

/0=p„(K-'p„+K7i'p„H-Kr,'pM)+î'«(Kr>i»+K7.Pu+K-î/)„) 
,J +P..(K-'p.,+ K7,V„+K-V„), 

)o=Pi.(Kû'Pi»+K7>„+Kïi'p„)+p„(Kr,'p„+Kr,'p„+Kr,'p») 

' -*-pH(Kri'Pi.+K7,'p„+K7,'/>„)- 

L'élimination de pn entre ces dernières donne, en exprimant 
les p et les K~' par leurs inverses js et K, la dernière des éga- 
lités suivantes, dont les autres -s'obtiennent par la même 
méthode : 



(48) 



Pu . ~ Ptl 

_ K„an + K«o„ + ^»|Oti 
\ ." . Pm 

On aurait de même 



(47) 



Pu ' 

d 'où, en vertu de la relation xsnPn + o«pii + tiuPn = 0, 

(48) a7t.\K/iv = 0; 
de même, 

(49) 55t<\WiV-0. 

Ceci démontre tout d'abord la proposition énoncée sans 
démonstration à la fin du § 8. 

En second lieu, si on égale à s tes quotients des rapports (46) 
et (47) et si on élimine o des résultats, on obtient une équa- 
tion cubique en s qui déterminera les grandeurs, puis les 
directions des diamètres conjugués communs. 

En troisième lieu, si on désigne par ai,a„... les coefficients 
du troisième ellipsoïde (44), il ftiut établir qu'ils vérifient les 
égalités 

(50) 5itr\o/ iv = ou SPi \a/-* P; = 0, 

où i diffère dej. 
Or, si on op^re-sur (48), (49) comme oiî a opéré sur (45), en 
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ssociant une équation du groupe (-48) à une du groupe (49), 
in retrouve des égalités du type (46), (47), dans lesquelles les 
oefficients K sont remplacés par ceux a du troisième ellipsoïde. 
Les équations (50) se trouvent démontrées. 

g 6. — Polarisation rotatoire. 

Je ne me propose pas de traiter les nombreux problèmes 
[ui se rattachent à la propagation des ondes. Pourtant, l'un 
['eux, celui relatif à la polarisation rotatoire, offre un exemple 
i simple de l'emploi que l'on peut faire des notations em- 
)loyées précédemment, que j'en dirai quelques mots. 

Je considérerai d'abord l'invariant signalé à la fin du § 7, 
chapitre I, 

S(X(,U .•.j = I,2,3, 
>îi 

_, dei dvj^ _ à%(v) â^{v) 

'^~ dxj"^ âxt' ^~ dxj "^ dx, ' 
,e coefficient de 8v, dans l'expression 



3 r Z(\r)da, 



lonnera les forces qui produisent la rotation du plan de pola- 
■isation soit dans les milieux isotropes, soit dans les milieux 
laturellement anisotropes ou rendus tels par des actions méca- 
niques. Dans le cas d'anisotropie, on doit substituer à u et a; 
les fonctions vectorielles symétriques qui, par un choix d'axes 
convenables, peuvent s'écrire ; 



des axes sont ou les axes principaux du cristal, ou ceux de 
l'ellipsoïde des pressions mécaniques que l'on exerce. La force 
jue l'on en déduit dépend des dérivées troisièmes en x et 
n'écrit explicitement 

' dxl ' dxi ' dxî 



:flnvGoogle 



FONCTION VECroniELLB. iéH 

On sait que l'existence de ces termes suffit pour expliquer la 
rotation de la vibration dans les cas énoncés plus haut. 

Mais la polarisation rotatoire peut se présenter dans d'autres 
circonstances. Elle peut être provoquée par des actions élec- 
triques ou magnétiques. Dans ce cas, l'invariant servant ù 
former l'énergie devra être une fonction de la direction de 
l'action magnétique. Désignons par pi„ les cosinus de cette 
direction, par Cm trois constantes, et soit posé 

à à â â 

dp "^'dx, "^'diP, dx^ 

L'énergie élémentaire pour un milieu quelconque soumis à 
une action magnétique sera 

dT ' dp "*" àT ' dp '*' dT ' dp ' 

La force que l'on en déduit a pour composantes 

à d%{i>) 
dT dp 

On retrouve la forme proposée par Maxwell lorsque le milieu 
devient isotrope et que les constantes c sont égales. 

Il sera, dans tous les cas, facile de passer de la variable v 
aux variables h et e, et de trouver ce que deviennent les deux 
groupes d'équations (7) et (8) existant entre ces variables. 

§ 7. — Comparaison avec la théorie de Freanel. 

La théorie des ondes électromagnétiques, telle qu'on vient de 
l'exposer, est-elle la même que celle des ondes lumineuses de 
Fresnel? La réponse à cette question s'obtiendra en examinant 
le cas particulier oii l'un des milieux fictifs considérés est 
isotrope. 

Supposons d'abord le milieu magnétique isotrope, c'est-à- 
dire 1*1 = 1*, ^ |*j. Les conséquences que l'on en tire au point 
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de vue de l'orientation des divere vecteurs peuvent se traduire 
par la âgure oi-cootre. 




Cette figure correspond à la théorie électromagnéticiuede 
Maxwell si on considère la vibration lumineuse comme repré- 
sentée par la force électrique e. Cette foKe est perpendiculaire 
au rayon lumineux m et oblique au plan de l'onde. Elle corres- 
pond à la théorie que Lamé a fondée sur l'élasticité si la force 
magnétique h représente la vibration lumineuse. Cette vibration 
est dans le plan d'onde, mais non dans celui du rayon m et de 
la normale à l'onde n. Ni l'une ni l'autre de ces hypothèses ne 
satisfait à la théorie de Fresnel, qui exige que la vibration soit 
dans le plan d'onde et dans le plan mOn. Si on suppose le 
milieu électrique isotrope, c'est-à-dire K, = K, = K„ il suffit 
de permuter les lettres e, h, fi. en h, e, K dans la %ure, et le 
résultat énoncé n'est pas essentiellement modifié. ' 

11 en est tout autrement si on considère la vibration lumi- 
neuse représentée non par les forces électrique et magnétique, 
mais par ce que Maxwell a appelé les déplacements électrique 
ou magnétique f ou b. Supposons, en effet, ^, = jt|=> ot, et 
remplaçons e par sa valeur = dans les équations (14), 

Elles deviennent 

ât' "" dic, Ui àse, K, dxj dx^ \K,(*a;, K. dxj' 
Les trois équations ainsi obtenues, différentiées et ajoutées, 
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àl' \<)is, 



dx. ox,/ 



et sont compatibles avec les conditions admises dans ta théorie 
de Fresnel. Mais alors deux modifications sont nécessaires dans 
l'expression de l'énergie. En premier lieu, l'énergie cinétique 
doit être le carré de la vitesse de la vibration et non une fonc- 
tion quadratique de ses composantes. En second lieu, l'énergie 
potentielle n'est plus une fonction quadratique des versions. 
On se souvient (chapitre 11) que, dans le passage d'un milieu 
isotrope à un milieu anisotropo, les versions n'avaient été 
conservées qu'en transformant v etx par des fonctions vecto- 
rielles l'une de l'autre. Pour obtenir les équations de Fresnel, 
il faudrait transformer seulement le vecteur u par une fonction 
vectorielle symétrique. Au point de vue mathématique, toutes 
ces hypothèses sont acceptables; au point de vue physique, 
elles correspondent à des façons très différentes de concevoii" 
comment s'établit l'anisotropie d'un milieu. En l'absence de 
faits d'expérience, nous ne poursuivrons pas plus loin cette 
discussion. Nous conclurons simplement du cas particulier que 
nous venons d'examiner que l'étude de l'électricité pourrait 
nous éclairer sur l'exactitude des assertions de Fresnel. En 
effet, si la vibration lumineuse est représentée par la force 
magnétique ou électrique, ce sont les formules de Maxwell ou 
de Lamé qu'il faut adopter; si, au contraire, cette vibration 
correspond au déplacement magnétique ou électrique, les hypo- 
thèses de Fresnel doivent être admises. 
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